Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,959,924 A→T 5.1% intergenic (‑111/+88) ppiC ← / ← yifN peptidyl‑prolyl cis‑trans isomerase C (rotamase C)/PemK toxin family pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,959,9240AT 5.1% 281.8 / 2.7 104intergenic (‑111/+88)ppiC/yifNpeptidyl‑prolyl cis‑trans isomerase C (rotamase C)/PemK toxin family pseudogene
Reads supporting (aligned to +/- strand):  ref base A (41/57);  new base T (2/4);  total (43/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 105 total aligned reads displayed.

TTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGCAA  >  NC_000913/3959776‑3960070
                                                                                                                                                    |                                                                                                                                                  
tttACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGaa                                                                                                                                                    >  1:421893/1‑149 (MQ=255)
   aCAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAAt                                                                                                                                                 >  1:570762/1‑149 (MQ=255)
   aCAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAAt                                                                                                                                                 >  2:1277835/1‑149 (MQ=255)
   aCAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAAt                                                                                                                                                 >  1:577251/1‑149 (MQ=255)
   aCAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAAt                                                                                                                                                 >  2:756406/1‑149 (MQ=255)
   aCAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAAt                                                                                                                                                 >  1:1351620/1‑149 (MQ=255)
        gctctGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGGGAATAATTGTCAGGCGGCTATTTTAGCATGAGCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCggg                                                                                                                                            <  1:252612/145‑1 (MQ=255)
         atatGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGc                                                                                                                                           <  2:41178/149‑1 (MQ=255)
             gCAGTGCTGCTGCTGTTTTTGCCCTGATCGTGCCTTGCGTGAATAATTGTCCGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGATTTTCCGGGCtttt                                                                                                                                       <  1:1322704/149‑3 (MQ=255)
             gCAGGGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGATTTTCCGGGCtttt                                                                                                                                       <  2:323291/149‑3 (MQ=255)
               aGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCa                                                                                                                                     >  1:205666/1‑149 (MQ=255)
                 tgctgctgcTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAgg                                                                                                                                   <  1:1569187/149‑1 (MQ=255)
                   ctgctgctgGTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTc                                                                                                                                 <  2:114945/149‑1 (MQ=255)
                     gctgctgTTTTTGCCATGATCGTGCCTTGGGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTccc                                                                                                                               >  2:64057/1‑149 (MQ=255)
                       tgctgTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTc                                                                                                                             <  1:620118/149‑1 (MQ=255)
                       tgctgTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTc                                                                                                                             >  1:818016/1‑149 (MQ=255)
                          tgTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAg                                                                                                                          <  1:813423/149‑1 (MQ=255)
                            tttttGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAggg                                                                                                                        <  1:1224053/149‑1 (MQ=255)
                               ttGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGgag                                                                                                                     <  2:897303/149‑1 (MQ=255)
                               ttGCCATGATCGAGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGgag                                                                                                                     <  2:85532/149‑1 (MQ=255)
                                   cATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGaaaa                                                                                                                 >  1:1445526/1‑149 (MQ=255)
                                   cATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGaaaa                                                                                                                 >  2:859454/1‑149 (MQ=255)
                                      gATCGTGCCATGAGTGAAAACTTGCCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACtt                                                                                                              <  1:209038/149‑1 (MQ=255)
                                       aTCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGGCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTa                                                                                                             <  2:592803/149‑1 (MQ=255)
                                        tCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                                                                            >  1:17769/1‑149 (MQ=255)
                                         cGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAgg                                                                                                           <  1:725020/149‑1 (MQ=255)
                                         cGTGCCTTGAGTGAATAATTCTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTACTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAgg                                                                                                           <  1:32575/149‑1 (MQ=255)
                                           tcccTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACTTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGc                                                                                                         <  2:198930/147‑1 (MQ=255)
                                           ggCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGc                                                                                                         <  2:1312712/148‑1 (MQ=255)
                                             ccTTGAGTGAATAATTGTCAGGGGGCGATTTTAGCATGATCCGCAGTTGCTAACGTAATTGTCAGTCTGGGAGTGGCTGGTACATATATGTCAAAAGCCCGTAAAATACGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCct                                                                                                       <  2:822619/149‑1 (MQ=255)
                                             ccTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCct                                                                                                       >  2:976756/1‑149 (MQ=255)
                                              cTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCctc                                                                                                      >  1:1512643/1‑149 (MQ=255)
                                               ttGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCACAGGTCCCTCAAGGGGAGAAAACTTAGGGCctct                                                                                                     <  1:515735/149‑1 (MQ=255)
                                                 gAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTCTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGCAAAATCCGGGCTTCGCAGGTCCCACAAGGGGAGAAAACTTAGGGCCTCTAt                                                                                                   >  2:259887/1‑149 (MQ=255)
                                                   gTGAAAAATTGTCAGGCGGCGATATTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGGGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCCAGGGGAGAAAACTTAGGGCCTCTATgg                                                                                                 <  2:401313/149‑1 (MQ=255)
                                                     gAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGc                                                                                               <  1:687454/149‑1 (MQ=255)
                                                     gAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGc                                                                                               <  1:344799/149‑1 (MQ=255)
                                                        taatTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTAc                                                                                            <  1:601273/149‑1 (MQ=255)
                                                            tGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGTTCCCTCAAGGGTAGAAAACTTAGGGCCTCTATGGGCTACAAAg                                                                                        <  2:602002/149‑1 (MQ=255)
                                                               cAGGCGGCGTTTTTAGCATGATCCGCAGATGAAAACGTAAGTGTAAGTCAGGGTGTGTCTGGTACATATAGGTCAAAAGCACGGATAATCCGTTCTTAGCAGGTCTCTCAAGGGGTGAAAACTTAGGGCCTCTATGGGCTACTAAGTCa                                                                                     <  2:1429310/149‑1 (MQ=255)
                                                                      cgATTTTAGCATGATCCGCAGATGCTACCGTAATTGTCAGTCAGGGAGGGGCCGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCgg                                                                              <  1:634937/149‑1 (MQ=255)
                                                                      cgATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGATTTTCCGGGCtttt                                                                                                                                       <  1:612007/92‑3 (MQ=255)
                                                                      cgATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGATTTTCCGGGCtttt                                                                                                                                       >  2:612007/1‑90 (MQ=255)
                                                                       gATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCggg                                                                             >  2:827462/1‑149 (MQ=255)
                                                                        aTTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCgggg                                                                            >  1:462577/1‑149 (MQ=255)
                                                                          tttAGCTTGATCCGAAGATGCTAACGTAATTGTCAGGCAGGGAGTGGCTGGTACATATATGTCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCGCTCAAGGTGAGAAAACTTAGGGGCTCTATGGGCTACAAAGTCATTATCGGGGAt                                                                          <  2:863736/149‑1 (MQ=255)
                                                                          tttAGCATGATCCGCAGATGCTAACGTACTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAACGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGag                                                                          <  2:953934/149‑2 (MQ=255)
                                                                          tttAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACATAGGGCCTCTATGGGCTACAAAGTCATTAACGGGGAt                                                                          >  2:137739/1‑149 (MQ=255)
                                                                             aGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGt                                                                       >  2:1256366/1‑149 (MQ=255)
                                                                              gCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAACACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGtt                                                                      >  1:1514229/1‑149 (MQ=255)
                                                                                aTGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttat                                                                    <  2:770000/149‑1 (MQ=255)
                                                                                   aTCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatttt                                                                 >  2:351530/1‑149 (MQ=255)
                                                                                    tcagCAGATGCTAACGCAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAACCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGCATGGttatttta                                                                >  1:1150843/4‑149 (MQ=255)
                                                                                    tCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatttta                                                                <  2:456805/149‑1 (MQ=255)
                                                                                    tCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatttta                                                                >  1:632328/1‑149 (MQ=255)
                                                                                      cGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAACTCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttatt                                                              <  2:1001328/149‑1 (MQ=255)
                                                                                      cGCAGATGCTAACGGAATTGTGCGTCAGGGCGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGGCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAACTCATTATCGGGGATGGttattttatt                                                              <  2:555033/149‑1 (MQ=255)
                                                                                       gCATATCCTAAAGTAAATGTCACTCAGGGAGTGGATGGTACATATATGGCAAAAGCACTGAAAATCCGGTCTTCGCAGTTCCCTCAAGGGGAGAAAACTTAGGTCCTCTATTGGCTACAAAGTCATTATCGGGGATTGTTATTTTATTg                                                             <  2:371489/149‑1 (MQ=255)
                                                                                       gCAGATGCTAACGTAATTTTGAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCATGTCCCGCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTg                                                             <  2:1382834/149‑1 (MQ=255)
                                                                                       gCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTg                                                             >  1:1081100/1‑149 (MQ=255)
                                                                                       gCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCGCAAGGGGCGAAAACTTAGGGCCTCGAGGGGCTACAACGTCATGATCGGGGATGGTTATTTAATTg                                                             >  2:873474/1‑149 (MQ=255)
                                                                                        cAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGt                                                            <  2:1512643/149‑1 (MQ=255)
                                                                                         agactcTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTc                                                           <  1:1056650/144‑1 (MQ=255)
                                                                                           aTGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAATGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCaa                                                         <  1:859454/149‑1 (MQ=255)
                                                                                             gCTAACGTAATTGTCAGTCACGGAGTGGCTGGTACATATAAGCCAAAAGCCCGCAAAATCACGGCTTCCCAGGTCCCTCTAGGGGAGAAAACTTAGGGCCTCTATCGGCGACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAAtt                                                       >  2:547892/1‑149 (MQ=255)
                                                                                               tAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTa                                                     <  2:455440/149‑1 (MQ=255)
                                                                                                accGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCAACGCAGGTCCCTCAAGTGGAGAAAACTTAGGGCCACTATGGGCTACAAAGGCATTATCGGGGATGTTTATTGTATTGTAAATTTAt                                                    >  2:924934/3‑149 (MQ=255)
                                                                                                 aCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATg                                                   >  1:559364/1‑149 (MQ=255)
                                                                                                  cGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGTAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGa                                                  <  2:507489/149‑1 (MQ=255)
                                                                                                  cGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGa                                                  <  2:205666/149‑1 (MQ=255)
                                                                                                         gtcagtcaGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGatat                                           <  2:421893/149‑1 (MQ=255)
                                                                                                            agtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGTAAAATCCGGGCTTCACAGGTCCCTCAAGGGAAGAAAACTTAGTTCCTCTATGGGCTAAAAACTCATTATCGGGGAGGCGTATTTTATTGACAATTTATGAGCGATATgtc                                        >  2:890951/1‑147 (MQ=255)
                                                                                                             gtcagGGAGTGGGTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCa                                       >  2:446718/1‑149 (MQ=255)
                                                                                                             gtcagGGAGTGGCTGGTACATATATGGCAAAACCCCGGAAAATCATGGCTTCGCAGGTCACTCAACGGCAGAAAACTTAGGGCCTCTATGGCATACAAAATCAGTATCGTGGTTGGATATTTTATTGTCAATTTATGAGCGATATGGCa                                       >  1:1437060/1‑149 (MQ=255)
                                                                                                              tcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGATTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGGTATTTTATTGTCAATTTATGAGCGATATGGCat                                      <  2:1202995/149‑1 (MQ=255)
                                                                                                               cagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCatt                                     <  2:58304/149‑1 (MQ=255)
                                                                                                                agGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCattg                                    >  1:540632/1‑149 (MQ=255)
                                                                                                                  ggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCattgat                                  <  2:570762/149‑1 (MQ=255)
                                                                                                                       ggCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACt                             <  1:756406/149‑1 (MQ=255)
                                                                                                                        gCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACtt                            <  2:1351620/149‑1 (MQ=255)
                                                                                                                        gCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACtt                            >  1:43054/1‑149 (MQ=255)
                                                                                                                        gCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCGCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACtt                            <  2:175042/149‑1 (MQ=255)
                                                                                                                         cTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTtg                           <  1:351530/149‑1 (MQ=255)
                                                                                                                               cATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCACTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATTTGGCATTGATTGACTTGATGAcc                     >  2:295315/1‑149 (MQ=255)
                                                                                                                                atatatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTtgatga                       >  1:1208378/1‑146 (MQ=255)
                                                                                                                                atatatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTtgatga                       <  2:1208378/146‑1 (MQ=255)
                                                                                                                                  atatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCCTTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGcc                  >  1:434056/1‑149 (MQ=255)
                                                                                                                                   tatGGCAAACGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGccc                 <  2:896910/149‑1 (MQ=255)
                                                                                                                                       gCAACAGCCCGGAAAATCCGGGCTTCTCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACCAAGTCATAATCGGGGATGGTTATTTTATTGTCAATTTATGCGCGATATGGCATTGATTGACTTGATGACCGCCCTCTg             <  1:1459135/149‑1 (MQ=255)
                                                                                                                                       gCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCATGGGGAGAATACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGTTGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTg             >  1:1047883/1‑149 (MQ=255)
                                                                                                                                       gCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTg             >  1:605757/1‑149 (MQ=255)
                                                                                                                                        aaaaaGCCCGGAAAATCGGGGCTTCGCAGGTCCCTGAAGGGGAGAAAACTTAGGGCGTCTATGGGCTACAAAGGCATTATCGGGGATGGTTATTTTTTTGTGAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGt            <  2:120821/148‑1 (MQ=255)
                                                                                                                                         aaaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCGATGGGCTACAATGTCATTATCGGGGTGGGGTATTTTATTGGCAATTTATGAGGGATTTGGGATTGATTGGCTTGATGCCCGGGCTCGGTa           >  2:1056554/1‑149 (MQ=255)
                                                                                                                                         aaaaGCCCGGAAAATCCGGGCTTCACAGGTCCCTCAAGGCGCGAACACTTAGGGCCGGTATGGGCTACAAAGTAATTATCGGTGATGGTTTTTGTATTGTCAATATATGAGCTATATGGCATTGATTGACTTGATGACCGCCCTCTGTa           >  1:1197038/1‑149 (MQ=255)
                                                                                                                                          aaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGTGAGAAAACTTAGCGCCTCTATGTGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTAAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTAt          <  2:818016/149‑1 (MQ=255)
                                                                                                                                            aGCCCGGATTTTCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGt        >  2:525652/1‑149 (MQ=255)
                                                                                                                                              cccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCt      <  1:1277835/149‑1 (MQ=255)
                                                                                                                                              cccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCt      <  1:1256366/149‑1 (MQ=255)
                                                                                                                                              cccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCt      <  1:750446/149‑1 (MQ=255)
                                                                                                                                                  gAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGcaa  <  1:256952/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGCAA  >  NC_000913/3959776‑3960070

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