Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,560,910 G→C 7.5% S181R (AGC→AGG glpG ← rhomboid intramembrane serine protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,560,9100GC 7.5% 296.3 / 13.2 94S181R (AGC→AGGglpGrhomboid intramembrane serine protease
Reads supporting (aligned to +/- strand):  ref base G (45/42);  new base C (5/2);  total (50/44)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.89e-01

ATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGCCAGAACTCAAA  >  NC_000913/3560764‑3561056
                                                                                                                                                  |                                                                                                                                                  
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                                                                                                                            <  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: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                                                                                 >  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:1058866/150‑1 (MQ=255)
                                                                              accacGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAgg                                                                   <  1:588253/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: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:1474729/1‑150 (MQ=255)
                                                                                                                   aTTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAAcg                              <  2:1627166/150‑1 (MQ=255)
                                                                                                                         gctgcACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTggg                        >  2:1531178/1‑150 (MQ=255)
                                                                                                                         gctgcACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTggg                        >  1:1274138/1‑150 (MQ=255)
                                                                                                                         gctgcACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTggg                        >  1:799217/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)
                                                                                                                                                  |                                                                                                                                                  
ATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGCCAGAACTCAAA  >  NC_000913/3560764‑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: