Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 304,088 G→T 7.2% T32N (ACT→AAT)  ykgJ ← putative zinc‑ or iron‑chelating domain‑containing protein YkgJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913304,0880GT 7.2% 254.4 / 5.8 76T32N (ACT→AAT) ykgJputative zinc‑ or iron‑chelating domain‑containing protein YkgJ
Reads supporting (aligned to +/- strand):  ref base G (39/31);  new base T (2/4);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATGTGGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTACATCGTATGTAAG  >  NC_000913/303940‑304237
                                                                                                                                                    |                                                                                                                                                     
gCATGTGGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGt                                                                                                                                                      >  1:7602/1‑150 (MQ=255)
 cATGTGGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCACGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTAACGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTa                                                                                                                                                     >  2:1266152/1‑150 (MQ=255)
   tgtgGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTAcc                                                                                                                                                   >  1:2499/1‑150 (MQ=255)
     tgGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGc                                                                                                                                                 >  1:1133603/1‑150 (MQ=255)
        aCGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGTACTATTCCCGCCAg                                                                                                                                              <  2:1409614/150‑1 (MQ=255)
         cgatcgatTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGTACTATTCCCGCCAgg                                                                                                                                             <  2:432379/150‑2 (MQ=255)
         cgatcgatTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCACTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGc                                                                                                                                             <  2:1108731/150‑1 (MQ=255)
           atcgatTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGTACTATTCCCGCCAggct                                                                                                                                           <  1:744879/150‑4 (MQ=255)
           atcgatTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCAt                                                                                                                                           >  1:282222/1‑150 (MQ=255)
           atcgatTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCAt                                                                                                                                           >  1:1404492/1‑150 (MQ=255)
                tttttATATTTCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACTTCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCAtcgtcg                                                                                                                                      <  2:583232/150‑1 (MQ=255)
                     atatatCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTc                                                                                                                                 <  1:1016994/150‑1 (MQ=255)
                     atatatCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTc                                                                                                                                 <  2:1425125/150‑1 (MQ=255)
                     aaaaaacGTGCAGCAGGCATTTTTGCCCGGGGTTCATGCAAGGGCAAAACATCGCGGGTATTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGCGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTc                                                                                                                                 >  1:1394052/7‑150 (MQ=255)
                           cGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATTGTACCGCCAGCATCGTCgg                                                                                                                                     <  2:1430616/140‑1 (MQ=255)
                           cGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATGGTACCGCCAGCATCGTCgg                                                                                                                                     >  1:1430616/1‑140 (MQ=255)
                             tGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGt                                                                                                                         >  2:345769/1‑150 (MQ=255)
                               cagcagGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTaa                                                                                                                       <  2:1191089/150‑1 (MQ=255)
                               cagcagGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTaa                                                                                                                       <  2:1186174/150‑1 (MQ=255)
                                  cagGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTaaaaa                                                                                                                    >  2:875076/1‑150 (MQ=255)
                                   agGCATTTTTGCCCGTGGTTCCTGCAAGGGCAATACTTCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAag                                                                                                                   <  2:104793/150‑1 (MQ=255)
                                   agGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATAGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAag                                                                                                                   >  2:797315/1‑150 (MQ=255)
                                      cATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAgaga                                                                                                                >  1:121683/1‑150 (MQ=255)
                                       aTTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGAc                                                                                                               <  1:93240/150‑1 (MQ=255)
                                             gCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGt                                                                                                                         <  1:1363072/134‑1 (MQ=255)
                                             gCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGt                                                                                                                         >  2:1363072/1‑134 (MQ=255)
                                                  gggTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATg                                                                                                    <  2:984382/150‑1 (MQ=255)
                                                  gggTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATg                                                                                                    >  1:256914/1‑150 (MQ=255)
                                                   ggTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGc                                                                                                   <  2:1275876/150‑1 (MQ=255)
                                                   ggTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGc                                                                                                   >  2:832844/1‑150 (MQ=255)
                                                     ttCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGcac                                                                                                 <  1:288963/150‑1 (MQ=255)
                                                       ccTGCAAGGGCAATACATCGGGGGTTTCTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCacaa                                                                                               >  2:1331332/1‑150 (MQ=255)
                                                         tGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCacaacc                                                                                             >  2:408985/1‑149 (MQ=255)
                                                          gCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAg                                                                                            <  2:761347/150‑1 (MQ=255)
                                                                         cGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTcatg                                                                             >  1:547959/1‑150 (MQ=255)
                                                                             ggTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTcatgcatg                                                                         <  1:258894/150‑1 (MQ=255)
                                                                               tttttCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGa                                                                       >  2:817430/1‑150 (MQ=255)
                                                                                   tCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCa                                                                   >  1:903327/1‑150 (MQ=255)
                                                                                    cTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCCCAACAGGCACCACACGTCATGCATGGATTCAg                                                                  >  2:1229760/1‑150 (MQ=255)
                                                                                    cTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAg                                                                  <  1:882548/150‑1 (MQ=255)
                                                                                         tGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGc                                                             <  2:1133603/150‑1 (MQ=255)
                                                                                                   cATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCa                                                   <  2:256914/150‑1 (MQ=255)
                                                                                                    aTGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCAc                                                  <  1:1331332/150‑1 (MQ=255)
                                                                                                     tGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCAcc                                                 <  2:136377/150‑1 (MQ=255)
                                                                                                       cATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAg                                               >  1:1087288/1‑150 (MQ=255)
                                                                                                             tGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGc                                         <  1:194143/150‑1 (MQ=255)
                                                                                                               aaaaGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTa                                       >  2:1287877/1‑150 (MQ=255)
                                                                                                                 aaGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACATGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTaaa                                     <  1:1266152/150‑1 (MQ=255)
                                                                                                                 aaGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTaaa                                     <  2:1404492/150‑1 (MQ=255)
                                                                                                                   ggggATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTaaaaa                                   >  1:1270208/1‑150 (MQ=255)
                                                                                                                   ggggATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTaaaaa                                   <  1:832844/150‑1 (MQ=255)
                                                                                                                     ggATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATc                                 >  2:340587/1‑150 (MQ=255)
                                                                                                                        tatTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTa                              <  2:1087288/150‑1 (MQ=255)
                                                                                                                          tttGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTaaa                            >  2:980789/1‑150 (MQ=255)
                                                                                                                           tttttaaGTGAGCCTGGCGGGATTATTAACGCCAGCATGGTCGGCTTCTGCCCAGTTAAAAGAGAGGCGTAAAAATGGACAACATGCACCACACGTCATGCATTGATTCAGATTGATCATAATTCACCAGTACAGCTAAAAATCGTAAag                           <  2:1270208/144‑1 (MQ=255)
                                                                                                                           ttGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAag                           <  1:727128/150‑1 (MQ=255)
                                                                                                                            tGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAaga                          >  1:857559/1‑150 (MQ=255)
                                                                                                                            tGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAaga                          >  1:1197474/1‑150 (MQ=255)
                                                                                                                             gTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAagag                         >  2:1014850/1‑150 (MQ=255)
                                                                                                                             gTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGGACAGCTAAAAATCGTAAagag                         >  1:925624/1‑150 (MQ=255)
                                                                                                                               tCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGaa                       >  1:516182/1‑150 (MQ=255)
                                                                                                                               tCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGaa                       <  2:121683/150‑1 (MQ=255)
                                                                                                                                cAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAAc                      <  1:875076/150‑1 (MQ=255)
                                                                                                                                 aGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTAGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATTCATGGATTAAGATTGCTCATAATTAACCAGTACAGCTAAAAATCGTAAAGAGAACa                     >  2:702174/1‑150 (MQ=255)
                                                                                                                                  gTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAg                    >  1:874131/1‑150 (MQ=255)
                                                                                                                                  gTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATCCAGAATGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAg                    >  2:1285203/1‑150 (MQ=255)
                                                                                                                                   tGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGc                   <  1:345769/150‑1 (MQ=255)
                                                                                                                                    gatgcTGGCGGTACTATTCCCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCg                  >  1:359966/5‑150 (MQ=255)
                                                                                                                                    gAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCg                  >  2:645134/1‑150 (MQ=255)
                                                                                                                                     aGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCgg                 >  2:430009/1‑150 (MQ=255)
                                                                                                                                       gcTGGCGGTACTATTCCCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTa               >  1:1374613/2‑150 (MQ=255)
                                                                                                                                          ggCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTACAt            >  2:170879/1‑150 (MQ=255)
                                                                                                                                             gggAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCCCCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTACATCGt         >  2:325787/1‑150 (MQ=255)
                                                                                                                                                aaTAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTACATCGTATg      <  2:2499/150‑1 (MQ=255)
                                                                                                                                                  tAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTACATCGTATGTa    <  2:320676/150‑1 (MQ=255)
                                                                                                                                                   aGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTACATCGTATGTaa   >  1:540880/1‑150 (MQ=255)
                                                                                                                                                    gTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTACATCGTATGTAAg  >  2:467364/1‑150 (MQ=255)
                                                                                                                                                    gTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTACATCGTATGTAAg  >  2:969061/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCATGTGGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTACATCGTATGTAAG  >  NC_000913/303940‑304237

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