Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,173,481 C→A 6.4% intergenic (+25/+23) ygiN → / ← parE quinol monooxygenase/DNA topoisomerase IV, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,173,4810CA 6.4% 395.2 / 11.0 133intergenic (+25/+23)ygiN/parEquinol monooxygenase/DNA topoisomerase IV, subunit B
Reads supporting (aligned to +/- strand):  ref base C (66/58);  new base A (5/4);  total (71/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 133 total aligned reads displayed.

GCAGTGGGAAAGCATCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATC  >  NC_000913/3173333‑3173623
                                                                                                                                                    |                                                                                                                                              
gCAGTGGGAAAGCATCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCTGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCt                                                                                                                                               >  1:748743/1‑150 (MQ=255)
   gTGGGAAAGCATCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGAc                                                                                                                                            >  2:1088671/1‑150 (MQ=255)
      ggAAAGCATCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCgg                                                                                                                                         >  1:1023010/1‑150 (MQ=255)
       gAAAGCATCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCggg                                                                                                                                        >  1:1141663/1‑150 (MQ=255)
       gAAAGCATCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGc                                                                                                                                       >  2:1471839/1‑150 (MQ=255)
             aTCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGTCAGGATGGCCGGGCaaaca                                                                                                                                  <  1:474940/150‑4 (MQ=255)
              tCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGtt                                                                                                                                 <  1:812130/150‑1 (MQ=255)
               cgcgCATCTTGACGCGCATCTGCAAACCCCGCACAGGAAGGGGGATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATTAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTc                                                                                                                                <  2:270149/150‑1 (MQ=255)
               cgcgCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTc                                                                                                                                >  2:171056/1‑150 (MQ=255)
                gcgcATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGTCAGGATGGCCGGGCaaacaagg                                                                                                                               <  2:613594/150‑7 (MQ=255)
                gcgcATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGAAATGTTCt                                                                                                                               >  1:1536487/1‑150 (MQ=255)
                gcgcATCTTGAAGCGCATCTGCAAACCCCGAACATGAAGGCGTATAGCGAATCCGTAAAATGTGAAGTGATGGAGATGAATATCCGTATTCTGCAGCCAGGGATGTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCTGGCAATGTTCt                                                                                                                               >  1:1543174/1‑150 (MQ=255)
                gcgcATCTAGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGACGCCGTAAAAGGTGACGTGCGGGAGAGGAATCTCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCt                                                                                                                               <  1:111353/150‑1 (MQ=255)
                  gcATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCttt                                                                                                                             <  2:1626491/150‑1 (MQ=255)
                  gcATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGAGGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCttt                                                                                                                             <  1:414247/150‑1 (MQ=255)
                   cATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTc                                                                                                                            >  2:1617692/1‑150 (MQ=255)
                   cATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTc                                                                                                                            <  2:1380176/150‑1 (MQ=255)
                     tCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCt                                                                                                                          <  1:989396/150‑1 (MQ=255)
                     ccTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCt                                                                                                                          <  2:792868/149‑1 (MQ=255)
                      cTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCtt                                                                                                                         <  2:936880/150‑1 (MQ=255)
                      cTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCtt                                                                                                                         >  1:499213/1‑150 (MQ=255)
                         gAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTaa                                                                                                                      <  2:311935/150‑1 (MQ=255)
                           aGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAAc                                                                                                                    <  2:1554593/150‑1 (MQ=255)
                           aGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAAc                                                                                                                    >  1:1167962/1‑150 (MQ=255)
                           aGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAAc                                                                                                                    <  2:680032/150‑1 (MQ=255)
                            gcgcATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAAcc                                                                                                                   >  2:1170268/1‑150 (MQ=255)
                            gcgcATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAAcc                                                                                                                   <  1:568869/150‑1 (MQ=255)
                             cgcATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCt                                                                                                                  >  2:1377547/1‑150 (MQ=255)
                              gcATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTc                                                                                                                 >  1:1649265/1‑150 (MQ=255)
                                 tCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAAt                                                                                                              >  1:1342484/1‑150 (MQ=255)
                                       aaCCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGc                                                                                                        <  1:233038/150‑1 (MQ=255)
                                            cGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAATCCTCTATCTCCGCCATGt                                                                                                   >  1:632943/1‑150 (MQ=255)
                                            cGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGt                                                                                                   >  1:639117/1‑150 (MQ=255)
                                             gCACATGAAGGAGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGACTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTc                                                                                                  <  1:1170268/150‑1 (MQ=255)
                                                caTGAAGGCTTATAGCGAAGCCGTAATAGGTGACGGGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATGTCCGCCATGTCGcc                                                                                               <  2:128902/150‑1 (MQ=255)
                                                 aTGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATTTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCt                                                                                              <  2:1691073/150‑1 (MQ=255)
                                                    aaGGCGGATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTtc                                                                                           <  2:134957/150‑1 (MQ=255)
                                                      ggCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTtctc                                                                                         <  1:356736/150‑1 (MQ=255)
                                                       gCGTATAGCGAGGCCGTAAAAGGTGACGTGCTGGAGATGAATCTCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTctct                                                                                        <  1:1457268/150‑1 (MQ=255)
                                                       gCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTctct                                                                                        <  1:1448732/150‑1 (MQ=255)
                                                        cGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCtt                                                                                       >  2:750431/1‑150 (MQ=255)
                                                         gTATAGCGAAGCCGTAACAGGTGACGTGCGGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTGGCCTTTCTCTTg                                                                                      <  2:1242098/150‑1 (MQ=255)
                                                         gTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCGTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTg                                                                                      <  1:327541/150‑1 (MQ=255)
                                                         gTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTg                                                                                      >  1:1005458/1‑150 (MQ=255)
                                                         gTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTg                                                                                      <  2:967065/150‑1 (MQ=255)
                                                                    ccGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTgc                                                                           <  1:487282/150‑1 (MQ=255)
                                                                     cGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTgct                                                                          >  1:261715/1‑149 (MQ=255)
                                                                       tAAAATGTGACGTGCTGGAGATGAATATCCGTATTCAGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTgcggc                                                                        <  1:1070868/150‑1 (MQ=255)
                                                                       tAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTgcggc                                                                        >  1:925141/1‑150 (MQ=255)
                                                                        aaaaGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTgcggcg                                                                       >  1:208394/1‑150 (MQ=255)
                                                                            ggTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCt                                                                   <  2:932373/150‑1 (MQ=255)
                                                                             gTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCtt                                                                  >  1:522971/1‑150 (MQ=255)
                                                                                   tGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGc                                                            >  1:748131/1‑150 (MQ=255)
                                                                                     cTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGt                                                          >  2:1468469/1‑150 (MQ=255)
                                                                                       ggAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGttt                                                        >  2:1671395/1‑150 (MQ=255)
                                                                                       ggAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGttt                                                        >  2:1454227/1‑150 (MQ=255)
                                                                                       ggAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCATGCCTTGTTTGCCCGGCCATCCTGAGCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGttt                                                        <  1:1542231/150‑1 (MQ=255)
                                                                                       ggAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATACTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGttt                                                        <  1:1535493/150‑1 (MQ=255)
                                                                                         agaTGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCAGGCAATCCGGACAGGCCAATGTTCTTTCCTTTAAAACTCAATCTCCGCAATGTCGACTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCCTttct                                                      >  1:823175/1‑150 (MQ=255)
                                                                                          gaTGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTttctt                                                     >  2:735753/1‑150 (MQ=255)
                                                                                           aTGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGCTCTTCCGAGCGTttcttc                                                    >  1:1674541/1‑150 (MQ=255)
                                                                                            tGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCg                                                   <  2:110800/150‑1 (MQ=255)
                                                                                            tGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCg                                                   >  2:1573491/1‑150 (MQ=255)
                                                                                            tGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCg                                                   <  2:440960/150‑1 (MQ=255)
                                                                                            tGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCg                                                   >  1:1647465/1‑150 (MQ=255)
                                                                                            tGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCATGACCGGTCAATGTTCTTTCCTTTAAACCTCAATATCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCg                                                   >  1:1172870/1‑150 (MQ=255)
                                                                                             gAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGAGGCGATCTTCCGAGCGTTTCTTCGc                                                  >  2:870768/1‑150 (MQ=255)
                                                                                              aaTATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGcc                                                 <  1:1400628/150‑1 (MQ=255)
                                                                                              aaTATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGcc                                                 <  1:1124579/150‑1 (MQ=255)
                                                                                                tatCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCAGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCcag                                               <  1:708254/150‑1 (MQ=255)
                                                                                                 atCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCcagc                                              >  2:1027284/1‑150 (MQ=255)
                                                                                                  tCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGACATGTCTCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCcagcc                                             >  2:496133/1‑149 (MQ=255)
                                                                                                   ccGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCcagcag                                            <  1:780381/150‑1 (MQ=255)
                                                                                                   ccGTATTCTGCAGCCAGGGATTTAATCATGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCAATGTCGCCTTTCTATTGCAACCAGTTTCGGAGATCTTCCTAGCGTTTCTTCGCcagcaa                                            >  1:316204/1‑149 (MQ=255)
                                                                                                     gTATTCTGCAGCCAGGGATTTACTCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCCgcagca                                          >  2:1069440/1‑150 (MQ=255)
                                                                                                     gTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCACCCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAgcggca                                          >  1:928328/1‑150 (MQ=255)
                                                                                                              atccAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTGCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCcatca                                 <  2:856876/148‑1 (MQ=255)
                                                                                                                ccAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCatcatc                               >  1:1322510/1‑150 (MQ=255)
                                                                                                                ccAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCatcatc                               <  1:1457601/150‑1 (MQ=255)
                                                                                                                 cAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATcg                              >  2:1562150/1‑150 (MQ=255)
                                                                                                                  aGGGATTTTATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATcgc                             >  1:557861/1‑150 (MQ=255)
                                                                                                                   gggATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATcgcg                            >  2:1640231/1‑150 (MQ=255)
                                                                                                                   gggATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATcgcg                            >  2:211790/1‑150 (MQ=255)
                                                                                                                   gggATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATcgcg                            >  2:321432/1‑150 (MQ=255)
                                                                                                                      aTTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCa                         >  1:1013925/1‑150 (MQ=255)
                                                                                                                          aaTCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTAc                     >  2:1243908/1‑150 (MQ=255)
                                                                                                                            tCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGc                   <  2:499213/150‑1 (MQ=255)
                                                                                                                              cTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTg                 >  1:939380/1‑150 (MQ=255)
                                                                                                                               tGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGa                >  1:1551728/1‑150 (MQ=255)
                                                                                                                                gCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGAt               >  2:682/1‑150 (MQ=255)
                                                                                                                                     gTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCt          <  2:261715/150‑1 (MQ=255)
                                                                                                                                     cattGCCCGGTCAGGATGGCCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCt          >  2:874135/3‑150 (MQ=255)
                                                                                                                                        tGCCCGGTCAGGATGGCCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTACTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTATCATCGCCAGCAGCATATCCAACATCGCGTCAGTACGCTGATCGTCTtca       >  1:608835/1‑150 (MQ=255)
                                                                                                                                        tGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTtca       <  2:1023010/150‑1 (MQ=255)
                                                                                                                                         gCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTtcat      >  1:1618171/1‑150 (MQ=255)
                                                                                                                                         gCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTtcat      >  2:1138577/1‑150 (MQ=255)
                                                                                                                                         gCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCCtcaa      >  2:125822/1‑149 (MQ=255)
                                                                                                                                          cccGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTtcatc     >  1:428584/1‑150 (MQ=255)
                                                                                                                                           ccGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTtcatca    >  2:492935/1‑150 (MQ=255)
                                                                                                                                             ggCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCatcatc  >  1:346347/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
GCAGTGGGAAAGCATCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATC  >  NC_000913/3173333‑3173623

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