Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,169,380 T→C 6.8% 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 6.8% 338.9 / 11.6 118intergenic (+50/+32)bhsA/ldtCbiofilm, cell surface and signaling protein/L,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand):  ref base T (63/46);  new base C (4/4);  total (68/50)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

Only 100 of 120 total aligned reads displayed.

GCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTG  >  NC_000913/1169231‑1169528
                                                                                                                                                     |                                                                                                                                                    
gcTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACAt                                                                                                                                                      >  2:568495/1‑150 (MQ=255)
   ggCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGcc                                                                                                                                                   >  1:581368/1‑150 (MQ=255)
   ggCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGcc                                                                                                                                                   >  2:875905/1‑150 (MQ=255)
    gcgcAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCg                                                                                                                                                  >  1:636430/1‑150 (MQ=255)
    gcgcAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACAGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCg                                                                                                                                                  >  1:1737989/1‑150 (MQ=255)
     cgcAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTAATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCgt                                                                                                                                                 >  1:1175360/1‑150 (MQ=255)
       cAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCATCCATGGAACAGAAGTAATTTATAAATAAGCATTAACCATCATTAATGCCTGATACTGCTGATTTTTTCCCCGCGACATGCCgtgt                                                                                                                                               >  2:1611177/1‑150 (MQ=255)
         aaaaGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTTTAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTcg                                                                                                                                             <  1:781055/150‑1 (MQ=255)
            aGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCgg                                                                                                                                          <  1:160162/150‑1 (MQ=255)
             gCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCggg                                                                                                                                         <  2:1280794/150‑1 (MQ=255)
             gCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCggg                                                                                                                                         <  1:51776/150‑1 (MQ=255)
             gCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCggg                                                                                                                                         <  1:1369767/150‑1 (MQ=255)
             gCGGATGAGATGGGCGCAAAATCTTTCCGGATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTATCCCCGCGACATGCCGTGTCGCggg                                                                                                                                         <  1:320578/150‑1 (MQ=255)
               ggATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGTGGa                                                                                                                                       >  1:645369/1‑150 (MQ=255)
                   gcgaTGGGCGCAAAAGCTGTTCGTATTACTTTTGTATCCGGTCCGTATACCCTCCATCGATCAGCAGTAATTTATAAATAAGCCTTATCCCTGATTATTGCCTGGGACTGGTGATTTTTTCCCCGCGCCATGCCGTGTCGCGGGGAtttt                                                                                                                                   <  2:1199047/148‑1 (MQ=255)
                     gaTGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCGTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAtttttt                                                                                                                                 >  2:327058/1‑150 (MQ=255)
                     gaTGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAtttttt                                                                                                                                 <  1:710070/150‑1 (MQ=255)
                       tGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGGCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTa                                                                                                                               <  2:284592/150‑1 (MQ=255)
                           cgaaaatTCTTTCCTTATTACTTCTGTAACCGGTCCTAATACCCTCCATGGAACTTCAGTATTTTAATATTATGCATTAACCCTCATTTATGCCTGCTACGGCTTATTTTTACGCCGCGACATGCCGTGTCGCGGTGATTTTTTTTTCCg                                                                                                                           <  2:1048130/147‑1 (MQ=255)
                            gcAAAATCTTTCCGTTTTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCgg                                                                                                                          <  1:274748/150‑1 (MQ=255)
                             cAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCAAGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCggg                                                                                                                         >  1:155768/1‑150 (MQ=255)
                              aaaaTCTTTCCGTATTACTTCTGTAACCGGTGCGAATACCCTCCATGGATCAGCATTAATTTATAAATTAGCCTTAACGCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGc                                                                                                                        <  2:1282845/150‑1 (MQ=255)
                              aaaaTCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGc                                                                                                                        >  2:1873267/1‑150 (MQ=255)
                                 aTCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCttt                                                                                                                     >  2:1830271/1‑150 (MQ=255)
                                    tttCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATGTTTTTATCCGGGCTTTACa                                                                                                                  >  2:711972/1‑150 (MQ=255)
                                      tCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGc                                                                                                                >  2:524571/1‑150 (MQ=255)
                                       ccGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCg                                                                                                               >  1:710102/1‑150 (MQ=255)
                                        cGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTTTCCGGGCTTTACAGCGt                                                                                                              >  1:1479215/1‑150 (MQ=255)
                                        cGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGt                                                                                                              >  1:213209/1‑150 (MQ=255)
                                        agTATTACTTCTGTAACCGTTACGAATAACCTCGATGGAACAGATGTAATTTATAAATAAGCATTAAACCTAATTAATGCCTTCTACTGCTGATTTTTTCCCCGAGACATGCCTTGTATCGTGTATTTTTTTATCAGGGCTTTACAGCGt                                                                                                              <  2:484362/149‑1 (MQ=255)
                                           aTTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAATCATTAACCCTCATTAATGCCTGCTACTGCTGTGTTTTTCCCCTCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGGTTTACAGCGTTtg                                                                                                           <  2:297353/150‑1 (MQ=255)
                                            ttACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTtgt                                                                                                          >  1:1423738/1‑150 (MQ=255)
                                              aCTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTgg                                                                                                        >  1:340884/1‑150 (MQ=255)
                                              aCTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTgg                                                                                                        >  1:1710915/1‑150 (MQ=255)
                                               cTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTggg                                                                                                       >  2:847636/1‑150 (MQ=255)
                                               cTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTggg                                                                                                       >  1:877250/1‑150 (MQ=255)
                                                  cTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTc                                                                                                    >  2:1550921/1‑150 (MQ=255)
                                                   tGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCa                                                                                                   >  1:1430220/1‑150 (MQ=255)
                                                     tAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCCATAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACt                                                                                                 <  2:286103/150‑1 (MQ=255)
                                                          ggTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgatg                                                                                            <  2:373365/150‑1 (MQ=255)
                                                          ggTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgatg                                                                                            <  2:1403294/150‑1 (MQ=255)
                                                                aaTACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTCCTGCTGATTTTTTCCCCGCGACATCCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGAcggcg                                                                                      <  2:336449/150‑1 (MQ=255)
                                                                 aTACCCTCCATGGAACAGCATTAATTTATAAATAAGCATTAACCCTCATTAATTCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGAcggcgg                                                                                     <  1:1210055/150‑1 (MQ=255)
                                                                     ccTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTGTTTCCCCGCGACATGCCGGGTCGCGGGGATGGTTTTATCCGGGCTTTACAGCGTGTGTGGGCTCACTTGATGCCGGCGGACAt                                                                                 >  1:1303754/1‑150 (MQ=255)
                                                                      cTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATc                                                                                >  1:437616/1‑150 (MQ=255)
                                                                      cTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATc                                                                                >  1:843627/1‑150 (MQ=255)
                                                                       tCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATcc                                                                               >  1:464039/1‑150 (MQ=255)
                                                                       tCCATGGAACAGCAGTAATGTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGCCATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGGTCCCTTGATGACGGCGGACATcc                                                                               >  2:343197/1‑150 (MQ=255)
                                                                         cATGGAACAGCATAAATTTATAAATAAGCATTAACCCTCATTAATGCCTGGTACTGCTGATTTTTTCCCCGCGACATGCCGAGGCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCAc                                                                             <  1:1744864/150‑1 (MQ=255)
                                                                            ggAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGCCATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCgg                                                                          >  2:1554863/1‑150 (MQ=255)
                                                                                 cgactgaaTTTATAACTAAGCATAACCCCTCATAAATGCCTGCGACTGCTTATTATTTCCCCGCGACATGCCGTGGCGCGGTGATTTTTATATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGCCGGTGGATGTCCGCCGTCATca                                                                     <  2:301355/144‑2 (MQ=255)
                                                                                 agcagTAATTTATAAATAAGCATTAACCCTCATTACTGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGCCGGTGGATGTCCGCCGTCATca                                                                     <  1:833887/150‑2 (MQ=255)
                                                                                  gcagTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGAGGGCTCACTTGATGACGGCGGACATCCACCGGCATccc                                                                    >  1:1754098/1‑150 (MQ=255)
                                                                                        aTTTATAAATACGCCCTACCCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGGCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGTCATCCCGGAACg                                                              <  1:15963/150‑1 (MQ=255)
                                                                                        aTTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACg                                                              <  1:493288/150‑1 (MQ=255)
                                                                                         tttATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACgg                                                             <  2:4430/150‑1 (MQ=255)
                                                                                              aaataaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatc                                                        <  2:636430/150‑1 (MQ=255)
                                                                                              aaataaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatc                                                        <  2:1737989/150‑1 (MQ=255)
                                                                                               aatacGCAATAACCCTCATTAATGCCTGCTCCTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatcc                                                       <  2:349564/150‑1 (MQ=255)
                                                                                               aataaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatcc                                                       >  2:860526/1‑150 (MQ=255)
                                                                                               aataaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatcc                                                       >  1:74489/1‑150 (MQ=255)
                                                                                                 taaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatccat                                                     >  2:395135/1‑150 (MQ=255)
                                                                                                  aaGCATTAACCATCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatccatc                                                    <  2:558371/150‑1 (MQ=255)
                                                                                                   aGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatccatcc                                                   >  2:1655534/1‑149 (MQ=255)
                                                                                                    gCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATcac                                                  >  2:14828/1‑150 (MQ=255)
                                                                                                     cATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGGCATCCACCGGCATCCCGGAACGGACATCCATcaca                                                 >  2:357017/1‑150 (MQ=255)
                                                                                                      aTTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCCCCGGCATCCCGGAACGGACCTCCATCACAt                                                >  2:1687585/1‑150 (MQ=255)
                                                                                                      aTTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGGACGGACATCCATCACAt                                                >  2:1738450/1‑150 (MQ=255)
                                                                                                        tAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTTTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGt                                              <  2:1782794/150‑1 (MQ=255)
                                                                                                         aaCCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGtt                                             <  2:1143459/150‑1 (MQ=255)
                                                                                                         aaCCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGtt                                             <  2:464039/150‑1 (MQ=255)
                                                                                                           cccTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGc                                           >  1:524750/1‑150 (MQ=255)
                                                                                                            ccTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCa                                          >  1:751235/1‑150 (MQ=255)
                                                                                                            ccTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGACGTGTCGCGGGGATTTTTTTATCAGGGCTTTACAGCTTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCa                                          >  1:354327/1‑150 (MQ=255)
                                                                                                                    aTGGCTGCTACTGCTGATTTTTTGCCCGCGAGATGCCGTGGCGCGGGGATTTTTTTATGCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCa                                  <  1:1258320/150‑1 (MQ=255)
                                                                                                                      gcatGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTAGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGc                                <  1:327058/147‑1 (MQ=255)
                                                                                                                      gCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACAGCCATCACATGTTGCATCCCTTCAGc                                >  1:62975/1‑150 (MQ=255)
                                                                                                                       ccTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCg                               <  1:1873267/150‑1 (MQ=255)
                                                                                                                       ccTGCTACTGCTGATTTTTTCCCCGCGACATGCCGGGGCGGGGGGAATTTTTTATCCGGGCGTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACCTACACCCGCCTACCGGAACGGACCTACATGACCTGTTTGATCACGTGAGCg                               >  2:1635640/1‑150 (MQ=255)
                                                                                                                                gctgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGATATCCATCACATGTTGCATCACTTCAGCGTCAGTTtgt                      >  1:677973/1‑150 (MQ=255)
                                                                                                                                gctgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTtgt                      >  2:1211222/1‑150 (MQ=255)
                                                                                                                                 ctgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTtgtg                     <  2:581368/150‑1 (MQ=255)
                                                                                                                                 ctgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGAGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTtgtg                     >  2:1463442/1‑150 (MQ=255)
                                                                                                                                  tgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTCTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGCCGGCGGACATCCACCGGCATCCCGGAACGGACCTCCATCACATGTTGCATCACTGCAGCGTCAGTTTGtgc                    >  1:959129/1‑150 (MQ=255)
                                                                                                                                   gATTTTTTCCCCGCGACATGGCGTGTAGCGGGGATTTTTTTATCCGGGATTTACAGCGTTAGTGGGCTCACTTGATGAAGGCGGACATCCACCGGCATCCCGGAACGGACATCCATAACATGTTGCATCACTTCAGCGTCAGTTTGtgct                   >  1:1803638/1‑150 (MQ=255)
                                                                                                                                   gATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGtgct                   >  1:1084425/1‑150 (MQ=255)
                                                                                                                                   gATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGtgct                   <  2:877250/150‑1 (MQ=255)
                                                                                                                                     aaaaatCCCCGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGtgctgc                 >  2:174928/6‑150 (MQ=255)
                                                                                                                                      tttttCCCCGCGACATGCCGTGGCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCa                <  1:183350/150‑1 (MQ=255)
                                                                                                                                        aatCCCCGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATc              >  1:1334394/3‑150 (MQ=255)
                                                                                                                                           ccccGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCttt           >  1:642024/1‑150 (MQ=255)
                                                                                                                                            cccGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTa          >  2:1292031/1‑150 (MQ=255)
                                                                                                                                            cccGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTa          <  1:870534/150‑1 (MQ=255)
                                                                                                                                            cccGCGACACGGCATGTCGCGGGTATTATTATATCCGGGCTTTACAGCGTTTGTTTGCTCACTTCATGACGGCGGACAACCACCCCCATCCCGGACCGGCCATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTa          >  2:1024176/1‑150 (MQ=255)
                                                                                                                                            cccGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCACGGAACGGACCTCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTa          >  1:320536/1‑150 (MQ=255)
                                                                                                                                             ccGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTaa         >  1:1655952/1‑150 (MQ=255)
                                                                                                                                              cgcgACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTaaa        >  2:1007759/1‑150 (MQ=255)
                                                                                                                                               gccaCATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAAt       <  2:155768/147‑1 (MQ=255)
                                                                                                                                                  aCATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGat    <  1:1550921/150‑1 (MQ=255)
                                                                                                                                                    aTGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGattg  <  1:1714390/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTG  >  NC_000913/1169231‑1169528

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: