Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,704,308 C→T 9.9% intergenic (‑245/+27) rnc ← / ← lepB RNase III/leader peptidase (signal peptidase I)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,704,3080CT 9.9% 451.5 / 30.9 151intergenic (‑245/+27)rnc/lepBRNase III/leader peptidase (signal peptidase I)
Reads supporting (aligned to +/- strand):  ref base C (64/72);  new base T (7/8);  total (71/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01

Only 100 of 151 total aligned reads displayed.

GTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAGCCC  >  NC_000913/2704161‑2704455
                                                                                                                                                   |                                                                                                                                                   
gTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACg                                                                                                                                                    >  2:231210/1‑149 (MQ=255)
  atatGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGc                                                                                                                                                  <  4:785341/149‑1 (MQ=255)
    atGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGa                                                                                                                                                <  1:1254053/149‑1 (MQ=255)
      ggggACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACa                                                                                                                                              <  4:786529/149‑1 (MQ=255)
      ggggACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACa                                                                                                                                              >  3:239439/1‑149 (MQ=255)
       gggACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGCTGTCGTTAGCCACGGGAGATTTATCGCATAAATAATTCACGTTGTCGCCATAACGGCGACaa                                                                                                                                             >  3:1390201/1‑149 (MQ=255)
       gggACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACaa                                                                                                                                             <  4:937738/149‑1 (MQ=255)
        ggACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACAAc                                                                                                                                            <  3:873904/149‑1 (MQ=255)
        ggACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAAc                                                                                                                                            <  3:1532506/149‑1 (MQ=255)
          aCCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACAACGt                                                                                                                                          <  2:788272/149‑1 (MQ=255)
            cAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGa                                                                                                                                        <  3:1299081/149‑1 (MQ=255)
             aaCGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGaa                                                                                                                                       >  4:489031/1‑149 (MQ=255)
             aaCGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTCCCAACCGTGAACTTCGGGGTAATATTCTATACACCACGACGCGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGCCGCCATAACGGCGACAACGTGaa                                                                                                                                       >  2:111049/1‑149 (MQ=255)
              aCGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAAc                                                                                                                                      <  4:685756/149‑1 (MQ=255)
                gCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGa                                                                                                                                                >  3:1078177/1‑137 (MQ=255)
                gCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGa                                                                                                                                                <  4:1078177/137‑1 (MQ=255)
                  tGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAg                                                                                                                                  >  2:1368556/1‑149 (MQ=255)
                    tttCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGAt                                                                                                                                <  4:239439/149‑1 (MQ=255)
                    tttCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGAt                                                                                                                                <  1:265876/149‑1 (MQ=255)
                       cGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAAGAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGc                                                                                                                             <  1:1368556/149‑1 (MQ=255)
                        gtgtGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCt                                                                                                                            >  3:101133/1‑149 (MQ=255)
                        gtgtGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCt                                                                                                                            <  1:568618/149‑1 (MQ=255)
                         tgtgCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTa                                                                                                                           >  2:234614/1‑149 (MQ=255)
                          gtgCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTAt                                                                                                                          >  3:1325290/1‑149 (MQ=255)
                             ccGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTaa                                                                                                                       >  3:828063/1‑149 (MQ=255)
                                   cAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGAc                                                                                                                                               <  2:673664/119‑1 (MQ=255)
                                   cAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGAc                                                                                                                                               >  1:673664/1‑119 (MQ=255)
                                    aaCCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATgc                                                                                                                >  4:139896/1‑149 (MQ=255)
                                     aCCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTaa                                                                                                                       >  4:8770/1‑141 (MQ=255)
                                     aCCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTaa                                                                                                                       <  3:8770/141‑1 (MQ=255)
                                     aCCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATgcc                                                                                                               >  1:390067/1‑149 (MQ=255)
                                       cTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGTCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATgccgc                                                                                                             <  2:635905/149‑1 (MQ=255)
                                       cTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCCCGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATgccgc                                                                                                             <  1:1041470/149‑1 (MQ=255)
                                       cTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGAt                                                                                                                                <  4:680021/130‑1 (MQ=255)
                                        tGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATgccgcc                                                                                                            >  1:259571/1‑149 (MQ=255)
                                          aggttGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCaa                                                                                                          <  4:874038/145‑1 (MQ=255)
                                           aGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATg                                                                                                                     <  2:366848/137‑1 (MQ=255)
                                           aGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATg                                                                                                                     >  1:366848/1‑137 (MQ=255)
                                               gCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGa                                                                                                                                 <  4:1004019/121‑1 (MQ=255)
                                               gCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGa                                                                                                                                 >  4:218662/1‑121 (MQ=255)
                                               gCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGa                                                                                                     <  2:504124/149‑1 (MQ=255)
                                               gCCAACCTTAAACTCCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGa                                                                                                     >  2:791241/1‑149 (MQ=255)
                                                 cAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGTATGCCGCCAATGCGACt                                                                                                   >  3:330451/1‑149 (MQ=255)
                                                    ccTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTaa                                                                                                <  1:1521991/149‑1 (MQ=255)
                                                     cTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAg                                                                                               >  3:1136815/1‑149 (MQ=255)
                                                           cTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTGGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGa                                                                                         <  2:1511248/149‑1 (MQ=255)
                                                            ttCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATgc                                                                                                                <  4:1268077/125‑1 (MQ=255)
                                                            ttCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGAc                                                                                        <  1:905429/149‑1 (MQ=255)
                                                                 gggAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGa                                                                                         <  4:240522/143‑1 (MQ=255)
                                                                 gggAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTc                                                                                   <  2:173848/149‑1 (MQ=255)
                                                                  ggAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCg                                                                                  <  2:221048/149‑1 (MQ=255)
                                                                   gAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCgg                                                                                 >  3:121616/1‑149 (MQ=255)
                                                                     atatTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGcc                                                                               <  3:1408521/149‑1 (MQ=255)
                                                                      tatTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCa                                                                              <  1:171469/149‑1 (MQ=255)
                                                                        ttCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATg                                                                                                                 <  4:1426982/112‑1 (MQ=255)
                                                                        ttCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAg                                                                                               >  3:686988/1‑130 (MQ=255)
                                                                              acacAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCctt                                                                      >  3:736204/1‑149 (MQ=255)
                                                                                acaACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAAc                                                                                                                                      >  1:1453244/1‑83 (MQ=255)
                                                                                acaACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAAc                                                                                                                                      <  2:1453244/83‑1 (MQ=255)
                                                                                acaACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCcttct                                                                    >  2:541407/1‑149 (MQ=255)
                                                                                  aaCGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTg                                                                  >  1:899527/1‑149 (MQ=255)
                                                                                   acgacgGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGc                                                                 <  2:119889/149‑1 (MQ=255)
                                                                                   acgacgGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGc                                                                 <  1:1191261/149‑1 (MQ=255)
                                                                                   acgacgGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGc                                                                 <  4:736204/149‑1 (MQ=255)
                                                                                    cgacgGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCt                                                                <  3:625581/149‑1 (MQ=255)
                                                                                     gacgGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGAc                                                                                        >  1:221908/1‑124 (MQ=255)
                                                                                      acgGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTa                                                              >  4:1476562/1‑149 (MQ=255)
                                                                                       cgGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTAt                                                             <  4:224501/149‑1 (MQ=255)
                                                                                        gggggATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATc                                                            <  1:683103/149‑1 (MQ=255)
                                                                                           ggATGTCGTTAGCCACGGGAGATTTATCTCATAAATAAGTCACGTTGTGGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTAGCCTTCTTGCTTATCTaa                                                         <  4:121616/149‑1 (MQ=255)
                                                                                              tGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCt                                                      >  1:1101004/1‑149 (MQ=255)
                                                                                                  gTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATc                                                  <  2:259571/149‑1 (MQ=255)
                                                                                                     aGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAg                                               >  3:1247220/1‑149 (MQ=255)
                                                                                                       ccACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGAt                                             <  4:97145/149‑1 (MQ=255)
                                                                                                        cACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATgg                                                                                                                    <  2:352948/77‑1 (MQ=255)
                                                                                                        cACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATgg                                                                                                                    >  1:352948/1‑77 (MQ=255)
                                                                                                              agaTTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGtt                                      <  3:229641/149‑1 (MQ=255)
                                                                                                                aTTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGc                                    <  2:655286/149‑1 (MQ=255)
                                                                                                                     tctcATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCccgac                               <  1:541407/149‑1 (MQ=255)
                                                                                                                     tctcATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCccgac                               <  2:1282413/149‑1 (MQ=255)
                                                                                                                     tctcATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCccgac                               <  2:839656/149‑1 (MQ=255)
                                                                                                                     tctcATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCccgac                               <  3:553970/149‑1 (MQ=255)
                                                                                                                       tcATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAAt                                                                                                                      <  4:1528706/60‑1 (MQ=255)
                                                                                                                       tcATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAAt                                                                                                                      >  3:1528706/1‑60 (MQ=255)
                                                                                                                       tcATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCccgaccg                             >  3:798464/1‑149 (MQ=255)
                                                                                                                          taaataaTTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCgaccgacc                          <  4:456872/149‑1 (MQ=255)
                                                                                                                           aaataaTTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCa                         >  4:1186918/1‑149 (MQ=255)
                                                                                                                                atTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTc                    <  4:101133/149‑1 (MQ=255)
                                                                                                                                 ttCACGTTGTCGCCGTTATGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCg                   >  4:733515/1‑149 (MQ=255)
                                                                                                                                   cACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGt                                       <  2:390703/127‑1 (MQ=255)
                                                                                                                                   cACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGt                                       >  1:390703/1‑127 (MQ=255)
                                                                                                                                    aCGTTGTCGCCATACCGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCtt                >  3:111834/1‑149 (MQ=255)
                                                                                                                                    aCGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCtt                >  2:1118098/1‑149 (MQ=255)
                                                                                                                                      gTTGTCGCCGTTATGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTcc              >  2:272207/1‑149 (MQ=255)
                                                                                                                                      gTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTcc              >  2:850711/1‑149 (MQ=255)
                                                                                                                                       ttGTCGCCGTTATGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCcgaccgac                           <  2:1239041/135‑1 (MQ=255)
                                                                                                                                       ttGTCGCCGTTATGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCcgaccgac                           >  1:1239041/1‑135 (MQ=255)
                                                                                                                                                tAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAGc    >  3:1030295/1‑149 (MQ=255)
                                                                                                                                                 aaCGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAGcc   <  3:355249/149‑1 (MQ=255)
                                                                                                                                                  aCGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAGccc  <  4:984045/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAGCCC  >  NC_000913/2704161‑2704455

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: