Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,560,914 A→T 5.1% I180N (ATC→AAC)  glpG ← rhomboid protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,560,9140AT 5.1% 252.7 / 2.9 79I180N (ATC→AAC) glpGrhomboid protease
Reads supporting (aligned to +/- strand):  ref base A (40/35);  new base T (2/2);  total (42/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

TAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGCCAGAACTCAAATTTCAG  >  NC_000913/3560766‑3561062
                                                                                                                                                    |                                                                                                                                                    
tAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGAt                                                                                                                                                     <  1:1971892/150‑1 (MQ=255)
    ccACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGc                                                                                                                                                 >  2:1992057/1‑150 (MQ=255)
     cACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTATCAGCGCCCTGTTAAGCg                                                                                                                                                <  1:458672/150‑1 (MQ=255)
     cACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCg                                                                                                                                                <  2:124716/150‑1 (MQ=255)
       cGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTATCAGCGCCCTGTTAAGCggc                                                                                                                                              <  2:1159163/150‑3 (MQ=255)
          tGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATg                                                                                                                                           >  1:1337126/1‑150 (MQ=255)
          tGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATg                                                                                                                                           >  1:1598980/1‑150 (MQ=255)
                               ggATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCg                                                                                                                      >  2:865415/1‑150 (MQ=255)
                                         cGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGttttt                                                                                                            >  2:1083670/1‑150 (MQ=255)
                                           ccACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGCATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTcc                                                                                                          >  2:719878/1‑150 (MQ=255)
                                              cGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCcacc                                                                                                       <  2:1504311/150‑1 (MQ=255)
                                                tAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGtttt                                                                                                             >  2:925043/1‑142 (MQ=255)
                                                tAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGtttt                                                                                                             <  1:925043/142‑1 (MQ=255)
                                                  gCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCcaccgcac                                                                                                   >  2:736844/1‑150 (MQ=255)
                                                  gCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCcaccgcac                                                                                                   >  2:743700/1‑150 (MQ=255)
                                                     aGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCaccgcaccgc                                                                                                <  2:1537931/150‑1 (MQ=255)
                                                       aCGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCAccgccg                                                                                              >  2:1111529/1‑150 (MQ=255)
                                                         gTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCgag                                                                                            <  1:152432/150‑1 (MQ=255)
                                                          tAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCgaga                                                                                           >  1:1845845/1‑150 (MQ=255)
                                                              ccATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATacc                                                                                       <  1:836620/150‑1 (MQ=255)
                                                              ccATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATacc                                                                                       <  2:1845845/150‑1 (MQ=255)
                                                              ccATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGAACTCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATacc                                                                                       >  2:981148/1‑150 (MQ=255)
                                                               cATCAGCGCATACACCACGCTAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATacca                                                                                      >  1:253129/1‑150 (MQ=255)
                                                               cATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCcac                                                                                                        <  2:1677817/132‑1 (MQ=255)
                                                               cATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCcac                                                                                                        >  1:1677817/1‑132 (MQ=255)
                                                               cATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATaccc                                                                                      >  2:842974/1‑149 (MQ=255)
                                                               cATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATacca                                                                                      >  2:486572/1‑150 (MQ=255)
                                                                   aGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATaccaccac                                                                                  >  1:1016373/1‑150 (MQ=255)
                                                                    gcgcATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATaccaccacc                                                                                 >  1:1931994/1‑150 (MQ=255)
                                                                     cgcATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATAccaccacca                                                                                <  2:1016373/150‑1 (MQ=255)
                                                                      gcATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCag                                                                               <  1:776360/150‑1 (MQ=255)
                                                                      gcATACACCACGCCAGAAAGCCCGCCAAACAACGGCCAGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCag                                                                               >  2:1615244/1‑150 (MQ=255)
                                                                      gcATACACAACGCAAGAAAGCCCGCCAAACCACGGCCAGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCag                                                                               <  2:531125/150‑1 (MQ=255)
                                                                        aTACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTAAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCagag                                                                             >  2:1885028/1‑150 (MQ=255)
                                                                             ccacGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGt                                                                        >  2:119504/1‑150 (MQ=255)
                                                                             ccacGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGt                                                                        >  1:879181/1‑150 (MQ=255)
                                                                             ccacGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGt                                                                        <  1:1199975/150‑1 (MQ=255)
                                                                                cGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTaa                                                                     <  1:1992057/150‑1 (MQ=255)
                                                                                 gCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTaaa                                                                    >  1:1604268/1‑150 (MQ=255)
                                                                                   cAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAaga                                                                  >  2:1334944/1‑150 (MQ=255)
                                                                                    aGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAagag                                                                 <  1:1111529/150‑1 (MQ=255)
                                                                                     gAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAgg                                                                <  1:486572/150‑1 (MQ=255)
                                                                                        aGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGata                                                             <  1:635303/150‑1 (MQ=255)
                                                                                        aGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGata                                                             <  2:1428957/150‑1 (MQ=255)
                                                                                            cGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCa                                                         <  1:345244/150‑1 (MQ=255)
                                                                                             gCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCAt                                                        >  1:74751/1‑150 (MQ=255)
                                                                                               cAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCa                                                      >  1:1172270/1‑150 (MQ=255)
                                                                                                   ccACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGCTACCACCACCAGAGCAGGTTAAAGAGGATATGCCTCAGCga                                                  >  2:441302/1‑150 (MQ=255)
                                                                                                     aCGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCgaga                                                <  1:1908095/150‑1 (MQ=255)
                                                                                                        gCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGt                                             <  2:382753/150‑1 (MQ=255)
                                                                                                             cTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCAtt                                        <  1:736844/150‑1 (MQ=255)
                                                                                                             cTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCAtt                                        <  1:743700/150‑1 (MQ=255)
                                                                                                             cTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCAtt                                        >  2:1803348/1‑150 (MQ=255)
                                                                                                               gAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTaa                                      >  2:1707778/1‑150 (MQ=255)
                                                                                                                 aTTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTg                                            >  1:1923265/1‑142 (MQ=255)
                                                                                                                 aTTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTg                                            <  2:1923265/142‑1 (MQ=255)
                                                                                                                 aTTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAAcg                                    <  1:1620803/150‑1 (MQ=255)
                                                                                                                 aTTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAAcg                                    <  1:222509/150‑1 (MQ=255)
                                                                                                                     ttGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTg                                <  2:31635/150‑1 (MQ=255)
                                                                                                                     ttGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTg                                <  2:24308/150‑1 (MQ=255)
                                                                                                                         tgcACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTg                            >  2:93553/1‑150 (MQ=255)
                                                                                                                           cactTAGCCGCTTAACAGGGCGGTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCAGCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGGGTGGGTGaa                          <  2:1122367/150‑1 (MQ=255)
                                                                                                                           cacaTAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGaa                          >  1:1441976/1‑150 (MQ=255)
                                                                                                                              cattacGCTTATCAGCGCCCTGTTAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTa                       >  1:1958694/6‑150 (MQ=255)
                                                                                                                              aTAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTa                       >  1:279908/1‑150 (MQ=255)
                                                                                                                                 gCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACg                    <  1:719878/150‑1 (MQ=255)
                                                                                                                                  acGCTTATCAGCGCCCTGTTAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGc                   >  1:1619435/2‑150 (MQ=255)
                                                                                                                                    gCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGaa                                               >  2:272488/1‑120 (MQ=255)
                                                                                                                                    gCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGaa                                               <  1:272488/120‑1 (MQ=255)
                                                                                                                                    gCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTggg                              <  2:344605/137‑1 (MQ=255)
                                                                                                                                    gCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTggg                              >  1:344605/1‑137 (MQ=255)
                                                                                                                                    gCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGCCa                 <  1:1334944/150‑1 (MQ=255)
                                                                                                                                        aaCAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGCCAGAAc             >  2:1892252/1‑150 (MQ=255)
                                                                                                                                         aCAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGCCAGAACt            >  2:1844643/1‑150 (MQ=255)
                                                                                                                                          cAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGCCAGAACTc           >  1:712845/1‑150 (MQ=255)
                                                                                                                                          cAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGCCAGAACTc           >  1:1588384/1‑150 (MQ=255)
                                                                                                                                           aGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGCCAGAACTCa          >  2:1681902/1‑150 (MQ=255)
                                                                                                                                               cgcTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGCCAGAACTCAAAtt      <  2:1206197/150‑1 (MQ=255)
                                                                                                                                                   gATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGCCAGAACTCAAATTTCAg  <  1:981148/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGCCAGAACTCAAATTTCAG  >  NC_000913/3560766‑3561062

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