Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,606,557 T→G 6.0% intergenic (‑156/+112) leuQ ← / ← rsmC tRNA‑Leu/16S rRNA m(2)G1207 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,606,5570TG 6.0% 288.6 / 6.2 93intergenic (‑156/+112)leuQ/rsmCtRNA‑Leu/16S rRNA m(2)G1207 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (38/49);  new base G (2/4);  total (40/53)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 113 total aligned reads displayed.

AATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGC  >  NC_000913/4606410‑4606703
                                                                                                                                                   |                                                                                                                                                  
aaTTCTTAAAAAGAATTGCTACCACGGAGGCGGATTCTAGTGTTTTTCAGCTTTTCGTCTATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcagat                                                                                                                                                    <  3:212441/148‑1 (MQ=255)
aaTTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcagat                                                                                                                                                    <  7:172794/148‑1 (MQ=255)
aaTTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcagat                                                                                                                                                    <  7:83328/148‑1 (MQ=255)
 aTTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcagatgc                                                                                                                                                  >  5:165515/1‑149 (MQ=255)
 aTTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcagatgc                                                                                                                                                  >  4:8624/1‑149 (MQ=255)
    cTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcagatgcca                                                                                                                                                <  3:426779/148‑1 (MQ=255)
     ttAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcagatgc                                                                                                                                                  >  5:445463/1‑145 (MQ=255)
     ttAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcagatgc                                                                                                                                                  >  2:82197/1‑145 (MQ=255)
     ttAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcagatgc                                                                                                                                                  >  2:318677/1‑145 (MQ=255)
      tAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcagatgccagat                                                                                                                                             >  2:352012/1‑149 (MQ=255)
       aaaaaGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcagatgccagatg                                                                                                                                            >  5:60499/1‑149 (MQ=255)
       aaaaaGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcagatgccagatg                                                                                                                                            >  2:207602/1‑149 (MQ=255)
       aaaaaGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcagatgccagat                                                                                                                                             <  5:478802/148‑1 (MQ=255)
        aaaaGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcagatgccagatgc                                                                                                                                           <  6:41195/149‑1 (MQ=255)
         aaaGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTTGTATGCATGCGGGGGcagatgccagatga                                                                                                                                           >  5:84612/1‑147 (MQ=255)
               ttGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTg                                                                                                                                    >  8:258436/1‑149 (MQ=255)
               ttGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTg                                                                                                                                    >  7:145879/1‑149 (MQ=255)
                   tACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGgcgc                                                                                                                                >  5:448588/1‑149 (MQ=255)
                     ttcaaggtgtatcgcgccattatgacccgccaggcgaagaaaggttaattatctcgccgaataccgtgtcggatgcggcgtgaacgccttatccgacctacaaaacctgcacgttagcccttcgtaGGCCAGATGCGACGCTGGCGCGt                                                                                                                              <  1:453375‑M1/23‑1 (MQ=255)
                     ccACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcagatgccag                                                                                                                                               >  7:360428/1‑132 (MQ=255)
                     ccACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcagatgccag                                                                                                                                               <  8:360428/132‑1 (MQ=255)
                      cACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTc                                                                                                                             <  1:384882/149‑1 (MQ=255)
                      cACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTc                                                                                                                             <  6:448588/149‑1 (MQ=255)
                             gcgcATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTg                                                                                                                      >  1:166611/1‑149 (MQ=255)
                                tcgcgccattatgacccgccaggcgaagaaaggttaattatctcgccgaataccgtgtcggatgcggcgtgaacgccttatccgacctacaaaacctgcacgttagcccttcgtaGGCCAGATGCGACGCTGGCGCGTCTTATCTGGc                                                                                                                    >  7:293270‑M1/116‑148 (MQ=255)
                                 aTTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGcc                                                                                                                   >  3:384051/1‑148 (MQ=255)
                                  ttCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATAAGACGCGCCAGCGTCGCATCTGGcatc                                                                                                                 <  3:8624/149‑4 (MQ=38)
                                  atcaactGGATTTCCGCGTTTCCTCAATAGATAATTATCGCCAGAGGTGTCATTGCTGGAAAAATGGCCCTCAGGAACCTAGCGAGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTa                                                                                                                 <  2:374770/143‑1 (MQ=255)
                                     tAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACg                                                                                                               >  2:252427/1‑148 (MQ=255)
                                          gTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAggg                                                                                                          <  1:252427/148‑1 (MQ=255)
                                          gTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAggg                                                                                                          <  7:388289/148‑1 (MQ=255)
                                          gTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAggg                                                                                                          <  1:459808/148‑1 (MQ=255)
                                          gTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAggg                                                                                                          <  3:360813/148‑1 (MQ=255)
                                          gTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGc                                                                                                         <  6:2481/149‑1 (MQ=255)
                                             ttCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTaa                                                                                                      >  4:287616/1‑149 (MQ=255)
                                               cAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAAc                                                                                                     <  1:162122/148‑1 (MQ=255)
                                               cAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAAc                                                                                                     <  8:145879/148‑1 (MQ=255)
                                               cAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAAGTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAAc                                                                                                     <  2:312915/148‑1 (MQ=255)
                                                aGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGt                                                                                                   >  8:241233/1‑149 (MQ=255)
                                                     ttCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAg                                                                                               <  5:446428/148‑1 (MQ=255)
                                                        gTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGttt                                                                                           <  5:176352/149‑1 (MQ=255)
                                                         tCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGtttt                                                                                          >  8:124090/1‑149 (MQ=255)
                                                           aaTAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTg                                                                                         <  6:445463/148‑1 (MQ=255)
                                                           aaTAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTg                                                                                         <  4:384051/148‑1 (MQ=255)
                                                           aaTAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGt                                                                                        <  5:190970/149‑1 (MQ=255)
                                                           aaTAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGt                                                                                        <  7:324198/149‑1 (MQ=255)
                                                             tAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAg                                                                                      >  5:64582/1‑149 (MQ=255)
                                                              aGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGa                                                                                                                                         >  6:134457/1‑97 (MQ=255)
                                                              aGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGa                                                                                                                                         <  5:134457/97‑1 (MQ=255)
                                                               gTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGt                                                                                    >  4:251905/1‑149 (MQ=255)
                                                               gTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGt                                                                                    >  6:282938/1‑149 (MQ=255)
                                                                ttAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTc                                                                                   >  8:216434/1‑149 (MQ=255)
                                                                ttAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTc                                                                                   >  2:137547/1‑149 (MQ=255)
                                                                  aaTTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCg                                                                                  <  6:60499/148‑1 (MQ=255)
                                                                       tCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAg                                                                            <  7:216434/149‑1 (MQ=255)
                                                                               ggTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCAc                                                                    >  1:29039/1‑149 (MQ=255)
                                                                                 tgtAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCAcgc                                                                  >  7:130867/1‑149 (MQ=255)
                                                                                  gtAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCAcgcc                                                                 >  4:169633/1‑149 (MQ=255)
                                                                                   gaaatGCTGGAAAAATCTCCATCAGGAAACAAGCGAGCAGGTTTGGTATGCATGCGGGGGCAGCTGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTACCGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCAcgcc                                                                 <  8:304252/144‑1 (MQ=255)
                                                                                     tttGCTGGAAAAATGTGCATCAGGAAACTAGCGTGCCGGTTTGGTATTCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCATGTTTTGTAGGTCGGATAAGGCGTTCAcgccgc                                                               <  5:127267/147‑1 (MQ=255)
                                                                                     aTTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTTGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCAcgccgc                                                               <  3:79098/148‑1 (MQ=255)
                                                                                     aTTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCAcgccgc                                                               <  6:64582/148‑1 (MQ=255)
                                                                                     aTTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCa                                                              >  2:443193/1‑149 (MQ=255)
                                                                                      ttGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCAt                                                             >  8:297726/1‑149 (MQ=255)
                                                                                        gCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATc                                                            <  4:113180/148‑1 (MQ=255)
                                                                                         cTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCg                                                          >  2:229296/1‑149 (MQ=255)
                                                                                          tGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCAcgccgc                                                               <  5:120567/143‑1 (MQ=255)
                                                                                          tGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGa                                                         >  1:222028/1‑149 (MQ=255)
                                                                                           ggAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGa                                                         >  7:414118/1‑148 (MQ=255)
                                                                                            gAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGaca                                                       <  7:278723/149‑1 (MQ=255)
                                                                                              aaaaTGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACg                                                     <  2:166611/149‑1 (MQ=255)
                                                                                                    tCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTAtt                                                <  1:82197/148‑1 (MQ=255)
                                                                                                    tCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTc                                               >  4:348264/1‑149 (MQ=255)
                                                                                                    tCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGGTCACGCCGCATCCGACACGGTATTc                                               <  1:182838/149‑1 (MQ=255)
                                                                                                      cATCAGGAAACTAGCGTGCAGGTTTGGTATGCAGGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCg                                              <  1:352012/148‑1 (MQ=255)
                                                                                                       aTCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGc                                            >  3:217310/1‑149 (MQ=255)
                                                                                                          aGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCga                                          <  2:222028/148‑1 (MQ=255)
                                                                                                           ggAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGc                                            >  8:381528/1‑145 (MQ=255)
                                                                                                               atccgacctacaaaacctgcacgttagcccttcgtaGGCCAGATAAGACGCGCCAGCGTCGCATCTGGCATACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAAt                                    <  2:462371‑M1/113‑1 (MQ=255)
                                                                                                                 tAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCt                                                                                                        >  5:81837/1‑79 (MQ=255)
                                                                                                                 tAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCt                                                                                                        <  6:81837/79‑1 (MQ=255)
                                                                                                                 tAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAAtt                                   >  4:163355/1‑148 (MQ=255)
                                                                                                                   gCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAAc                                >  8:408851/1‑149 (MQ=255)
                                                                                                                      tGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGATATTCGGCGAGATAATTAACCtt                             >  1:352/1‑149 (MQ=255)
                                                                                                                        acaaaacctgcacgttagcccttcgtaGGCCAGATAAGACGCGCCAGCGTCGCATCTGGCATACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTttc                           <  1:231992‑M1/122‑1 (MQ=255)
                                                                                                                           gTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTttcttc                        >  7:435783/1‑149 (MQ=255)
                                                                                                                           gTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTATTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTttcttc                        <  2:8484/149‑1 (MQ=255)
                                                                                                                               ggTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGcc                     <  4:217310/148‑1 (MQ=255)
                                                                                                                               ggTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGcc                     <  7:297726/148‑1 (MQ=255)
                                                                                                                               ctgcacgttagcccttcgtaGGCCAGATAAGACGCGCCAGCGTCGCATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCt                    <  2:453375‑M1/129‑1 (MQ=255)
                                                                                                                                gTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTg                   <  8:414118/149‑1 (MQ=255)
                                                                                                                                 tATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTgg                  >  8:137037/1‑149 (MQ=255)
                                                                                                                                   tgcatgcGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCg                <  8:130867/149‑1 (MQ=255)
                                                                                                                                     gttagcccttcgtaGGCCAGATAAGACGCGCCAGCGTCGCATCTGGCATACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCggg              >  1:462371‑M1/15‑149 (MQ=255)
                                                                                                                                     catgcGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCgg               <  8:435783/148‑1 (MQ=255)
                                                                                                                                     catgcGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCATCGCCTGGCgg               <  4:431367/148‑1 (MQ=255)
                                                                                                                                      atgcGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGt             <  1:207602/149‑1 (MQ=255)
                                                                                                                                       tgcGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCAGCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTc            <  1:137547/149‑1 (MQ=255)
                                                                                                                                        gcGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCa           >  6:215959/1‑149 (MQ=255)
                                                                                                                                                 gatgccagatgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGgcgc  <  6:396396/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGC  >  NC_000913/4606410‑4606703

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: