Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,729,182 T→C 6.0% A32A (GCT→GCC lhr → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,729,1820TC 6.0% 254.3 / 4.5 83A32A (GCT→GCClhrputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base T (37/41);  new base C (2/3);  total (39/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

TATAGTGAGGCTATTCCACGCATCCTGCATGATATTCACGGGGAATAGCGTTAATGGCAGATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCTCCGGGAAAACGCTGGCAGCATTTCTCTACGCCCTCGATCGGCTCTTCCGCGAAGGCGGCGAAGATACCCGCGAGGCGCA  >  NC_000913/1729034‑1729328
                                                                                                                                                    |                                                                                                                                                  
tataGTGAGGCTATTCCACGCATCCTGCATGATATTCACGGGGAATAGCGTTAATGGCAGATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCt                                                                                                                                                    >  1:426172/1‑149 (MQ=255)
 ataGTGAGGCTATTCCACGCATCCTGCATGATATTCACGGGGAATAGCGTTAATGGCAGATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTg                                                                                                                                                   >  4:30729/1‑149 (MQ=255)
        ggCTATTCCACGCATCCTGCATGATATTCACGGGGAATAGCGTTAATGGCAGATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGcc                                                                                                                                            <  4:90518/149‑1 (MQ=255)
           tATTCCACGCATACTGAATGAAATTAACGGGGAATAGCGTTAATGCCATATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGCTTTCTTCGCGCCGTTAAACAGCCGACCGCTGTCCcgccgc                                                                                                                                          >  1:226609/1‑148 (MQ=255)
             ttCCACGCATCCTGCATGATATTCACGGGGAATAGCGTTAATGGCAGATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCaaa                                                                                                                                       >  1:85169/1‑149 (MQ=255)
                cACGCATCCTGCATGATATTCCCCCGGAAGAGCGTTAACGGCAGATAATCCAGACCCTTCAGCGCTCCCCCCCGACGTGTTTCCCCCGGCGACCCGCGACCGGCCTCACCGCGCCTCTAAACCCCCGCCCGCCGCCCAGCCGCAAACCt                                                                                                                                    <  6:306371/149‑1 (MQ=255)
                cACGCATCCTGCATGATATTCACGGGGAATAGCGTTAATGGCAGATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCt                                                                                                                                    <  6:137904/149‑1 (MQ=255)
                cACGCATCCTGCATGATATTCACGGGGAATAGCGTTAATGGCAGATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCt                                                                                                                                    <  2:426172/149‑1 (MQ=255)
                 aCGCATCCTGCATGATATTCCCGGGGAAGAGCGTTAAAGGCAGATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTATTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAATCAGCCGCCCGCTGTCCAGCCGCAAACCTg                                                                                                                                   >  6:177278/1‑149 (MQ=255)
                 aCGCATCCTGCATGATATTCATGGGGAATAGCGTTAATGGCAGATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCa                                                                                                                                               <  5:157126/137‑1 (MQ=255)
                 aCGCATCCTGCATGATATTCATGGGGAATAGCGTTAATGGCAGATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCa                                                                                                                                               >  6:157126/1‑137 (MQ=255)
                 aCGCATCCTGCATGATATTCACGGGGAATAGCGTTAATGGCAGATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTg                                                                                                                                   <  5:322945/149‑1 (MQ=255)
                 aCGCATCCTGCATGATATTCACGGGGAATAGCGTTAATGGCAGATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCACGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTg                                                                                                                                   >  3:278498/1‑149 (MQ=255)
                 aCGCATCCTGCATGATATTCACGGGGAATAGCGTTAATGGCAGATAATCCAGACCCTTCACCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTg                                                                                                                                   <  3:97909/149‑1 (MQ=255)
                  cGCATCCTGCATGATATTCACGGGGAATAGCGTTAATGGCAGATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTgg                                                                                                                                  <  5:319712/149‑1 (MQ=255)
                   gCATCCTGCATGATATTCACGGGGAATAGCGTTAATGGCAGATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGc                                                                                                                                 <  2:30425/149‑1 (MQ=255)
                    cATCCTGCATGATATTCACGGGGAATAGCGTTAATGGCAGATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCa                                                                                                                                <  7:393015/149‑1 (MQ=255)
                         tGCATGATATTCACGGGGAATAGCGTTAATGGCAGATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTgg                                                                                                                           >  3:251345/1‑149 (MQ=255)
                           cATGATATTCACGGGGAATAGCGTTAATGGCAGATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTgg                                                                                                                           <  5:107608/147‑1 (MQ=255)
                               atatTAACTGGGAATAGCTTTAATGGCAGATTATCAAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACGGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGgcgc                                                                                                                     <  6:100755/149‑1 (MQ=255)
                                  ttCACGGGGAATAGCGTTAATGGCAGATACTCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCCGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCCTGTGGCGGCgcgaa                                                                                                                  <  6:134712/149‑1 (MQ=255)
                                  ttCACGGGGAATAGCGTTAATGGCAGATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGcgcg                                                                                                                    >  8:333452/1‑147 (MQ=255)
                                  ttCACGGGGAATAGCGTTAATGGCAGATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGcgcg                                                                                                                    >  4:136602/1‑147 (MQ=255)
                                     aCGGGGAATAGCGTTAATGGCAGATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCgcgaagcg                                                                                                               <  3:358434/149‑1 (MQ=255)
                                     aCGGGGAATAGCGTTAATGGCAGATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCgcgaagc                                                                                                                >  1:254423/1‑148 (MQ=255)
                                        gggAATAGCGGTAATGGCAGATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAAc                                                                                                            <  5:311139/149‑1 (MQ=255)
                                         ggAATAGCGTTAATGGCAGATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCAGGTGGCGGCGCGACGCGAACa                                                                                                           >  3:385128/1‑149 (MQ=255)
                                               gCGTTAATGGCAGATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCt                                                                                                     <  8:264582/149‑1 (MQ=255)
                                                cGTTAATGGCAGATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTg                                                                                                    >  4:16885/1‑149 (MQ=255)
                                                 gTTAATGGCAGATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTgg                                                                                                   >  3:292054/1‑149 (MQ=255)
                                                 gTTAATGGCAGATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTgg                                                                                                   <  1:254011/149‑1 (MQ=255)
                                                    aaTGGCAGATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTg                                                                                                 <  7:333452/148‑1 (MQ=255)
                                                     aTGGCAGATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCAACGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGa                                                                                                <  8:154241/148‑1 (MQ=255)
                                                      tGGCAGATAATCCAGCCCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGAtt                                                                                              >  8:304054/1‑149 (MQ=255)
                                                          aGATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACACCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCac                                                                                          >  3:225553/1‑149 (MQ=255)
                                                               aTCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCaccgac                                                                                      <  8:46829/148‑1 (MQ=255)
                                                               aTCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCaccgac                                                                                      <  7:8282/148‑1 (MQ=255)
                                                                  cAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGc                                                                                  >  5:273925/1‑149 (MQ=255)
                                                                   aGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCt                                                                                 >  1:431098/1‑149 (MQ=255)
                                                                       ccTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCacc                                                                                         >  6:176158/1‑137 (MQ=255)
                                                                       ccTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCacc                                                                                         <  5:176158/137‑1 (MQ=255)
                                                                         ttCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCTCCggg                                                                            >  5:97533/1‑148 (MQ=255)
                                                                                cTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCCGTCCAGCCGCAAACCt                                                                                                                                    <  6:18783/85‑1 (MQ=255)
                                                                                cTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCCGTCCAGCCGCAAACCt                                                                                                                                    <  6:13016/85‑1 (MQ=255)
                                                                                cTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCCGTCCAGCCGCAAACCt                                                                                                                                    >  5:13016/1‑85 (MQ=255)
                                                                                cTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCCGTCCAGCCGCAAACCt                                                                                                                                    >  5:18783/1‑85 (MQ=255)
                                                                                   cTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCTCCGGGAAAACGCTgg                                                                  <  2:431098/148‑1 (MQ=255)
                                                                                   cTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCTCCGGGAAAACGCTg                                                                   <  4:292054/147‑1 (MQ=255)
                                                                                    tGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCTCCGGGAAAACGCTgg                                                                  >  1:218835/1‑147 (MQ=255)
                                                                                      ccGGACGTGTTTTCACCGGCGACCCGCGACTGGTGTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCTCCGGGAAAACGCTGgcagc                                                              <  2:172703/149‑1 (MQ=255)
                                                                                        ggACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCTCCGGGAAAACGCTgg                                                                  <  2:139123/143‑1 (MQ=255)
                                                                                        ggACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCTCCGGGAAAACGCTgg                                                                  >  1:139123/1‑143 (MQ=255)
                                                                                        ggACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCTCCGGGAAAACGCTGGCAGCAt                                                            >  8:195457/1‑149 (MQ=255)
                                                                                         gACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCTCCGGGAAAACGCTGGCAGCAtt                                                           >  8:267458/1‑149 (MQ=255)
                                                                                            gtgtTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCTCCGGGAAAACGCTGGCAGCATTtct                                                        >  3:225980/1‑149 (MQ=255)
                                                                                              gtTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCTCCGGGAAAACGCTGGCAGCATTtctc                                                       >  5:21214/1‑148 (MQ=255)
                                                                                               ttttCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGACGACCGCTGTCCAGCCGCAAAACTGGCTTGTTGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCTGCTCCGGGAAAACGCTGGCAGCATTTCTCTa                                                     <  2:226609/149‑1 (MQ=255)
                                                                                                   cACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGc                                                                                                                                 <  5:429809/69‑1 (MQ=255)
                                                                                                   cACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGc                                                                                                                                 >  6:429809/1‑69 (MQ=255)
                                                                                                             ccGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGCCCGGCTCCGGGAAAACGCTGGCAGCATTTCTCTACGCCCTCGATCGGc                                       >  1:339309/1‑149 (MQ=255)
                                                                                                                       tttCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCTCCGGGAAAACGCTGGCAGCATTTCTCTACGCCCTCGATCGGCTCTTCCGCGa                             <  6:97533/149‑1 (MQ=255)
                                                                                                                         tcttcGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCTCCGGGAAAACGCTGGCAGCATTTCTCTACGCCCTCGATCGGCTCTTCCGCGaa                            <  4:225553/148‑1 (MQ=255)
                                                                                                                          cttcGCGCCTTTAAACATCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGAGAACATGCGCTGGAGATAGCACCGCCCGGCTCCGGGAAAACGCTGGCAGCATTTCACGACGCACTCGATCGGCTCTACCGCGAAgg                          >  1:82838/1‑149 (MQ=255)
                                                                                                                           ttcGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCTCCGGGAAAACGCTGGCAGCATTTCTCTACGCCCTTGATCGGCTCTTCCGCGAAggc                         >  4:172615/1‑149 (MQ=255)
                                                                                                                              gcgcCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATgc                                                                                                        >  4:76426/1‑67 (MQ=255)
                                                                                                                              gcgcCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATgc                                                                                                        <  3:76426/67‑1 (MQ=255)
                                                                                                                                gcCTTTAACCAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCTCCGGGAAAACGCTGGCAGCATTTCTCTACGCCCTCGATCGGCTCTTCCGCGAAGGCGGCGa                    <  1:150381/149‑1 (MQ=255)
                                                                                                                                 ccTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCGCCGGGAAAACGCTGGCAGCATTTCTCTACGCCCTCGATCGGCTCTTCCGCGAAGGCGGCGaa                   <  7:195457/149‑1 (MQ=255)
                                                                                                                                  cTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCTCCGGGAAAACGCTGGCAGCATTTCTCTACGCCCTCGATCGGCTCTTCCGCGAAGgcggcg                     >  7:75882/1‑146 (MQ=255)
                                                                                                                                   tttAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCTCCGGGAAAACGCTGGCAGCATTTCTCTACGCCCTCGATCGGCTCTTCCGCGAAGgcggcg                     >  1:440102/1‑145 (MQ=255)
                                                                                                                                    ttAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCTCCGGGAAAACGCTGGCAGCATTTCTCTACGCCCTCGATCGGCTCTTCCGCGAAGGCGGCGAAGa                 >  7:123490/1‑148 (MQ=255)
                                                                                                                                     tAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCTCCGGGAAAACGCTGGCAGCATTTCTCTACGCCCTCGATCGGCTCTTCCGCGAAGGCGGCGAAGa                 <  4:225980/147‑1 (MQ=255)
                                                                                                                                      aaaCAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCTCCGGGAAAACGCTGGCAGCATTTCTCTACGCCCTCGATCGGCTCTTCCGCGAAGGCGGCGAAGATa               <  4:202242/148‑1 (MQ=255)
                                                                                                                                       aaCAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCTCCGGGAAAACGCTGGCAGCATTTCTCTACGCCCTCGATCGGCTCTTCCGCGAAGGCGGCGAAGATAcc             <  3:136602/149‑1 (MQ=255)
                                                                                                                                       aaCAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCTCCGGGAAAACGCTGGCAGCATTTCTCTACGCCCTCGATCGGCTCTTCCGCGAAGGCGGCGAAGATAcc             >  5:62036/1‑149 (MQ=255)
                                                                                                                                        aCAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCTCCGGGAAAACGCTGGCAGCATTTCTCTACGCCCTCGATCGGCTCTTCCGCGaa                            >  4:295035/1‑133 (MQ=255)
                                                                                                                                        aCAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCTCCGGGAAAACGCTGGCAGCATTTCTCTACGCCCTCGATCGGCTCTTCCGCGaa                            <  3:295035/133‑1 (MQ=255)
                                                                                                                                            ccGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCTCCGGGAAAACGCTGGCAGCATTTCTCTACGCCCTCGATCGGCTCTTCCGCGAAGGCGGCGAAGATACCCGCGa        <  7:88254/149‑1 (MQ=255)
                                                                                                                                             cGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCTCCGGGAAAACGCTGGCAGCATTTCTCTACGCCCTCGATCGGCTCTTCCGCGAAGgcggcg                     <  7:353735/135‑1 (MQ=255)
                                                                                                                                             cGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCTCCGGGAAAACGCTGGCAGCATTTCTCTACGCCCTCGATCGGCTCTTCCGCGAAGgcggcg                     >  8:353735/1‑135 (MQ=255)
                                                                                                                                             cGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCTCCGGGAAAACGCTGGCAGCATTTCTCTACGCCCTCGATCGGCTCTTCCGCGAAGGCGGCGAAGATACCCGCGAg       <  7:267458/149‑1 (MQ=255)
                                                                                                                                               aCCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCTCCGGGAAAACGCTGGCAGCATTTCTCTACGCCCTCGATCGGCTCTTCCGCGAAGGCGGCGAAGATACCCGCGAGgc     <  2:339309/149‑1 (MQ=255)
                                                                                                                                               aCCGCTGTCCAGCCGCAAACCTGGCATGGGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCGCCGGGAAAACGCTGGCAGCATTTCTCTACGCCCTCGATCGGCTCTTCCGCGAAGGCGGCGAAGATACCCGCGAGgc     <  3:16885/149‑1 (MQ=255)
                                                                                                                                                   cTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGCCCGGCTCCGGGAAAACGCTGGCAGCATTTCTCTACGCCCTCGATCGGCTCTTCCGCGAAGGCGGCGAAGATACCCGCGAGGCGCa  >  4:57795/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TATAGTGAGGCTATTCCACGCATCCTGCATGATATTCACGGGGAATAGCGTTAATGGCAGATAATCCAGACCCTTCATCGCTCCTGCCGGACGTGTTTTCACCGGCGACCCGCGACTGGTTTCTTCGCGCCTTTAAACAGCCGACCGCTGTCCAGCCGCAAACCTGGCATGTGGCGGCGCGAAGCGAACATGCGCTGGTGATTGCACCGACCGGCTCCGGGAAAACGCTGGCAGCATTTCTCTACGCCCTCGATCGGCTCTTCCGCGAAGGCGGCGAAGATACCCGCGAGGCGCA  >  NC_000913/1729034‑1729328

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: