Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,169,380 T→C 5.5% intergenic (+50/+32) bhsA → / ← ldtC biofilm, cell surface and signaling protein/L,D‑transpeptidase linking Lpp to murein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,169,3800TC 5.5% 368.6 / 8.6 109intergenic (+50/+32)bhsA/ldtCbiofilm, cell surface and signaling protein/L,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand):  ref base T (50/53);  new base C (2/4);  total (52/57)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

Only 100 of 111 total aligned reads displayed.

CTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGC  >  NC_000913/1169232‑1169529
                                                                                                                                                    |                                                                                                                                                     
cTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATg                                                                                                                                                      <  2:591706/150‑1 (MQ=255)
   gcgcAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCg                                                                                                                                                   >  1:966465/1‑150 (MQ=255)
   gcgcAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCg                                                                                                                                                   >  1:1481301/1‑150 (MQ=255)
      cAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCgtgt                                                                                                                                                <  1:1901375/150‑1 (MQ=255)
      cAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCgtgt                                                                                                                                                <  1:413888/150‑1 (MQ=255)
      cAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCgtgt                                                                                                                                                <  1:2501103/150‑1 (MQ=255)
      cAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCgtgt                                                                                                                                                <  1:2326499/150‑1 (MQ=255)
         aaaGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTcgc                                                                                                                                             <  1:1098508/150‑1 (MQ=255)
           aGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCgg                                                                                                                                           >  1:2728468/1‑150 (MQ=255)
           aGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCgg                                                                                                                                           >  1:2728491/1‑150 (MQ=255)
           aGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCgg                                                                                                                                           <  2:2845613/150‑1 (MQ=255)
           aGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCgg                                                                                                                                           <  1:1722647/150‑1 (MQ=255)
           aGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCgg                                                                                                                                           <  2:1211963/150‑1 (MQ=255)
            gCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATCCCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTCCTGCTGATTTTTTCCCCGCGACACGGCATGTCGCggg                                                                                                                                          <  1:2719899/150‑1 (MQ=255)
            gCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCggg                                                                                                                                          <  2:1194659/150‑1 (MQ=255)
                 tGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAttt                                                                                                                                     >  1:2784136/1‑150 (MQ=255)
                     aTGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAttttttt                                                                                                                                 >  1:115287/1‑150 (MQ=255)
                     aTGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATGTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAttttttt                                                                                                                                 >  1:115290/1‑150 (MQ=255)
                       gggCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGCTTTGTTTAt                                                                                                                               >  2:2008329/1‑150 (MQ=255)
                       gggCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTAt                                                                                                                               >  2:919514/1‑150 (MQ=255)
                        ggCGCAAAATCTTTCCGTATTACGTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATc                                                                                                                              <  1:1289605/150‑1 (MQ=255)
                             aaaaTCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGc                                                                                                                         >  1:1035686/1‑150 (MQ=255)
                                 tCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAGGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTa                                                                                                                     <  1:2167059/150‑1 (MQ=255)
                                 tCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTa                                                                                                                     <  2:2665168/150‑1 (MQ=255)
                                 tCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTa                                                                                                                     <  1:1261974/150‑1 (MQ=255)
                                  cTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTAc                                                                                                                    <  1:480518/150‑1 (MQ=255)
                                  cTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTAc                                                                                                                    <  1:480517/150‑1 (MQ=255)
                                    ttCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGAAATGCCGTGTCGCGGGGATTTTTTTATCCGGGATTTACAg                                                                                                                  >  1:2340842/1‑150 (MQ=255)
                                     tCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGc                                                                                                                 >  2:601598/1‑150 (MQ=255)
                                     tCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGc                                                                                                                 <  2:304696/150‑1 (MQ=255)
                                       cGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGt                                                                                                               >  2:2196678/1‑150 (MQ=255)
                                       cGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATAAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATATTTTTATCCGTGCTTTACAGCGt                                                                                                               >  2:191941/1‑150 (MQ=255)
                                        gTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGtt                                                                                                              <  2:2344250/150‑1 (MQ=255)
                                                tCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCt                                                                                                      >  2:2306432/1‑150 (MQ=255)
                                                  tGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAGTAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCa                                                                                                    <  2:1449289/150‑1 (MQ=255)
                                                  tGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCa                                                                                                    <  1:1045825/150‑1 (MQ=255)
                                                  tGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCa                                                                                                    <  1:1045684/150‑1 (MQ=255)
                                                  tGTAACCGGTCAGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTACCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCa                                                                                                    <  2:1354296/150‑1 (MQ=255)
                                                   gTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCAc                                                                                                   <  2:971099/150‑1 (MQ=255)
                                                   gTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCAc                                                                                                   <  1:363608/150‑1 (MQ=255)
                                                    tAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACt                                                                                                  >  2:1843180/1‑150 (MQ=255)
                                                    tAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACt                                                                                                  >  2:2140675/1‑150 (MQ=255)
                                                     aaCCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACtt                                                                                                 <  2:807224/150‑1 (MQ=255)
                                                      aCCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTGTGTGGGCTCACTtg                                                                                                >  2:2160576/1‑150 (MQ=255)
                                                       ccGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtga                                                                                               <  1:834435/150‑1 (MQ=255)
                                                       ccGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtga                                                                                               >  1:2680210/1‑150 (MQ=255)
                                                        cGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgat                                                                                              >  2:1566803/1‑150 (MQ=255)
                                                        cGGTCCGAAGACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgat                                                                                              <  1:1686493/150‑1 (MQ=255)
                                                                 tACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGa                                                                                     <  2:1235785/150‑1 (MQ=255)
                                                                   cccTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACa                                                                                   <  2:1943778/150‑1 (MQ=255)
                                                                    ccTCCATGGAACAGCATGAATTTATAAATAAGCATTAACCCTCATTAATTCCTGCTACTGCTGATGTTTTCCCCGCGACATGCCGTGTCGCGGGGATGTTTATATCCGGGCTGTGCAGCGGGTGGGGGCTCTCTTGATGACGGCGGAcgt                                                                                  >  1:2837220/1‑148 (MQ=255)
                                                                     cTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATc                                                                                 <  2:4647/150‑1 (MQ=255)
                                                                       ccATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGAATTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCa                                                                               <  1:1993088/150‑1 (MQ=255)
                                                                         ctGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCAcc                                                                             <  2:1902121/149‑1 (MQ=255)
                                                                         aTGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCAcc                                                                             <  2:1902049/150‑1 (MQ=255)
                                                                         aTGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCAcc                                                                             <  2:2394254/150‑1 (MQ=255)
                                                                            gAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGCCGGTGGATGTccgccgtc                                                                          <  1:860382/150‑3 (MQ=255)
                                                                            gAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGCCGGTGGATGTccgccgtc                                                                          <  1:860381/150‑3 (MQ=255)
                                                                              aCAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCAt                                                                        >  1:2045251/1‑150 (MQ=255)
                                                                                   agTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACAGCGGACATCCACCGGCATCCCgg                                                                   >  1:1952891/1‑150 (MQ=255)
                                                                                     tAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGac                                                                 >  2:2073173/1‑149 (MQ=255)
                                                                                      aaTTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAAc                                                                <  2:1501531/150‑1 (MQ=255)
                                                                                         ttATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGa                                                             >  1:511254/1‑150 (MQ=255)
                                                                                           ataaataaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAca                                                           <  2:234053/150‑1 (MQ=255)
                                                                                           ataaataaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAca                                                           <  2:376864/150‑1 (MQ=255)
                                                                                           ataaataaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAca                                                           <  2:1895960/150‑1 (MQ=255)
                                                                                               ataaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatcca                                                       <  1:524070/150‑1 (MQ=255)
                                                                                                     aTTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACAt                                                 <  2:2942105/150‑1 (MQ=255)
                                                                                                         aCCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTg                                             >  1:780436/1‑150 (MQ=255)
                                                                                                          cccTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACAACCATCACATGTTGc                                            >  1:742543/1‑150 (MQ=255)
                                                                                                           ccTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCa                                           <  1:2196678/150‑1 (MQ=255)
                                                                                                            cTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATc                                                                       >  1:959183/1‑121 (MQ=255)
                                                                                                            cTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATc                                                                       <  2:959183/121‑1 (MQ=255)
                                                                                                             tCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATc                                         >  1:2124334/1‑150 (MQ=255)
                                                                                                                  aaTGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTc                                    <  1:873538/150‑1 (MQ=255)
                                                                                                                    tGCCTGCTCCTGCGGATTTATTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTGATCCGGCCTTTACAGCGGTTGTGGGCTCACTTGAGGACGGCGGCCATCCACCGGCCTCCCGGAACGCACATCCAGCGCGTGGGGCATCCCTTccg                                  >  2:2169429/1‑148 (MQ=255)
                                                                                                                    tGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAg                                  <  1:1970459/150‑1 (MQ=255)
                                                                                                                    tGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAg                                  >  2:1204053/1‑150 (MQ=255)
                                                                                                                    tGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAg                                  <  1:1765764/150‑1 (MQ=255)
                                                                                                                    tGCCTGCTACTGCTGATTTTGTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAg                                  <  1:411530/150‑1 (MQ=255)
                                                                                                                    tGCCTGCTACTGCTCATTTTTTCCCCGCGACATGCCGTGACGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTTATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCAGCACTTCAg                                  >  2:622518/1‑150 (MQ=255)
                                                                                                                       cTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGGGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGTCGGCGGACATCCACCGGCATCCCGGAACGGTCATCCATCACATGTTGCATCACTTCAGCGt                               >  1:1156568/1‑150 (MQ=255)
                                                                                                                       cTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACCTCCATCACATGTTGCATCACTTCAGCGt                               >  1:1344381/1‑150 (MQ=255)
                                                                                                                       cTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGt                               >  1:1156616/1‑150 (MQ=255)
                                                                                                                            aCTGCTGATTTTGTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGtt                          <  1:479094/150‑1 (MQ=255)
                                                                                                                               gctgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTtgt                       >  1:1513525/1‑150 (MQ=255)
                                                                                                                               gctgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTtgt                       >  1:1637328/1‑150 (MQ=255)
                                                                                                                               gctgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTtgt                       >  1:1513457/1‑150 (MQ=255)
                                                                                                                               gctgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTtgt                       <  1:612507/150‑1 (MQ=255)
                                                                                                                                ctgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTtgtg                      <  1:955737/150‑1 (MQ=255)
                                                                                                                                ctgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACAGCCATCACATGTTGCATCACTTCAGCGTCAGTTtgtg                      >  1:2717517/1‑150 (MQ=255)
                                                                                                                                 tgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCGTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGtgc                     >  2:1680822/1‑150 (MQ=255)
                                                                                                                                  gATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGtgct                    <  2:2680210/150‑1 (MQ=255)
                                                                                                                                           cccGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTa           <  2:2215783/150‑1 (MQ=255)
                                                                                                                                            ccGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTaa          >  2:1594924/1‑150 (MQ=255)
                                                                                                                                            ccGCGACATGCCGGGTCGCGGGGATTTTTTTATCCGGGCGTTACCGCGTGTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCAACCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCACTTTGTGCTGCATCTTTaa          >  2:2359630/1‑150 (MQ=255)
                                                                                                                                            ccGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTaa          >  2:2490526/1‑150 (MQ=255)
                                                                                                                                            ccGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTAAAGCGTTTGTGGGCTCACTTGATGACGGCGGACAACCACCCGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTaa          >  1:1014118/1‑150 (MQ=255)
                                                                                                                                                    tGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGattgc  >  2:2697922/1‑150 (MQ=255)
                                                                                                                                                    tGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGattgc  >  2:2697930/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGC  >  NC_000913/1169232‑1169529

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: