Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,270,195 A→C 11.7% intergenic (+593/+21) yhaC → / ← rnpB uncharacterized protein YhaC/RNase P catalytic RNA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,270,1950AC11.7% 225.5 / 10.9 107intergenic (+593/+21)yhaC/rnpBuncharacterized protein YhaC/RNase P catalytic RNA component
Reads supporting (aligned to +/- strand):  ref base A (61/32);  new base C (9/4);  total (71/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 107 total aligned reads displayed.

TCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCC  >  NC_000913/3270053‑3270343
                                                                                                                                              |                                                                                                                                                      
tCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTAAAGACTCGTAGAAAGTCAATCAAGATAATGCACCAACAAAAACAAAAACCCGCCGAAGCgcgg                                                                                                                                               >  2:1955070/1‑147 (MQ=11)
         ccACAAGCTCAAACCCGCCGTATCGGTTTATTACGTAAAGCAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCTTCGGCGGGtttttgcttt                                                                                                                                      <  1:1665954/150‑6 (MQ=11)
           ttacgtaaaaaCCCGCTTCGGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGTGTATAGAGTATTGTAAAGTCAATCAAAATAAAACACCAACAAAAACAAAAAACCGCCGAACCGGGTTTTTACATaaa                                                                                                                                    >  2:254255/7‑150 (MQ=2)
              aGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTGTGCCCCTGCAAAAGCAAAAACCCGCCGGAGCGGGTTTTTAAGTAAAtga                                                                                                                                 >  1:2216969/1‑148 (MQ=255)
               gtaaaaaCCCGCTTCGGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCAGTCAGACCCAACAAAAACAAAAACCCGCAGAAGCGCCGTTTTACGTAAATCaa                                                                                                                                >  2:163578/3‑149 (MQ=11)
                cAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGCAAGTCAATCATAGTTCTGCGCCTACAAAAGCAAAAACCCGCAGAAGAGGGTTTTTTAGAAAATCacg                                                                                                                               >  2:1835586/1‑148 (MQ=11)
                cAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGAGAGTCATTCAACTTTCGGCACCTACAAAAGCAAAAACCCGCAGAAGCGGGTTTTTAAGTAAATCAgg                                                                                                                               >  2:933932/1‑150 (MQ=11)
                      cccGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGAAAGTCAATCAAAATTAGACACCAACAAAAACAAAAACCCCCCTAAGACCGGTATTTACGAAATCAAGAGaaaa                                                                                                                         >  2:591221/1‑149 (MQ=2)
                        cgccgAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCAATCATCTTTAAGCGCCAACAAAAGCAAAAACCCGCCGAAGAACGTTTTTAAGTAAATCAGGTGaaaatt                                                                                                                       >  2:1511106/1‑147 (MQ=11)
                                cGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGCCAGTCATTCATATTTGTGCCCCGGCAAAAGCAAAAACCCGCCGAAGAGGGGTTTTAAGGAAATGAGGGGAAAATTGCCGATaa                                                                                                               >  1:264766/1‑150 (MQ=11)
                                  ggTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGTGTATAGCGTCGTGGAAAGTCAATCATAATTTGGCGCCTACAAAAGCAAAAACCCGCGGAAGGGGGGTTTTTACTAAATGAGGGGAAAATTGCCGATGAGc                                                                                                             >  2:293752/1‑150 (MQ=2)
                                  ggTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTGGGGGAGAGGCATTCATGTTGGAGCGCCAGCAAAAGCAAAAACCCGCCGAAGGGCGTGGGTACGGTAATCAGGGGCAAATGACCGATAAGc                                                                                                             >  2:2298216/1‑150 (MQ=2)
                                      tttACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGCAGCGGGTTTTTACGTAAATCAGGTGTAACTGACCCATAAGCCggg                                                                                                         >  1:1145350/1‑150 (MQ=255)
                                        tACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGAGTCGTGGAAAGTCATTCATCATTGAGCCCCAACAAAAGCAAAAACCCGCCGAAGGGGGGTTTTACGTAAATCAGGAGAAAATGACCCATAATCCGGGtt                                                                                                       >  1:1572645/1‑150 (MQ=11)
                                            tAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCAATCATAATTAAGCACCAACAAAACCAAAAACCCGCCGCAGCGCGGTTTTACGAAAATCAGGGGAAAATTCCCGATAAGCCGGGTTCAGt                                                                                                   >  2:198240/1‑150 (MQ=14)
                                                   ggTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGCAGCGGGTTTTTACGTAAATCAGGTCCAACTGACCGATAAGCCGGGTTATGTCGTGGAc                                                                                            >  1:1349333/1‑150 (MQ=255)
                                                   ggTGAAACTGACCGATAAGCCGCGTTCTTTTGGGTATAGCGTCTTGGACAGTCAATCAACATGGCCCGCCAACAAAAACAAAAACCCGCAGAAGCGCCGTGATACGTAAATCAGGTGAAACTGACCAAGAATCCCGGTTCTGTCggggcc                                                                                            >  2:1429242/1‑148 (MQ=9)
                                                    gTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCAGCCCCTACAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAAGTAAATCAGTTGAAACTGACCGATAAGCCGGGTTATGTCGTgggga                                                                                           >  2:1173472/1‑147 (MQ=14)
                                                     tGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTAAATCATCATTCAGCCCCAACAAAAACAAAAACCCCCCGAAGACCCATATTAACTAAATCAGCTGAAACTGACCCATAAGCCCGGTAATGTAGTGGACac                                                                                          >  2:2142991/1‑149 (MQ=2)
                                                          cTGACCGTTTAGCCGCATTCTTTTGTGTTTAGCGGCGTGGACAGTCAGTCATCATTGTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtcatt                                                                                     <  1:1834597/150‑1 (MQ=17)
                                                            gACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGCAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGAAAGtcattca                                                                                   >  1:1655021/1‑150 (MQ=255)
                                                              ccGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGAGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATc                                                                                 >  2:1442392/1‑150 (MQ=255)
                                                              ccGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCGCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtcattgatg                                                                                 >  1:2409521/1‑149 (MQ=255)
                                                              ccGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCAACAAAAACAAAAACCCGCCGAATCCGCTTTTTACGTAAATCAGGTGAAACTGACCGATAAACCAGGTTATGTAGTTGACAGtcaatcatg                                                                                 >  2:2060703/1‑149 (MQ=14)
                                                                gATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTg                                                                                                                             <  1:1378405/104‑1 (MQ=35)
                                                                   aaGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGc                                                                            >  1:416002/1‑150 (MQ=255)
                                                                   aaGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCCGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGc                                                                            >  2:2064159/1‑150 (MQ=255)
                                                                    aGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtcattca                                                                                   >  2:698248/1‑142 (MQ=255)
                                                                    aGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtcattca                                                                                   <  1:698248/142‑1 (MQ=255)
                                                                     gCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGCAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCa                                                                          >  2:704427/1‑150 (MQ=255)
                                                                     gCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAACCTGCCCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGccc                                                                          >  2:1351257/1‑149 (MQ=255)
                                                                     gCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCa                                                                          >  1:797669/1‑150 (MQ=255)
                                                                       cGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCAGGGTTATGTAGTGGACAGTCATTCATCTAGGCCAGc                                                                        >  2:1725478/1‑150 (MQ=255)
                                                                           ttcttaTGGGTTTTGTGGCGTTGTCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATc                                                                    <  1:1325414/150‑1 (MQ=18)
                                                                            tctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGtt                                                                                                       >  2:1231991/1‑114 (MQ=18)
                                                                            tctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGtt                                                                                                       <  1:1231991/114‑1 (MQ=18)
                                                                             ctttTGGGTATAGCGTCGTGGACAGTCATCCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGc                                                                  >  1:1338122/1‑150 (MQ=255)
                                                                                   ggTATAGCGCCGTTGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAACAACCCTCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTg                                                            <  1:1172161/150‑1 (MQ=255)
                                                                                    gTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTgg                                                           >  1:1396238/1‑150 (MQ=255)
                                                                                            tctttgACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATTATGGGTTTTGTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAagca                                                   <  2:248335/145‑1 (MQ=18)
                                                                                             cgtGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAagcag                                                  >  1:1852632/1‑150 (MQ=255)
                                                                                               tGGACAGTCATTCATCTTTCTGCCCATCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGcc                                                <  1:254255/150‑1 (MQ=255)
                                                                                                ggACAGTCATTCATCTTTTTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATAGCTCACTGGCTCAAGCAGCCt                                               >  1:646036/1‑150 (MQ=255)
                                                                                                    aGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACt                                                                                                                        <  2:699066/73‑1 (MQ=34)
                                                                                                    aGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTAccc                                           >  2:1117627/1‑150 (MQ=255)
                                                                                                      tcattcatCTTTCATCCCCTCCAAAACCAAAAACCCGCCGAATCTGGTTTTTACGTAAATCAGGTGAACCTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCgg                                         <  1:295782/150‑1 (MQ=37)
                                                                                                       cattcatCTTTCTGCCCCTCCAAAAGCCACACGCCGCCGCAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCggg                                        <  1:1429242/150‑1 (MQ=37)
                                                                                                       cattcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCggg                                        >  2:1633959/1‑150 (MQ=255)
                                                                                                       cattcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCggg                                        >  1:25306/1‑150 (MQ=255)
                                                                                                        tctcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGt                                       <  1:1786430/148‑1 (MQ=255)
                                                                                                        attcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTg                                                                                                                       <  2:524936/70‑1 (MQ=35)
                                                                                                        attcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTCACCGATAAGCAGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGt                                       >  2:1547071/1‑150 (MQ=255)
                                                                                                          tcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTc                                     >  1:69267/1‑150 (MQ=255)
                                                                                                          tcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTc                                     >  1:260859/1‑150 (MQ=255)
                                                                                                          tcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTc                                     >  2:29245/1‑150 (MQ=255)
                                                                                                           catCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGt                                                                                         <  2:186145/97‑1 (MQ=255)
                                                                                                           catCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGt                                                                                         >  1:186145/1‑97 (MQ=255)
                                                                                                           catCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCa                                    >  2:1708003/1‑150 (MQ=255)
                                                                                                           catCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCa                                    >  2:1721710/1‑150 (MQ=255)
                                                                                                           catCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCa                                    >  2:1052570/1‑150 (MQ=255)
                                                                                                           catCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCa                                    >  1:975644/1‑150 (MQ=255)
                                                                                                            atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGACGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAg                                   <  1:1724131/150‑1 (MQ=255)
                                                                                                            atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAg                                   <  1:2060703/150‑1 (MQ=255)
                                                                                                             tCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGt                                  >  1:1357704/1‑150 (MQ=255)
                                                                                                             tCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGt                                  >  1:1366428/1‑150 (MQ=255)
                                                                                                              cTTTCTGACCATCGAAAAGCAAAAACCCGCAGAAGCTGGTTTTTACGTAAATCATGTGAATCTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTa                                 <  2:1655021/150‑1 (MQ=37)
                                                                                                                    gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAg                                                                         >  1:864614/1‑104 (MQ=255)
                                                                                                                    gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAg                                                                         <  2:864614/104‑1 (MQ=255)
                                                                                                                    gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCggg                                                                                                         <  1:949222/72‑1 (MQ=38)
                                                                                                                    gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGcc                           <  1:1719381/150‑1 (MQ=255)
                                                                                                                     ccccTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTTAATCTGACCGATAAGCGGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCg                          <  2:574412/150‑1 (MQ=255)
                                                                                                                      cccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGt                         >  1:2481593/1‑150 (MQ=255)
                                                                                                                      cccTCCAAAAGAAACAACCCGCCGAAGCGGGTTTTTACGTAAACCACGGGCTCCTGACCGCTAAGGCGGCTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGt                         <  1:313864/150‑1 (MQ=37)
                                                                                                                          ccAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCt                     >  2:1236138/1‑150 (MQ=255)
                                                                                                                           cAAAATCAAAACCCCCCCCCAGCCTTTTTTTACGTACAGCCCCTGAACCTGACCGCTCTGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCtt                    <  1:1045506/150‑1 (MQ=255)
                                                                                                                            aaaaccaaaaccCCGCCGAAGCGGGTTTTTAGGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTa                   <  1:277613/150‑1 (MQ=255)
                                                                                                                            aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTg                                                                                               <  1:849194/74‑1 (MQ=255)
                                                                                                                            aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTg                                                                                               >  2:849194/1‑74 (MQ=255)
                                                                                                                             tacgtaaaaaCCCGCTTCGGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTAt                  >  2:28419/6‑150 (MQ=38)
                                                                                                                             aaaTCACAAACCCGCCGCACCGGGTTTTTACGTACATCAGGTCAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTAt                  <  1:1128023/150‑1 (MQ=39)
                                                                                                                             aaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGCGCCGTACCTTAt                  <  1:2391706/150‑1 (MQ=255)
                                                                                                                              aaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATg                 >  1:580185/1‑150 (MQ=255)
                                                                                                                               cgtaaaaaCCCGCTTCGGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGa                >  2:2078909/4‑150 (MQ=39)
                                                                                                                               cgtaaaaaCCCGCTTCGGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGAAAACGCTCACAGGCTCAAGCAGCCTACCCGGGATCAGTACGGGCCGTACCGTAAGa                >  1:2276258/4‑150 (MQ=37)
                                                                                                                                gCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGaa               >  2:2333280/1‑150 (MQ=255)
                                                                                                                                ataaaaaCCCGCTTCGGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTCATGaa               >  1:406100/3‑150 (MQ=38)
                                                                                                                                 cAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAAc              >  2:723114/1‑150 (MQ=255)
                                                                                                                                  ttaaaCCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAAcc             <  1:1251549/148‑1 (MQ=255)
                                                                                                                                  aaaaaCCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAAcc             <  1:591221/150‑1 (MQ=255)
                                                                                                                                     aaCCCGCAGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCt          >  1:1657624/1‑150 (MQ=255)
                                                                                                                                       cccGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCt                                                                                                     <  1:2471036/57‑1 (MQ=255)
                                                                                                                                       cccGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTAt        >  2:1826806/1‑150 (MQ=255)
                                                                                                                                       cccGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTAt        >  2:1229164/1‑150 (MQ=255)
                                                                                                                                       cccGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGGGTTCt                                                                                                     >  2:2471036/1‑57 (MQ=255)
                                                                                                                                        ccgccgAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTAtt       <  1:1276780/150‑1 (MQ=255)
                                                                                                                                        ccgccgAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTAtt       >  1:894627/1‑150 (MQ=255)
                                                                                                                                         cgccgAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTAttt      >  1:1723576/1‑150 (MQ=255)
                                                                                                                                         cgccgAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTAttt      >  2:698794/1‑150 (MQ=255)
                                                                                                                                           ccgAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTgg    <  1:204298/150‑1 (MQ=255)
                                                                                                                                             gAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGcc  <  1:606135/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                      
TCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCC  >  NC_000913/3270053‑3270343

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: