Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,560,907 G→C 7.5% A182A (GCC→GCG glpG ← rhomboid intramembrane serine protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,560,9070GC 7.5% 286.6 / 12.7 94A182A (GCC→GCGglpGrhomboid intramembrane serine protease
Reads supporting (aligned to +/- strand):  ref base G (44/43);  new base C (5/2);  total (49/45)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.56e-01

AGATAATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGCCAGAACTCAAA  >  NC_000913/3560759‑3561056
                                                                                                                                                    |                                                                                                                                                     
aGATAATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAggg                                                                                                                                                      <  2:957134/150‑1 (MQ=255)
aGATAATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAggg                                                                                                                                                      <  2:472956/150‑1 (MQ=255)
 gATAATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGgc                                                                                                                                                     <  1:1719447/150‑1 (MQ=255)
 gATAATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGgc                                                                                                                                                     <  2:403337/150‑1 (MQ=255)
     aTTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTg                                                                                                                                                 <  1:630572/150‑1 (MQ=255)
       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                                                                                                                            >  1:1515170/1‑150 (MQ=255)
                          aTGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGc                                                                                                                            <  2:1037242/150‑1 (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                                                                                 >  1:21860/1‑150 (MQ=255)
                                                                     ccATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATacc                                                                                 >  2:839938/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                                                                   <  1:1058866/150‑1 (MQ=255)
                                                                                   accacGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAgg                                                                   <  2:1622559/150‑1 (MQ=255)
                                                                                   accacGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAgg                                                                   <  2:1616083/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:859160/1‑150 (MQ=255)
                                                                                            gAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAgg                                                          >  1:887014/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:887014/150‑1 (MQ=255)
                                                                                                                 ccGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGc                                     <  2:859160/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                        >  1:1274138/1‑150 (MQ=255)
                                                                                                                              gctgcACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTggg                        >  1:799217/1‑150 (MQ=255)
                                                                                                                              gctgcACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTggg                        >  2:1531178/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                       <  2:705726/139‑1 (MQ=255)
                                                                                                                                     cattacGCTTATCAGCGCCCTGTTAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGt                       >  1:705726/6‑144 (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)
                                                                                                                                                    |                                                                                                                                                     
AGATAATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGCCAGAACTCAAA  >  NC_000913/3560759‑3561056

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: