Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,280,219 A→C 8.8% E80D (GAA→GAC yjcE → putative transporter YjcE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,280,2190AC 8.8% 250.1 / 13.8 89E80D (GAA→GACyjcEputative transporter YjcE
Reads supporting (aligned to +/- strand):  ref base A (49/32);  new base C (5/3);  total (54/35)
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

CCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCTGTGGC  >  NC_000913/4280070‑4280368
                                                                                                                                                     |                                                                                                                                                     
ccgcttATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGaa                                                                                                                                                       >  2:1094957/1‑150 (MQ=255)
 cgcttATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAAc                                                                                                                                                      <  1:1562066/150‑1 (MQ=255)
 cgcttATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAAc                                                                                                                                                      <  2:1003661/150‑1 (MQ=255)
       tGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTc                                                                                                                                                >  2:1401366/1‑150 (MQ=255)
        gCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCg                                                                                                                                               >  1:12636/1‑150 (MQ=255)
         cAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTGTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCga                                                                                                                                              >  2:1485865/1‑150 (MQ=255)
             tCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCgagaga                                                                                                                                          >  2:353440/1‑150 (MQ=255)
               gCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGACCATGTTCGAGAAAtt                                                                                                                                        <  2:1014001/150‑1 (MQ=255)
               gCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGAtt                                                                                                                                        >  2:947549/1‑150 (MQ=255)
                  aTCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGGTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTc                                                                                                                                     <  2:415889/150‑1 (MQ=255)
                   tCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCg                                                                                                                                    <  1:363556/150‑1 (MQ=255)
                   tCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCg                                                                                                                                    <  2:1466532/150‑1 (MQ=255)
                    cGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGACCATGTTCGAGAAAttcacgg                                                                                                                                   <  1:747521/150‑6 (MQ=255)
                          gcTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTcgc                                                                                                                             >  1:1463631/1‑150 (MQ=255)
                          gcTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTcgc                                                                                                                             <  1:541815/150‑1 (MQ=255)
                            tACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCgcgc                                                                                                                           >  1:354103/1‑150 (MQ=255)
                            tACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCgcgc                                                                                                                           >  1:1322657/1‑150 (MQ=255)
                                ggCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCgcgctggc                                                                                                                       >  1:575675/1‑150 (MQ=255)
                                      gCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCtggt                                                                                                                 >  1:1481661/1‑150 (MQ=255)
                                      gCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCtggt                                                                                                                 <  1:792398/150‑1 (MQ=255)
                                       ccGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCtggtg                                                                                                                >  1:1024241/1‑150 (MQ=255)
                                           cGTTTGGTTTGCATGTGGATTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCtggtggtgg                                                                                                            <  1:1478655/150‑1 (MQ=255)
                                           cGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCtggtggtgg                                                                                                            <  1:1468370/150‑1 (MQ=255)
                                            gtttggtttgCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTggtggtggt                                                                                                           >  2:1301309/1‑150 (MQ=255)
                                                     gCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGt                                                                                                  >  2:828672/1‑150 (MQ=255)
                                                         gtgGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGc                                                                                              >  1:1195257/1‑150 (MQ=255)
                                                         gtgGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGCCGCTCGTGGTGCTCACCCTGGTCGGc                                                                                              >  2:1455533/1‑150 (MQ=255)
                                                         gtgGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCCAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGc                                                                                              >  1:912907/1‑150 (MQ=255)
                                                          tgGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCa                                                                                             >  2:1226888/1‑150 (MQ=255)
                                                           ggAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCAt                                                                                            <  2:247690/150‑1 (MQ=255)
                                                           ggAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCAt                                                                                            >  2:269931/1‑150 (MQ=255)
                                                             aGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTg                                                                                          <  1:1197035/150‑1 (MQ=255)
                                                                    tCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCt                                                                                   >  1:450466/1‑150 (MQ=255)
                                                                     ccTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCtt                                                                                  >  2:1368917/1‑150 (MQ=255)
                                                                       tGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTAt                                                                                <  1:559824/150‑1 (MQ=255)
                                                                       tGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTAt                                                                                >  2:606650/1‑150 (MQ=255)
                                                                        gAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTAtt                                                                               <  1:838233/150‑1 (MQ=255)
                                                                         aaCTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGGCACCGTGGTCGGCATTGGCTTCCTTAttt                                                                              <  1:1588891/150‑1 (MQ=255)
                                                                             ctTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTg                                                                          <  1:79435/150‑1 (MQ=255)
                                                                                tttAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGggt                                                                       <  1:1094957/150‑1 (MQ=255)
                                                                                tttAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGACTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGggt                                                                       >  1:1264338/1‑150 (MQ=255)
                                                                                 ttAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGggtg                                                                      <  1:1401366/150‑1 (MQ=255)
                                                                                     tCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGGCGGCATTGGCTTCCTTATTTACTGGGTGGTGc                                                                  <  1:387845/150‑1 (MQ=255)
                                                                                      cTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATATACTGGGTGGTGcc                                                                 >  1:635653/1‑150 (MQ=255)
                                                                                       ttGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCg                                                                >  1:1348114/1‑150 (MQ=255)
                                                                                        tGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCgg                                                               >  1:1564394/1‑150 (MQ=255)
                                                                                        tGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCgg                                                               <  1:490119/150‑1 (MQ=255)
                                                                                         gTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCggg                                                              <  2:12636/150‑1 (MQ=255)
                                                                                         gTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCggg                                                              <  1:245179/150‑1 (MQ=255)
                                                                                          tttATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGc                                                             <  1:450300/150‑1 (MQ=255)
                                                                                             aTCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCAtt                                                          >  1:164625/1‑150 (MQ=255)
                                                                                               cccGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTcc                                                        >  2:341218/1‑150 (MQ=255)
                                                                                               cccGCCGGTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGGCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGGGCCGGGCATTcc                                                        <  1:766131/150‑1 (MQ=255)
                                                                                                  gccgTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCt                                                     <  2:1481661/150‑1 (MQ=255)
                                                                                                   ccgTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTg                                                    <  1:1301309/150‑1 (MQ=255)
                                                                                                        gCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGCTccc                                               >  1:1440587/1‑150 (MQ=255)
                                                                                                        gCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATccc                                               >  2:1475738/1‑150 (MQ=255)
                                                                                                         cTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCt                                              >  1:1275851/1‑150 (MQ=255)
                                                                                                           gTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGAGTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGc                                            >  1:833148/1‑150 (MQ=255)
                                                                                                             tCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCt                                          >  1:70946/1‑150 (MQ=255)
                                                                                                              cGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCtt                                         >  2:128996/1‑150 (MQ=255)
                                                                                                                  gATGGCTGGAAACCGCCGACCCGGGAATTTCTCGACCATGGGCGAGAGATTTTCGGCCTCGCGCGGGCGCTGGTGGGGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGGGGTGCCGGGCATTCCGGGGATCCCTGCCTTTgcg                                     <  2:912907/150‑1 (MQ=255)
                                                                                                                  gATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTgcg                                     <  2:1024241/150‑1 (MQ=255)
                                                                                                                  gATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTgca                                     >  2:902092/1‑149 (MQ=255)
                                                                                                                    tGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTcc                                                        <  1:1116593/129‑1 (MQ=255)
                                                                                                                    tGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTcc                                                        >  2:1116593/1‑129 (MQ=255)
                                                                                                                       cTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTggc                                <  2:1195257/150‑1 (MQ=255)
                                                                                                                        tGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTggcg                               <  1:1368917/150‑1 (MQ=255)
                                                                                                                          gAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATATTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATAGGCTTCCTTATTTACTGGGAGGTGCCGGGCATTCCGCTGATCCCTGCCTTAGCGCTggcggc                             >  1:1589268/1‑150 (MQ=255)
                                                                                                                           aaaaCGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGCTCCCTGCCTTTGCGCTGgcggcg                            >  1:588423/1‑150 (MQ=255)
                                                                                                                           aaaaCGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGgcggcg                            >  2:716229/1‑150 (MQ=255)
                                                                                                                           aaaaCGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGgcggcg                            >  1:744798/1‑150 (MQ=255)
                                                                                                                               cGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGc                        >  1:152150/1‑150 (MQ=255)
                                                                                                                               cGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGc                        >  2:1121093/1‑150 (MQ=255)
                                                                                                                               cGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGc                        >  2:146950/1‑150 (MQ=255)
                                                                                                                                gCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCt                       >  1:219156/1‑150 (MQ=255)
                                                                                                                                    aCCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTtct                   <  1:1370062/150‑1 (MQ=255)
                                                                                                                                     cccGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTtctc                  >  2:1331353/1‑150 (MQ=255)
                                                                                                                                      ccgaaaaTCTCTCGCCCATGTTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTcc                 >  1:943805/6‑150 (MQ=255)
                                                                                                                                      ccgaaaaTCTCTCGACCATGTTCGAGAGATTTACGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTcc                 >  1:188944/6‑150 (MQ=255)
                                                                                                                                      ccgaaaaTCTCTCGACCATGTTCGAGAAATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTcc                 >  1:1178036/6‑150 (MQ=255)
                                                                                                                                       cGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTccg                >  1:1234094/1‑150 (MQ=255)
                                                                                                                                       cGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTccg                >  2:1341290/1‑150 (MQ=255)
                                                                                                                                            atctctcGACCATGTTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTccgaccga           >  2:934130/4‑150 (MQ=255)
                                                                                                                                            atctctcGACCATGTTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCACTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTccgaccga           >  2:518770/4‑150 (MQ=255)
                                                                                                                                             tttCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGAt          >  1:774612/1‑150 (MQ=255)
                                                                                                                                               tctcGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGc        >  1:380873/1‑150 (MQ=255)
                                                                                                                                                ctcGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCt       <  2:1463631/150‑1 (MQ=255)
                                                                                                                                                  cGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCtgt     >  1:126190/1‑150 (MQ=255)
                                                                                                                                                  cGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCGTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCtgt     <  2:1162500/150‑1 (MQ=255)
                                                                                                                                                    aaCATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCTGTgg   <  1:341218/150‑1 (MQ=255)
                                                                                                                                                     aCATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCTGTGgc  >  1:1229815/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCTGTGGC  >  NC_000913/4280070‑4280368

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