Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,280,219 A→C 5.1% 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 5.1% 297.1 / 7.4 99E80D (GAA→GACyjcEputative transporter YjcE
Reads supporting (aligned to +/- strand):  ref base A (48/46);  new base C (2/3);  total (50/49)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.31e-01

CGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCTGTGG  >  NC_000913/4280071‑4280367
                                                                                                                                                    |                                                                                                                                                    
cgcttATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAAc                                                                                                                                                     >  1:1050008/1‑150 (MQ=255)
 gcttATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACa                                                                                                                                                    <  2:465099/150‑1 (MQ=255)
  cttATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACAt                                                                                                                                                   <  1:682006/150‑1 (MQ=255)
       gCAAATCGCCCTCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCg                                                                                                                                              <  1:1537394/150‑1 (MQ=255)
        cAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCga                                                                                                                                             >  2:1249873/1‑150 (MQ=255)
         aaaTCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGACCATGTTCgag                                                                                                                                            <  2:116691/150‑1 (MQ=255)
         aaaTCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCgag                                                                                                                                            <  1:625598/150‑1 (MQ=255)
            tCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCgagaga                                                                                                                                         <  1:510831/150‑1 (MQ=255)
            tCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCgagaga                                                                                                                                         <  1:396709/150‑1 (MQ=255)
             cGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGAt                                                                                                                                        >  2:1312989/1‑150 (MQ=255)
               ccATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGAttt                                                                                                                                      >  1:557233/1‑150 (MQ=255)
                 aTCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGACCATGTTCGAGAAAttcac                                                                                                                                    <  2:199671/150‑4 (MQ=255)
                  tCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGAAAATCTCTCGACCATGTTCGAGAAAttcacg                                                                                                                                   <  1:521225/150‑5 (MQ=255)
                   cGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCgg                                                                                                                                  >  1:739757/1‑150 (MQ=255)
                       gcgcTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTc                                                                                                                              >  2:445536/1‑150 (MQ=255)
                         gcTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCTGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTcgc                                                                                                                            >  2:800810/1‑150 (MQ=255)
                         gcTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGTTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTCGGCCTcgcg                                                                                                                           >  1:388012/1‑150 (MQ=255)
                              tGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCgcgctgg                                                                                                                       <  2:176772/150‑1 (MQ=255)
                               ggCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCgcgctggc                                                                                                                      >  2:1135051/1‑150 (MQ=255)
                                 cGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCgcgctggcgc                                                                                                                    <  2:114101/150‑1 (MQ=255)
                                   tGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCgcgctggcgctg                                                                                                                  >  2:637097/1‑150 (MQ=255)
                                     gCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCACGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCtggt                                                                                                                >  2:1275644/1‑150 (MQ=255)
                                      ccGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCtggtg                                                                                                               >  2:1506051/1‑150 (MQ=255)
                                       cGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCtggtgg                                                                                                              <  2:528352/150‑1 (MQ=255)
                                       cGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCtggtgg                                                                                                              <  1:615658/150‑1 (MQ=255)
                                         aCGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCtggtggtg                                                                                                            >  2:760428/1‑150 (MQ=255)
                                          cGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCtggtggtgg                                                                                                           <  1:1152994/150‑1 (MQ=255)
                                          cGTTTGGTGTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGGCGAGAGATTTTCGGCCTCGCGCTGGCGCtggtggtgg                                                                                                           <  2:806307/150‑1 (MQ=255)
                                           gtttggtttgCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGTCCTCGCGCTGGCGCTggtggtggt                                                                                                          >  2:710898/1‑150 (MQ=255)
                                              tggtttgCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCAc                                                                                                       <  1:1506051/150‑1 (MQ=255)
                                               ggtttgCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCAcc                                                                                                      >  2:1089506/1‑150 (MQ=255)
                                               ggtttgCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCAcc                                                                                                      <  1:414681/150‑1 (MQ=255)
                                                 tttgCATCTGGAGTTTGATCCTGAAATCATTTTAGTCTTATTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTAGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGt                                                                                                    >  2:1400843/1‑150 (MQ=255)
                                                   tgCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTgg                                                                                                  <  2:1076794/150‑1 (MQ=255)
                                                    gCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGt                                                                                                 >  2:708551/1‑150 (MQ=255)
                                                      aTGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCTGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCg                                                                                               <  1:800810/150‑1 (MQ=255)
                                                          ggAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCAt                                                                                           >  1:621656/1‑150 (MQ=255)
                                                          ggAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCAt                                                                                           >  1:464081/1‑150 (MQ=255)
                                                            aGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGTTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTgg                                                                                        <  2:388012/150‑1 (MQ=255)
                                                            aGTTTGATCCGTACCTCTTTTTTGTCTTGTTTATTCCTCCGTTGCTGTTCGCTGATTGCTGTAAAACGCCGACCCGTGAATTTCTCGAAAATTTTCTAGATATTTTCGTCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTg                                                                                         <  2:1395059/150‑1 (MQ=255)
                                                             gTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATTGTCTAGAGATTTTCGGCCTCGTGCTGTCGCTGTTGGTGGTGACCTTGGTCGGGATTgg                                                                                        <  2:373591/150‑1 (MQ=255)
                                                             gTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTgg                                                                                        >  1:678288/1‑150 (MQ=255)
                                                             gTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTgg                                                                                        >  1:268530/1‑150 (MQ=255)
                                                               ttGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCt                                                                                      >  1:458117/1‑150 (MQ=255)
                                                               ttGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCt                                                                                      <  1:546789/150‑1 (MQ=255)
                                                                    ccTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCtt                                                                                 <  1:637097/150‑1 (MQ=255)
                                                                     cTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTa                                                                                >  1:588424/1‑150 (MQ=255)
                                                                        aaCTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTc                                                                                                <  2:1190976/131‑1 (MQ=255)
                                                                        aaCTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTc                                                                                                >  1:1190976/1‑131 (MQ=255)
                                                                        aaCTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTAttt                                                                             <  2:866312/150‑1 (MQ=255)
                                                                         aCTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATCGCTGGAAAACGCCGACCCGTGAGTTTCTCGAACATGCTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTa                                                                            >  1:1129476/1‑150 (MQ=255)
                                                                                   gTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGgtggtg                                                                  <  2:464081/150‑1 (MQ=255)
                                                                                       tGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGATTCCTTATTTACTGGGTGGTGCCgg                                                              >  2:1007419/1‑150 (MQ=255)
                                                                                        gTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCggg                                                             <  2:173677/150‑1 (MQ=255)
                                                                                        gTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCggg                                                             <  1:484204/150‑1 (MQ=255)
                                                                                        gTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCggg                                                             <  1:1364042/150‑1 (MQ=255)
                                                                                         tttATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGGGCCGGGc                                                            <  2:1286122/150‑1 (MQ=255)
                                                                                          ttATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCa                                                           <  1:1523212/150‑1 (MQ=255)
                                                                                           tATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCAt                                                          >  2:967371/1‑150 (MQ=255)
                                                                                            aTCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCAtt                                                         >  2:188822/1‑150 (MQ=255)
                                                                                                cgccgTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGTTGTTTGTTGTCTCCGTGGTCGTCATTGTGTTACTTATTTACTGTGTGGTGCCGTGCATTCCGc                                                     <  2:1050008/150‑1 (MQ=255)
                                                                                                 gccgTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCt                                                    >  2:734004/1‑150 (MQ=255)
                                                                                                 gccgTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGGGAATTTCTCGAACCTGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGGGGGCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCt                                                    <  1:1249873/150‑1 (MQ=255)
                                                                                                       gCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATccc                                              >  1:1017010/1‑150 (MQ=255)
                                                                                                        cTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCt                                             >  1:1242778/1‑150 (MQ=255)
                                                                                                        cTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGCGCCGGGCATTCCGCTGATCCCt                                             <  1:709010/150‑1 (MQ=255)
                                                                                                          gTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGc                                           >  2:202543/1‑150 (MQ=255)
                                                                                                          gTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGc                                           <  1:1135051/150‑1 (MQ=255)
                                                                                                           ttCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGcc                                          <  1:1106346/150‑1 (MQ=255)
                                                                                                               cTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTg                                      >  2:1046650/1‑150 (MQ=255)
                                                                                                               cTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTg                                      >  2:1041318/1‑150 (MQ=255)
                                                                                                                     gCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTgg                                >  2:540480/1‑150 (MQ=255)
                                                                                                                     gCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTgg                                >  1:1022144/1‑150 (MQ=255)
                                                                                                                     gCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTgg                                <  2:1118128/150‑1 (MQ=255)
                                                                                                                      cTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTggc                               >  2:685401/1‑150 (MQ=255)
                                                                                                                      cTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTggc                               >  1:1150496/1‑150 (MQ=255)
                                                                                                                      cTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTggc                               <  1:943487/150‑1 (MQ=255)
                                                                                                                      cTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTggc                               >  1:929956/1‑150 (MQ=255)
                                                                                                                      cTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTAGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTggc                               >  1:1357444/1‑150 (MQ=255)
                                                                                                                       tGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTggcg                              <  2:108835/150‑1 (MQ=255)
                                                                                                                        ggAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGTGCTggcgg                             >  1:994634/1‑150 (MQ=255)
                                                                                                                        ggAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTggcgg                             >  2:1107458/1‑150 (MQ=255)
                                                                                                                            aaCGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCAATGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGt                         >  2:948843/1‑150 (MQ=255)
                                                                                                                             aCGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGTTGATCCCTGCCTTTGCGCTGGCGGCGGTg                        <  1:227758/150‑1 (MQ=255)
                                                                                                                             aCGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGGGGGCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTg                        <  2:1367968/150‑1 (MQ=255)
                                                                                                                              cGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGc                       <  2:588424/150‑1 (MQ=255)
                                                                                                                              cGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGc                       >  2:648442/1‑150 (MQ=255)
                                                                                                                               gCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCt                      <  2:557233/150‑1 (MQ=255)
                                                                                                                                ccGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCtt                     >  1:886228/1‑150 (MQ=255)
                                                                                                                                   aCCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGTGCTGGCGGCGGTGCTTtct                  <  2:994634/150‑1 (MQ=255)
                                                                                                                                       gTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTccga              >  2:628282/1‑150 (MQ=255)
                                                                                                                                       gTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTccga              >  1:690148/1‑150 (MQ=255)
                                                                                                                                         gAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCGGGCGCTGGGGGTGGGCCCCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTccgacc            <  1:445536/150‑1 (MQ=255)
                                                                                                                                           atctctcGACCATGTTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTccgaccga          >  1:613603/4‑150 (MQ=255)
                                                                                                                                           atctctcGACCATGTTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTccgaccga          >  1:722544/4‑150 (MQ=255)
                                                                                                                                              tctcGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGc       <  1:1089506/150‑1 (MQ=255)
                                                                                                                                                tcGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCtg     <  2:1017010/150‑1 (MQ=255)
                                                                                                                                                 cGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCtgt    <  2:1242778/150‑1 (MQ=255)
                                                                                                                                                  gAACATCGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGAGGTGGTCACCGTCGTCGGAATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCtgtg   >  2:1134856/1‑150 (MQ=255)
                                                                                                                                                   aaCATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCTGTgg  <  2:1022144/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCTGTGG  >  NC_000913/4280071‑4280367

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: