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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,560,9140AT 7.3% 268.5 / 8.3 96I180N (ATC→AAC) glpGrhomboid intramembrane serine protease
Reads supporting (aligned to +/- strand):  ref base A (46/43);  new base T (5/2);  total (51/45)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.89e-01

TAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGCCAGAACTCAAATTTCAGT  >  NC_000913/3560766‑3561063
                                                                                                                                                    |                                                                                                                                                     
tAACCCAAGCTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACCCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGAt                                                                                                                                                      <  1:246408/150‑1 (MQ=255)
 aaCCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATa                                                                                                                                                     <  1:1169104/150‑1 (MQ=255)
 aaCCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATa                                                                                                                                                     <  1:1503114/150‑1 (MQ=255)
 aaCCCACGTTGCAGGTAAATGCCACTTTGCGCATCGAGTTCGCCACGTAGCCAGAAGTAGCCCATCAGCGCATACCCCACGCCAGAAAGCCCGCCAAACCACGGACCGCTGAATTTTTGATCCACATAGACGCTTACCAGGGCGCTGATa                                                                                                                                                     >  2:431209/1‑150 (MQ=255)
    ccACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGc                                                                                                                                                  >  2:1184957/1‑150 (MQ=255)
       cGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACATAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTATCAGCGCCCTGTTAAGCggc                                                                                                                                               <  1:1029656/150‑3 (MQ=255)
            cAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGAc                                                                                                                                          >  2:1475423/1‑150 (MQ=255)
              ggTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACaa                                                                                                                                        >  1:1622559/1‑150 (MQ=255)
              ggTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACaa                                                                                                                                        >  1:1616083/1‑150 (MQ=255)
                  aaTGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATGAGCTCATACACCACGCCAGAAAGCCCGCCATTCCACGTCCCGCTTAATTTTTGCTTCACATAGGCGCTTAACAGGGCGCTGATAAGCGTAATGACATTTAg                                                                                                                                    <  1:458949/150‑1 (MQ=255)
                   aTGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGc                                                                                                                                   <  2:1037242/150‑1 (MQ=255)
                   aTGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGc                                                                                                                                   >  1:1515170/1‑150 (MQ=255)
                    tGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCt                                                                                                                                  <  2:1418914/150‑1 (MQ=255)
                      ccACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTa                                                                                                                                >  1:1627166/1‑150 (MQ=255)
                       cACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTAc                                                                                                                               <  2:223139/150‑1 (MQ=255)
                        aCTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTAcc                                                                                                                              >  2:1019837/1‑150 (MQ=255)
                         cTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCg                                                                                                                             <  1:726370/150‑1 (MQ=255)
                              cGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCACACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATACGCGTAATGACAATTAGCTTACCGCTAcc                                                                                                                        >  1:250561/1‑150 (MQ=255)
                              cGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTAcc                                                                                                                        <  1:725400/150‑1 (MQ=255)
                               ggATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCg                                                                                                                       >  1:673596/1‑150 (MQ=255)
                               ggATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAAATTTTGCTGCACAAAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCg                                                                                                                       <  1:596325/150‑1 (MQ=255)
                                gATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGa                                                                                                                      >  1:52483/1‑150 (MQ=255)
                                 aTCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAg                                                                                                                     >  1:1136271/1‑150 (MQ=255)
                                  tCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTAAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAgg                                                                                                                    <  2:746662/150‑1 (MQ=255)
                                    attgtCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCg                                                                                                                  <  2:1266937/146‑1 (MQ=255)
                                    attgtCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCAAGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCg                                                                                                                  <  2:872127/146‑1 (MQ=255)
                                      gTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGCCAATTAGCTTACCGCTACCGAGGCGtt                                                                                                                >  1:195288/1‑150 (MQ=255)
                                        tCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCTGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGtttt                                                                                                              <  1:1286278/150‑1 (MQ=255)
                                         cGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGttttt                                                                                                             <  1:481721/150‑1 (MQ=255)
                                           ccACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTcc                                                                                                           <  1:918357/150‑1 (MQ=255)
                                           ccACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTcc                                                                                                           <  1:825900/150‑1 (MQ=255)
                                           ccACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCGGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACAGCTACCGAGGCGTTTTTcc                                                                                                           >  2:1186481/1‑150 (MQ=255)
                                                    cAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGAAATGACAATTAGCTTACCGCTAACGAGGCGTTTTTCcaccgccccg                                                                                                  >  2:253018/1‑150 (MQ=255)
                                                        cGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCga                                                                                              <  1:1389200/150‑1 (MQ=255)
                                                           aGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGAt                                                                                           >  1:1290771/1‑150 (MQ=255)
                                                             cccATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAATAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATac                                                                                         >  2:54429/1‑150 (MQ=255)
                                                             cccATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATac                                                                                         >  2:1058866/1‑150 (MQ=255)
                                                              ccATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATacc                                                                                        >  2:839938/1‑150 (MQ=255)
                                                              ccATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATacc                                                                                        >  1:21860/1‑150 (MQ=255)
                                                               cATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATacca                                                                                       >  2:957923/1‑150 (MQ=255)
                                                                aTCAGCGCATACACCAAGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATaccac                                                                                      <  2:381316/150‑1 (MQ=255)
                                                                 tCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATaccacc                                                                                     <  2:1658180/150‑1 (MQ=255)
                                                                   aGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATaccaccac                                                                                   >  1:1630891/1‑150 (MQ=255)
                                                                    gcgcATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATaccaccacc                                                                                  <  2:52483/150‑1 (MQ=255)
                                                                         tACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGCGATACCACCACCAGAGc                                                                             >  1:1529404/1‑150 (MQ=255)
                                                                          acacCACGCCAGAACGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTACCAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCa                                                                            <  2:832521/150‑1 (MQ=255)
                                                                           caccacGCCAGAAAGCCCGCCAAACCACGGCCAGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGCCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGCGGCGTTTTTCCACCGCACCGCCGAGATACCCCCACCAGAGCAg                                                                           >  1:208969/1‑150 (MQ=255)
                                                                            accacGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAgg                                                                          <  1:588253/150‑1 (MQ=255)
                                                                            accacGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAgg                                                                          <  2:1622559/150‑1 (MQ=255)
                                                                            accacGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAgg                                                                          <  2:1616083/150‑1 (MQ=255)
                                                                            accacGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAgg                                                                          <  1:1058866/150‑1 (MQ=255)
                                                                                cGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTaa                                                                      >  2:1380624/1‑150 (MQ=255)
                                                                                cGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTaa                                                                      >  1:795352/1‑150 (MQ=255)
                                                                                  ccAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAag                                                                    >  2:397893/1‑150 (MQ=255)
                                                                                     gAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAgg                                                                 >  1:887014/1‑150 (MQ=255)
                                                                                     gAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAgg                                                                 >  1:859160/1‑150 (MQ=255)
                                                                                       aaccccGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGat                                                               <  2:494858/147‑1 (MQ=255)
                                                                                       aaGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTATCAGGGCGCTGATAAGCGTAATGACAATTAGCTAACCGCTTCCGAGGCGTTTTTCCACCTCACCTCCGAGATACAACCACCAGAGCAGGTTTAAGAGGat                                                               >  2:390136/1‑150 (MQ=255)
                                                                                            cGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCa                                                          >  2:1521853/1‑150 (MQ=255)
                                                                                            cGCCAAACCACGGACCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGATAATTCGCTTACCGCTACCGAGGAGTGTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCa                                                          <  1:98633/150‑1 (MQ=255)
                                                                                              ccAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATc                                                        >  1:445529/1‑150 (MQ=255)
                                                                                               cAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTATCAGGGCGCTGATAAGCGTAATGACAATTAGATTACCGCTACCGAGGCGTGTTTCCACCGCACCGCCGAGATTCCACCACCAGAGCAGGTTAAAGAGGTTATGCATCa                                                       >  2:815962/1‑150 (MQ=255)
                                                                                                   ccACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCga                                                   <  1:1184957/150‑1 (MQ=255)
                                                                                                          ccGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGc                                            <  2:859160/150‑1 (MQ=255)
                                                                                                          ccGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGc                                            <  2:887014/150‑1 (MQ=255)
                                                                                                           cGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCCGAGCAGGTGAAAGAGGATATGCATCAGCGAGAAGTGCa                                           >  2:1276915/1‑150 (MQ=255)
                                                                                                           cGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAACTAGCTTACCGCTACCGAGGCGTTTTTCCAACGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCa                                           >  2:215001/1‑150 (MQ=255)
                                                                                                               gAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTACACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGAAATGCATCAGCGAGAAGTGCATTaa                                       <  2:673596/150‑1 (MQ=255)
                                                                                                                 aTTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAAcg                                     <  2:1627166/150‑1 (MQ=255)
                                                                                                                 aTTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAAcg                                     >  2:1474729/1‑150 (MQ=255)
                                                                                                                       gctgcACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTggg                               >  2:1531178/1‑150 (MQ=255)
                                                                                                                       gctgcACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTggg                               >  1:799217/1‑150 (MQ=255)
                                                                                                                       gctgcACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTggg                               >  1:1274138/1‑150 (MQ=255)
                                                                                                                        ctgcACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGt                              >  2:1463463/1‑150 (MQ=255)
                                                                                                                        ctgcACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGt                              <  1:397893/150‑1 (MQ=255)
                                                                                                                        ctgcACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGCATATGCATCAGCGAGAAGTGCATTAACGCGTGGGt                              >  2:1438402/1‑150 (MQ=255)
                                                                                                                         ttcacaTATCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTg                             <  1:431209/148‑1 (MQ=255)
                                                                                                                           cacaTAGCCGCTTAATAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGaa                           <  1:54429/150‑1 (MQ=255)
                                                                                                                            gtcattacGCTTATCAGCGCCCTGTTAAGCGTAATGACAATTAGCTTACCGTTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAg                          >  2:19897/8‑150 (MQ=255)
                                                                                                                            gtcattacGCTTATCAGCGCCCTGTTAAGCGTAATGACAATTAGCTTACCGTTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAg                          >  2:13881/8‑150 (MQ=255)
                                                                                                                            acaTAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAg                          <  1:1019837/150‑1 (MQ=255)
                                                                                                                              cattacGCTTATCAGCGCCCTGTTAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGt                              >  1:705726/6‑144 (MQ=255)
                                                                                                                              cattacGCTTATCAGCGCCCTGTTAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGt                              <  2:705726/139‑1 (MQ=255)
                                                                                                                              cattacGCTTATCAGCGCCCTGTTAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTa                        >  2:411082/6‑150 (MQ=255)
                                                                                                                                aGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAAc                      >  1:1323807/1‑150 (MQ=255)
                                                                                                                                 gCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCATACCGCTACCGAGGCGTTTTTCAACCGCACAGCCGAGATTCCACCACCAGAGAAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCAATAACGCGTGGGTGAAGTacag                     >  2:161224/1‑147 (MQ=255)
                                                                                                                                  ccGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGGTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGc                    <  2:1428716/150‑1 (MQ=255)
                                                                                                                                   cGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGcc                   <  2:21860/150‑1 (MQ=255)
                                                                                                                                    gCTTATCAGCGCCCTGTTAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGCCa                  >  2:1128156/1‑150 (MQ=255)
                                                                                                                                    gCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGCCa                  >  1:1386703/1‑150 (MQ=255)
                                                                                                                                         aCAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGCCAGAACt             <  2:799217/150‑1 (MQ=255)
                                                                                                                                            gggCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGCCAGAACTCaa          <  2:942453/150‑1 (MQ=255)
                                                                                                                                             ggCGCGGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTATCCACCGCCCCGCCGAGATACCACCACCAGAGCAGGGTAAATAGGATATCCCTCCGCGAGAACTTCATAAACGCGTGGGTGAAGTACCGCCAGCACGCaaa         >  1:51104/1‑150 (MQ=255)
                                                                                                                                                   gATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTGCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGCAGAGCATTAACGCGTGGGTGATGTAACGCCAGAACTCGAAATTCAg   >  1:324520/1‑150 (MQ=255)
                                                                                                                                                   gATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGCCAGAACTCAAATTTCAg   <  2:1630891/150‑1 (MQ=255)
                                                                                                                                                   gATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGCCAGAACTCAAATTTCAg   <  1:839938/150‑1 (MQ=255)
                                                                                                                                                    aTAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGCCAGAACTCAAATTTCAgt  <  1:718180/150‑1 (MQ=255)
                                                                                                                                                    aTAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGCCAGAACTCAAATTTCAgt  <  1:1411176/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGCCAGAACTCAAATTTCAGT  >  NC_000913/3560766‑3561063

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: ‑