Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,382,145 T→C 5.1% intergenic (+35/‑55) degQ → / → degS serine endoprotease, periplasmic/serine endoprotease, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,382,1450TC 5.1% 419.2 / 9.0 139intergenic (+35/‑55)degQ/degSserine endoprotease, periplasmic/serine endoprotease, periplasmic
Reads supporting (aligned to +/- strand):  ref base T (67/63);  new base C (4/3);  total (71/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.19e-01

Only 100 of 143 total aligned reads displayed.

CCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCGGCGCTATTCTGCTGGTTGCCATGCCTTCGCTGCGCAGCCTTAACCCGCT  >  NC_000913/3381996‑3382294
                                                                                                                                                     |                                                                                                                                                     
ccGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGt                                                                                                                                                       >  1:173923/1‑150 (MQ=255)
      tCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGtt                                                                                                                                                 <  1:478255/150‑1 (MQ=255)
      tCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGtt                                                                                                                                                 <  1:213564/150‑1 (MQ=255)
       cGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTa                                                                                                                                                <  2:1416845/150‑1 (MQ=255)
       cGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCACACGTAAGCCTGATGCCCGGttt                                                                                                                                                <  2:27245/150‑2 (MQ=255)
         ggTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCAAATTGTACGCGTCAATGAAAGCATCTATCTGCAGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAAc                                                                                                                                              <  2:1216574/150‑1 (MQ=255)
          gTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACt                                                                                                                                             <  2:1749945/150‑1 (MQ=255)
            taaCTCGATTGCTGAAATGCGTAAAGTGCTGTCGGCAAAACCGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCg                                                                                                                                           <  2:89172/149‑1 (MQ=255)
             aaCTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGt                                                                                                                                          >  2:1427049/1‑150 (MQ=255)
               cTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCAACTTGTACGCGGCAATGAAAGCATCTATCTGCTGCTGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTgg                                                                                                                                        <  1:1949037/150‑1 (MQ=255)
                tCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGt                                                                                                                                       <  2:1815433/150‑1 (MQ=255)
                   aTTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTAtg                                                                                                                                    >  1:822635/1‑150 (MQ=255)
                   aTTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCACCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTAtg                                                                                                                                    <  1:1290925/150‑1 (MQ=255)
                      gCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCGTACGTGGGATGTCCGGTTAACTCGTGGTAtgctg                                                                                                                                 >  1:763334/1‑150 (MQ=255)
                         gAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTAtgctgctg                                                                                                                              >  1:537333/1‑150 (MQ=255)
                         gAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTAtgctgctg                                                                                                                              >  1:719707/1‑150 (MQ=255)
                              gCGTAAAGTGCTGGCTGCAAAACAGTCCATCATAGCCCTGCAAATTGTACGCGGCAATTAAAGAATCTATCTGCTGATGCGTTAATGTCGTAAACAGGGAATCATGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGtt                                                                                                                         <  1:1305562/150‑1 (MQ=255)
                              gCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGtt                                                                                                                         >  1:1534115/1‑150 (MQ=255)
                              gCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTAGCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCggt                                                                                                                         >  1:1333663/1‑148 (MQ=255)
                               cGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAATCCGGGCATCAGGCTTACGTGTGATGTGCGGTTAACTCGTGGTATGCTGCTGCCGTTc                                                                                                                        >  2:1360341/1‑150 (MQ=255)
                                 cacaGTGCGGGCGGCAAAACCGGCCATCATCGCCCTGCAAATTGTACGCGGCACTGAAAGCATCTATCTGCTGCTGCGTTCATGTCGTAAACCGGGCAGCAGGCTTACGGGTGCAGTCCGGTTAACTCGTGGTATGCTGCAGCCGTTccc                                                                                                                      <  1:244017/147‑1 (MQ=255)
                                       gCTGGCGGCAAAACCGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCtttttt                                                                                                                <  2:782235/150‑1 (MQ=255)
                                        cTGGCGGCAAAACCGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGAATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACACGTGGTATGCTGCTGCCGTTCCCTATTTTa                                                                                                               >  2:570248/1‑150 (MQ=255)
                                         tGGCGGCAAAACCGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTaa                                                                                                              >  1:526749/1‑150 (MQ=255)
                                          ggcggcAAAACCGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAAt                                                                                                             >  2:1076935/1‑150 (MQ=255)
                                           gcggcAAAACCGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCAGCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATg                                                                                                            <  1:681045/150‑1 (MQ=255)
                                            ctgcAAAACCGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGa                                                                                                           <  1:1862316/148‑1 (MQ=255)
                                             ggcAAAACCGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGAc                                                                                                          >  2:712739/1‑150 (MQ=255)
                                             ggcAAAACCGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGAc                                                                                                          >  1:1263005/1‑150 (MQ=255)
                                             ggcAAAACCGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGATGCCGTTCCCTTTTTTAATGAc                                                                                                          >  1:1182023/1‑150 (MQ=255)
                                              gcAAAACCGGCCATCCTCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGGGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACg                                                                                                         <  2:132757/150‑1 (MQ=255)
                                              gcAAAACCGGCCATCAACGCCCTGCAAATTCTACGCGGCAATGAAAGCATCTAACTGCTGATGCGGTAATTTCGTAAACCGGTCATCAGGCTTATGTGTGATGTCCGGTTAACTCGAGCTATGCTGCTGCCGTTCCCTTATTTAATGagg                                                                                                         >  1:1330630/1‑148 (MQ=255)
                                               cAAAACCGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACgg                                                                                                        >  1:1216305/1‑149 (MQ=255)
                                               cAAAACCGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGc                                                                                                        >  1:736737/1‑150 (MQ=255)
                                                aaaaCCGGCCATCATCGCACTGCAAATTGTACGCGGCAAGGAAAGCATCGATCTGCTGATGCGTTAATGGCGTAAACCGGTCAGCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGcc                                                                                                       <  1:1820471/150‑1 (MQ=255)
                                                  aaCCGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTc                                                                                                     >  2:740978/1‑150 (MQ=255)
                                                    ccGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAACCCGGGCATCCGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCca                                                                                                   <  1:436072/150‑1 (MQ=255)
                                                    ccGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCca                                                                                                   <  1:428730/150‑1 (MQ=255)
                                                     cGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCcat                                                                                                  >  1:795050/1‑150 (MQ=255)
                                                     cGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTACCCTTTTTTAATGACGCCTCcat                                                                                                  >  1:962036/1‑150 (MQ=255)
                                                      ggCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCATGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCcatc                                                                                                 <  1:835748/150‑1 (MQ=255)
                                                       gCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCcatca                                                                                                <  1:557041/150‑1 (MQ=255)
                                                       gCCATCAGCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCGGGCTTACGTGTGATGTCCGGTTAACTCGGGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCcatca                                                                                                <  1:984857/150‑1 (MQ=255)
                                                        ccATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAAATCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCcatcat                                                                                               >  2:731032/1‑150 (MQ=255)
                                                          atcatcGCCCTGCAAATTGTACGCGGCAAAGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGt                                                                                             <  2:339588/150‑1 (MQ=255)
                                                           tcatcGCCCTGCAAATTGTACGCGGCAATGAAAGCTTCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGtt                                                                                            <  2:500996/150‑1 (MQ=255)
                                                            catcGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGttt                                                                                           >  1:1639841/1‑150 (MQ=255)
                                                            catcGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGttt                                                                                           >  1:1954364/1‑150 (MQ=255)
                                                             atcGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGTGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTtg                                                                                          <  1:414836/150‑1 (MQ=255)
                                                             atcGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTtg                                                                                          >  1:25941/1‑150 (MQ=255)
                                                             atcGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTtg                                                                                          <  2:537333/150‑1 (MQ=255)
                                                               cGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTATTTAATGACGCCTCCATCATGTTtgtg                                                                                        >  2:704970/1‑150 (MQ=255)
                                                                gCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGa                                                                                       >  2:1495588/1‑150 (MQ=255)
                                                                 cccTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGaa                                                                                      <  2:567126/150‑1 (MQ=255)
                                                                  ccTGCAAATTGTACGCGGAAATGAAAGCATATATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGTGTTACGTGTGATGTCCTGTTAACTCGTGGTATTCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAg                                                                                     <  2:956523/150‑1 (MQ=255)
                                                                    tGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGct                                                                                   >  2:904515/1‑150 (MQ=255)
                                                                     gCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGctc                                                                                  >  2:1361945/1‑150 (MQ=255)
                                                                       aaaTTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCtt                                                                                >  1:399335/1‑150 (MQ=255)
                                                                         aTTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGGCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTAc                                                                              <  1:1457156/150‑1 (MQ=255)
                                                                                     aaTGAAAGAATATATCTGCTGATGCTTTAATGTCGTTAACCGGGCATCAGGCTTACGTGTGATGTCCTGTTAACTCGTTGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGa                                                                  <  2:426998/150‑1 (MQ=255)
                                                                                        gAAAGCATCTATCTGCTGATGCTTTAATGTCTGAAACCGGGCATCAGGCTTACGTGTGATGTGCGGTTAACTCGTGGTATGCTGCTGACGTTCCATTTTTTAATGTCGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTg                                                               <  2:1003811/150‑1 (MQ=255)
                                                                                        gAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGGGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTg                                                               <  1:1107776/150‑1 (MQ=255)
                                                                                          aaGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGa                                                             <  2:1411945/150‑1 (MQ=255)
                                                                                           agcagcTATCTGCTGATGCGTTAAGTGCGTCACCCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTAGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGTTTGGAt                                                            <  1:106686/150‑1 (MQ=255)
                                                                                            gCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGAtt                                                           <  2:1986519/150‑1 (MQ=255)
                                                                                              atctttctGTTGTTGCGTTAATGTCGTAAACCGGGGATCATGCTTACGGGTGATGGCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGAATaa                                                         <  1:209704/150‑1 (MQ=255)
                                                                                                    ctgctgATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTTCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCg                                                   <  2:194048/150‑1 (MQ=255)
                                                                                                    ctgctgATGCGTTAATGTCATAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCg                                                   >  2:443213/1‑150 (MQ=255)
                                                                                                     tgctgATGCGTTAATGTCGTAAACCGGGCAGCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCgg                                                  <  1:1410462/150‑1 (MQ=255)
                                                                                                          aTGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCGGCGCTa                                             <  2:1660082/150‑1 (MQ=255)
                                                                                                             cGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCGGCGCTATTc                                          >  1:374620/1‑150 (MQ=255)
                                                                                                             cGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCGGCGCTATTc                                          >  1:740616/1‑150 (MQ=255)
                                                                                                                   tGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCGGCGCTATTCTGCTgg                                    <  2:1768058/150‑1 (MQ=255)
                                                                                                                      cGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCGGCGCTATTCTGCTGGTTg                                 >  2:940819/1‑150 (MQ=255)
                                                                                                                      cGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCGGCGCTATTCTGCTGGTTg                                 >  1:1782717/1‑150 (MQ=255)
                                                                                                                      cGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCGGCGCTATTCTGCTGGTTg                                 >  1:884472/1‑150 (MQ=255)
                                                                                                                      agttaaCCGGACATCACACGTAAGCCTGATGCCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCGGCGCTATTCTGCTGGTTg                                 >  2:1410852/5‑150 (MQ=38)
                                                                                                                      agttaaCCGGACATCACACGTAAGCCTGATGCCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTGTTTTAATGAAGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGAGATTGGATTAATTGTCGTCGCTATTCTGCTGGTTg                                 >  1:1005549/5‑150 (MQ=37)
                                                                                                                       gttaaCCGGACATCACACGTAAGCCTGATGCCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCGGCGCTATTCTGCTGGTTGc                                >  2:696661/4‑150 (MQ=38)
                                                                                                                       gTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCGGCGCTATTCTGCTGGTTGc                                >  1:1119009/1‑150 (MQ=255)
                                                                                                                        tAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCGGCGCTATTCTGCTGGTTGcc                               >  2:348308/1‑150 (MQ=255)
                                                                                                                         aaaCCGGGCATCAGGCTTACGTGTGAGGTCCGGTTAACTCGTGGTATGCTGCAGCCGTGCCCTTTTTTAATGACGCCTCAATCATGTTTGTGACGCTCTTACGTTCCGTTGCGATTTGATTAATTGTCGGCGCTATTCTGCTGGTTGCCa                              >  1:1531671/1‑150 (MQ=255)
                                                                                                                          aaCCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCGGCGCTATTCTGCTGGTTGCCAt                             <  1:175085/150‑1 (MQ=255)
                                                                                                                            ccGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTGTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCGGCGCTATTCTGCTGGTTGCCATGc                           >  2:1390751/1‑150 (MQ=255)
                                                                                                                             cGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCGGCGCTATTCTGCTGGTTGCCATGcc                          <  1:231302/150‑1 (MQ=255)
                                                                                                                             cGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCGGCGCTATTCTGCTGGTTGCCATGcc                          <  2:1609191/150‑1 (MQ=255)
                                                                                                                             cGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTAAGTTCCGTTGCGATTGGATTAATTGTCGGCGCTATTCTGCTGGTTGCCATGcc                          >  2:79107/1‑150 (MQ=255)
                                                                                                                                  aTCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCGGCGCTATTCTGCTGGTTGCCATGCCTTCGc                     >  2:474107/1‑150 (MQ=255)
                                                                                                                                   tCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCGGCGCTATTCTGCTGGTTGCCATGCCTTCGCt                    >  2:635091/1‑150 (MQ=255)
                                                                                                                                    cAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCGGCGCTATTCTGCTGGTTGCCATGCCTTCGCTg                   <  1:1866266/150‑1 (MQ=255)
                                                                                                                                      ggCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCGGCGCTATTCTGCTGGTTGCCATGCCTTCGCTgca                 >  1:731580/1‑149 (MQ=255)
                                                                                                                                       gCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCGGCGCTATTCTGCTGGTTGCCATGCCTTCGCTgcgc                <  1:1076935/150‑1 (MQ=255)
                                                                                                                                           aCGTGGGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATGGGATTAATTGTCGGCGCTATTCTGCTGGTTGCCATGCCTTCGCTGCGCAGcc            <  2:1442826/150‑1 (MQ=255)
                                                                                                                                              tgtgATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCTTTGAGATTGGATTAATTGTCGTCGCTATTCTGCTGGTTGACATGCCTTCGCTGCGCAGCCtaa         >  1:340152/1‑148 (MQ=255)
                                                                                                                                              tgtgATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGACTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCGGCGCTATTCTGCTGGTTGCCATGCCTTCGCTGCGCAGCCTTa         >  1:592995/1‑150 (MQ=255)
                                                                                                                                                tgATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCGGCGCTATTCTGCTGGTTGCCATGCCTTCGCTGCGCAGCCTTAAc       <  2:1867677/150‑1 (MQ=255)
                                                                                                                                                 gATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCGGCGCTATTCTGCTGGTTGCAATGCCTTCGCTGCGCAGCCTTAAcc      >  2:800028/1‑150 (MQ=255)
                                                                                                                                                    gTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCGGCGCTATTCTGCTGGTTGCCATTCCTTCGCTGCGCAGCCTTAACCCGc   >  2:1954933/1‑150 (MQ=255)
                                                                                                                                                     tCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCGGCGCTATTCTGCTGGTTGCCATGCCTTCGCTGCGCAGCCTTAACCCGCt  <  1:1860568/150‑1 (MQ=255)
                                                                                                                                                     tCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCGGCGCTATTCTGCTGGTTGCCATGCCTTCGCTGCGCAGCCTTAACCCGCt  <  2:1159333/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCAAATTGTACGCGGCAATGAAAGCATCTATCTGCTGATGCGTTAATGTCGTAAACCGGGCATCAGGCTTACGTGTGATGTCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCGGCGCTATTCTGCTGGTTGCCATGCCTTCGCTGCGCAGCCTTAACCCGCT  >  NC_000913/3381996‑3382294

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