Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,560,903 A→T 9.0% L184I (TTA→ATA)  glpG ← rhomboid intramembrane serine protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,560,9030AT 9.0% 244.0 / 9.9 97L184I (TTA→ATA) glpGrhomboid intramembrane serine protease
Reads supporting (aligned to +/- strand):  ref base A (42/44);  new base T (7/2);  total (50/47)
Fisher's exact test for biased strand distribution p-value = 1.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

GCAAAGATAATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGCCAGAACT  >  NC_000913/3560755‑3561052
                                                                                                                                                    |                                                                                                                                                     
gcAAAGATAATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAAc                                                                                                                                                      >  1:1037242/1‑150 (MQ=255)
  aaaGATAATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAg                                                                                                                                                    <  2:468412/150‑1 (MQ=255)
    aGATAATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAggg                                                                                                                                                  <  2:957134/150‑1 (MQ=255)
    aGATAATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAggg                                                                                                                                                  <  2:472956/150‑1 (MQ=255)
     gATAATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGgc                                                                                                                                                 <  2:403337/150‑1 (MQ=255)
     gATAATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGgc                                                                                                                                                 <  1:1719447/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:1616083/1‑150 (MQ=255)
                         ggTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACaa                                                                                                                             >  1:1622559/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:825900/150‑1 (MQ=255)
                                                      ccACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTcc                                                                                                <  1:918357/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                                                               <  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                    >  2:1531178/1‑150 (MQ=255)
                                                                                                                                  gctgcACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTggg                    >  1:1274138/1‑150 (MQ=255)
                                                                                                                                  gctgcACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTggg                    >  1:799217/1‑150 (MQ=255)
                                                                                                                                   ctgcACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGt                   <  1:397893/150‑1 (MQ=255)
                                                                                                                                   ctgcACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGt                   >  2:1463463/1‑150 (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:13881/8‑150 (MQ=255)
                                                                                                                                       gtcattacGCTTATCAGCGCCCTGTTAAGCGTAATGACAATTAGCTTACCGTTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAg               >  2:19897/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)
                                                                                                                                                    |                                                                                                                                                     
GCAAAGATAATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGCCAGAACT  >  NC_000913/3560755‑3561052

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: