Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,169,384 G→A 6.9% intergenic (+54/+28) 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,3840GA 6.9% 349.5 / 8.9 115intergenic (+54/+28)bhsA/ldtCbiofilm, cell surface and signaling protein/L,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand):  ref base G (59/47);  new base A (4/4);  total (63/52)
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.00e+00

Only 100 of 119 total aligned reads displayed.

GCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCAT  >  NC_000913/1169235‑1169531
                                                                                                                                                     |                                                                                                                                                   
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                                                                                                                                            <  1:51776/150‑1 (MQ=255)
         gCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCggg                                                                                                                                            <  1:1369767/150‑1 (MQ=255)
         gCGGATGAGATGGGCGCAAAATCTTTCCGGATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTATCCCCGCGACATGCCGTGTCGCggg                                                                                                                                            <  1:320578/150‑1 (MQ=255)
         gCGGATGAGAAGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTGCGATTACCCTCCATGGAACAGCAGTAATTTATACATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCggg                                                                                                                                            <  2:362763/150‑1 (MQ=255)
           ggATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGTGGa                                                                                                                                          >  1:645369/1‑150 (MQ=255)
               gcgaTGGGCGCAAAAGCTGTTCGTATTACTTTTGTATCCGGTCCGTATACCCTCCATCGATCAGCAGTAATTTATAAATAAGCCTTATCCCTGATTATTGCCTGGGACTGGTGATTTTTTCCCCGCGCCATGCCGTGTCGCGGGGAtttt                                                                                                                                      <  2:1199047/148‑1 (MQ=255)
                agaTGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAttttt                                                                                                                                     <  1:774080/150‑1 (MQ=255)
                 gaTGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCGTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAtttttt                                                                                                                                    >  2:327058/1‑150 (MQ=255)
                 gaTGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAtttttt                                                                                                                                    <  1:710070/150‑1 (MQ=255)
                   tGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGGCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTa                                                                                                                                  <  2:284592/150‑1 (MQ=255)
                       cgcAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCg                                                                                                                              <  2:48707/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:1869612/150‑1 (MQ=255)
                                  tCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGc                                                                                                                   >  2:524571/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)
                                             tCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCt                                                                                                        >  2:686777/1‑150 (MQ=255)
                                               tGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCa                                                                                                      >  1:1430220/1‑150 (MQ=255)
                                                   aCCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGCCAGGCAGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGGTTGTGGGCTCACTtg                                                                                                  >  1:904873/1‑150 (MQ=255)
                                                   aCCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtg                                                                                                  <  2:131951/150‑1 (MQ=255)
                                                     cGGTCCGAATAACCTCCATTGAACAGCAGTAATTTATAAATAATCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgat                                                                                                <  2:305230/150‑1 (MQ=255)
                                                      ggTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgatg                                                                                               <  2:1403294/150‑1 (MQ=255)
                                                      ggTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgatg                                                                                               <  2:373365/150‑1 (MQ=255)
                                                       gTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgatga                                                                                              >  2:1720287/1‑150 (MQ=255)
                                                          cGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGAcgg                                                                                           >  2:1556101/1‑150 (MQ=255)
                                                            aaTACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTCCTGCTGATTTTTTCCCCGCGACATCCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGAcggcg                                                                                         <  2:336449/150‑1 (MQ=255)
                                                             aTACCCTCCATGGAACAGCATTAATTTATAAATAAGCATTAACCCTCATTAATTCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGAcggcgg                                                                                        <  1:1210055/150‑1 (MQ=255)
                                                                 ccTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTGTTTCCCCGCGACATGCCGGGTCGCGGGGATGGTTTTATCCGGGCTTTACAGCGTGTGTGGGCTCACTTGATGCCGGCGGACAt                                                                                    >  1:1303754/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:1737989/150‑1 (MQ=255)
                                                                                           aataaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatcc                                                          >  1:74489/1‑150 (MQ=255)
                                                                                           aataaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatcc                                                          >  2:860526/1‑150 (MQ=255)
                                                                                             taaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatccat                                                        >  2:395135/1‑150 (MQ=255)
                                                                                              aaGCATTAACCATCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatccatc                                                       <  2:558371/150‑1 (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)
                                                                                                     aaCCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGtt                                                <  2:1143459/150‑1 (MQ=255)
                                                                                                       cccTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGc                                              >  1:524750/1‑150 (MQ=255)
                                                                                                        ccTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCa                                             >  1:1793767/1‑150 (MQ=255)
                                                                                                        ccTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCa                                             >  1:751235/1‑150 (MQ=255)
                                                                                                        ccTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGACGTGTCGCGGGGATTTTTTTATCAGGGCTTTACAGCTTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCa                                             >  1:354327/1‑150 (MQ=255)
                                                                                                            aTTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCAc                                         >  1:661185/1‑150 (MQ=255)
                                                                                                                aTGGCTGCTACTGCTGATTTTTTGCCCGCGAGATGCCGTGGCGCGGGGATTTTTTTATGCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCa                                     <  1:1258320/150‑1 (MQ=255)
                                                                                                                  gcatGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTAGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGc                                   <  1:327058/147‑1 (MQ=255)
                                                                                                                  gCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGc                                   <  1:568495/150‑1 (MQ=255)
                                                                                                                   ccTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCg                                  <  1:1873267/150‑1 (MQ=255)
                                                                                                                   ccTGCTACTGCTGATTTTTTCCCCGCGACATGCCGGGGCGGGGGGAATTTTTTATCCGGGCGTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACCTACACCCGCCTACCGGAACGGACCTACATGACCTGTTTGATCACGTGAGCg                                  >  2:1635640/1‑150 (MQ=255)
                                                                                                                            gctgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGATATCCATCACATGTTGCATCACTTCAGCGTCAGTTtgt                         >  1:677973/1‑150 (MQ=255)
                                                                                                                            gctgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTtgt                         >  1:359777/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)
                                                                                                                               gATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGtgct                      <  2:877250/150‑1 (MQ=255)
                                                                                                                               gATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGtgct                      >  1:1084425/1‑150 (MQ=255)
                                                                                                                                 aaaaatCCCCGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGtgctgc                    >  2:174928/6‑150 (MQ=255)
                                                                                                                                  tttttCCCCGCGACATGCCGTGGCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCa                   <  1:183350/150‑1 (MQ=255)
                                                                                                                                    tttCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATc                 <  1:1313763/150‑1 (MQ=255)
                                                                                                                                       ccccGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCttt              >  1:642024/1‑150 (MQ=255)
                                                                                                                                        cccGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTa             >  2:1292031/1‑150 (MQ=255)
                                                                                                                                        cccGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTa             <  1:870534/150‑1 (MQ=255)
                                                                                                                                         ccGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTaa            >  1:1655952/1‑150 (MQ=255)
                                                                                                                                          cgcgACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTaaa           >  2:1007759/1‑150 (MQ=255)
                                                                                                                                           gccaCATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAAt          <  2:155768/147‑1 (MQ=255)
                                                                                                                                            cgACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTTCATCACTTCCGCGTCAGTTTGTGCTGCATCTTTAAATg         >  1:284709/1‑150 (MQ=255)
                                                                                                                                               cctGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTGGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGatt      <  2:1175360/148‑1 (MQ=255)
                                                                                                                                                aTGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGattg     >  1:1504737/1‑150 (MQ=255)
                                                                                                                                                aTGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGattg     <  1:1714390/150‑1 (MQ=255)
                                                                                                                                                  gCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGattgca   >  2:1814326/1‑150 (MQ=255)
                                                                                                                                                   ccGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGattgcat  <  1:1556101/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCAT  >  NC_000913/1169235‑1169531

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: