Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,270,182 C→T 6.6% intergenic (+580/+34) 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,1820CT 6.6% 210.8 / 3.7 89intergenic (+580/+34)yhaC/rnpBuncharacterized protein YhaC/RNase P catalytic RNA component
Reads supporting (aligned to +/- strand):  ref base C (57/24);  new base T (3/3);  total (60/29)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAAC  >  NC_000913/3270045‑3270331
                                                                                                                                         |                                                                                                                                                       
ttcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGAGTCGTGGAGAGTCAAGCATGGTGTGGCGCCGACAAAAGCAAAAACccgcgg                                                                                                                                             >  2:606135/1‑148 (MQ=11)
        tCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTAAAGACTCGTAGAAAGTCAATCAAGATAATGCACCAACAAAAACAAAAACCCGCCGAAGCgcgg                                                                                                                                   >  2:1955070/1‑147 (MQ=11)
                 ccACAAGCTCAAACCCGCCGTATCGGTTTATTACGTAAAGCAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCTTCGGCGGGtttttgcttt                                                                                                                          <  1:1665954/150‑6 (MQ=11)
                   ttacgtaaaaaCCCGCTTCGGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGTGTATAGAGTATTGTAAAGTCAATCAAAATAAAACACCAACAAAAACAAAAAACCGCCGAACCGGGTTTTTACATaaa                                                                                                                        >  2:254255/7‑150 (MQ=2)
                      cgtaaaaaCCCGCTTCGGCGGGTTTTTGCGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGTGTATAGCGTAGTGGACAGTCATTCAAATTTATGCACCAACAAAAACAAAAACCCGCCGAAGAGGGTTTTTACGTAAATCa                                                                                                                     >  2:1195395/4‑150 (MQ=11)
                      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)
                                                       aCCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGTCTAGAGAGTCGAGGAAAGTCCATCAACAAAAACCACCATCAAAAACAAAAACCCGCATAAGAACGTTGTTAAGTTAATCAGGTGAAAATTACCGATAATACGGGTTgtgtatt                                                                                    >  2:1801759/1‑147 (MQ=2)
                                                           ggTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGCAGCGGGTTTTTACGTAAATCAGGTCCAACTGACCGATAAGCCGGGTTATGTCGTGGAc                                                                                >  1:1349333/1‑150 (MQ=255)
                                                           ggTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTACAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGaa                                                                                >  1:1462377/1‑149 (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                                                                       <  1:698248/142‑1 (MQ=255)
                                                                            aGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtcattca                                                                       >  2:698248/1‑142 (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)
                                                                                           tgTATAGCGTCGTTGAGAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATTAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTg                                                <  1:536601/149‑1 (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                            >  1:25306/1‑150 (MQ=255)
                                                                                                               cattcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCggg                            >  2:1633959/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:260859/1‑150 (MQ=255)
                                                                                                                  tcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTc                         >  2:29245/1‑150 (MQ=255)
                                                                                                                  tcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTc                         >  1:69267/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                        >  1:975644/1‑150 (MQ=255)
                                                                                                                   catCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCa                        >  2:1052570/1‑150 (MQ=255)
                                                                                                                    atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGACGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAg                       <  1:1724131/150‑1 (MQ=255)
                                                                                                                    atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAg                       <  1:2060703/150‑1 (MQ=255)
                                                                                                                    atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCAGAAGCCGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGACAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAg                       >  2:277217/1‑150 (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)
                                                                                                                            gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCggg                                                                                             >  2:949222/1‑72 (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)
                                                                                                                                     aaaTCACAAACCCGCCGCACCGGGTTTTTACGTACATCAGGTCAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTAt      <  1:1128023/150‑1 (MQ=39)
                                                                                                                                     aaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGCGCCGTACCTTAt      <  1:2391706/150‑1 (MQ=255)
                                                                                                                                      aaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATg     >  1:580185/1‑150 (MQ=255)
                                                                                                                                      aaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATg     <  1:2249788/150‑1 (MQ=255)
                                                                                                                                        gCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGaa   >  2:2333280/1‑150 (MQ=255)
                                                                                                                                         cAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAAc  >  2:723114/1‑150 (MQ=255)
                                                                                                                                         |                                                                                                                                                       
TTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAAC  >  NC_000913/3270045‑3270331

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: