Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,403,148 A→C N371H (AAC→CAC)  hflK → modulator for HflB protease specific for phage lambda cII repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,403,1480AC99.6% 850.5 / NA 238N371H (AAC→CAC) hflKmodulator for HflB protease specific for phage lambda cII repressor
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (130/107);  total (130/108)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 239 total aligned reads displayed.

CCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGG  >  NC_000913/4403050‑4403244
                                                                                                  |                                                                                                
cccGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATaa                                                                                                 >  2:2974361/1‑100 (MQ=255)
  cGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACg                                                                                               >  1:2229725/1‑100 (MQ=255)
   gCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACgg                                                                                              >  2:690232/1‑100 (MQ=255)
   gCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACgg                                                                                              >  2:2918904/1‑100 (MQ=255)
   gCAAAGTGCTGGTTAACGATAAAGGTGGAAACCTGATGGTTCTGACGTTAGACCAGATGATGAAAGGTGGTCACGCCCCTGCGGCGAAAAGCGATCACgg                                                                                              >  1:1380622/1‑100 (MQ=255)
     aaaGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTg                                                                                            <  1:4922373/100‑1 (MQ=255)
     aaaGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTg                                                                                            <  1:2641874/100‑1 (MQ=255)
        gTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCa                                                                                         >  1:591842/1‑100 (MQ=255)
        gTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCa                                                                                         >  1:133538/1‑100 (MQ=255)
             ggTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAAt                                                                                    >  2:583599/1‑100 (MQ=255)
             ggTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAAt                                                                                    >  1:238843/1‑100 (MQ=255)
             ggTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAAt                                                                                    >  2:4715020/1‑100 (MQ=255)
             ggTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAAt                                                                                    >  1:373378/1‑100 (MQ=255)
             ggTTAACGATAAAGGTGGCAACCTGATGCTTCTGCCGTTAGACCAGAGGCTGAAAGGTGGTAACGCCCCTGCGCCGAAGAGCGGTCACGGTGCCAGCAAt                                                                                    >  2:4067944/1‑100 (MQ=255)
               ttAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATct                                                                                  >  1:4087751/1‑100 (MQ=255)
                naaCGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATctg                                                                                 >  1:309788/2‑100 (MQ=255)
                 aaCGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATctgc                                                                                <  2:62996/100‑1 (MQ=255)
                   cGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGCGCGATCACGGTGCCAGCAATctgctg                                                                              >  2:4479554/1‑100 (MQ=255)
                   cGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATctgctg                                                                              <  2:4994899/100‑1 (MQ=255)
                   cGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATctgctg                                                                              >  1:3826199/1‑100 (MQ=255)
                   cGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATctgctg                                                                              >  1:3776975/1‑100 (MQ=255)
                   cGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATctgctg                                                                              >  2:3738010/1‑100 (MQ=255)
                     aTAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCg                                                                            >  1:1603978/1‑100 (MQ=255)
                       aaaGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTc                                                                          >  1:713264/1‑100 (MQ=255)
                          ggTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTgc                                                                       >  1:152664/1‑100 (MQ=255)
                          ggTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTgc                                                                       >  1:881359/1‑100 (MQ=255)
                          ggTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTgc                                                                       >  1:800323/1‑100 (MQ=255)
                          ggTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTgc                                                                       >  1:5512721/1‑100 (MQ=255)
                           gTGGCAACCTGATGGTTCTGCCGTTAGTCCAGATGCTGAACGGTGGTAGCGCCTATGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTgac                                                                      >  2:5478080/1‑98 (MQ=255)
                           gTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTgcc                                                                      >  2:3159023/1‑100 (MQ=255)
                             ggCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAACCTGCTGCGTCTgccgc                                                                    >  1:65104/1‑100 (MQ=255)
                             ggCAACCAAAAGGATCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTTCGGCGAAGAGCGATCACGGTGCCCGCAATCTGCTGCGTCTgccgc                                                                    >  2:1637663/1‑100 (MQ=255)
                               cAAGCTGATGGTTCCGCCGATAGCGCAGATGCTGAAAGGTTGCAGCGCCTCTGCGGCGAAGAGCGATCACGGTACCAGCAACCTGCTGCGTCTACCGCCa                                                                  <  2:1015022/100‑1 (MQ=255)
                               cAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCa                                                                  >  1:1418243/1‑100 (MQ=255)
                                 aCCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGc                                                                >  2:5066129/1‑100 (MQ=255)
                                  ccTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGcc                                                               <  2:4965818/100‑1 (MQ=255)
                                   cTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCct                                                              >  2:760465/1‑100 (MQ=255)
                                   cTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCct                                                              >  1:1862630/1‑100 (MQ=255)
                                    tGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCctc                                                             >  2:5387392/1‑100 (MQ=255)
                                    tGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCctc                                                             <  2:3954976/100‑1 (MQ=255)
                                     gATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGACCCTGCGACGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCctct                                                            >  1:5286764/1‑100 (MQ=255)
                                      aTGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCtt                                                           <  1:3961686/100‑1 (MQ=255)
                                       tGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTtc                                                          >  1:3551398/1‑100 (MQ=255)
                                        ggTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTtcc                                                         <  1:4859130/100‑1 (MQ=255)
                                           tCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTtcctcc                                                      <  1:2767568/100‑1 (MQ=255)
                                            cTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCca                                                     >  2:5203498/1‑100 (MQ=255)
                                            cTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCca                                                     <  2:1180255/100‑1 (MQ=255)
                                            cTGCCGTTAGACCAGATGCTGAAAGGTGCTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCca                                                     <  1:5057808/100‑1 (MQ=255)
                                             tGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCCGCAATCTGCAGCGTCTGCTGCCAGCACCTtccccacc                                                    >  2:4594464/1‑95 (MQ=255)
                                             tGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCcan                                                    <  1:1749431/100‑2 (MQ=255)
                                             tGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCcac                                                    <  1:3565555/100‑1 (MQ=255)
                                             tGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCcac                                                    <  2:2392345/100‑1 (MQ=255)
                                              gCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCcaca                                                   <  1:205788/100‑1 (MQ=255)
                                                 gTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAAcc                                                <  1:473608/100‑1 (MQ=255)
                                                    aGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGt                                             >  2:96783/1‑100 (MQ=255)
                                                    aGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGt                                             >  2:4861798/1‑100 (MQ=255)
                                                    aGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGt                                             >  2:2792842/1‑100 (MQ=255)
                                                     gACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTg                                            >  1:5226394/1‑100 (MQ=255)
                                                       ccAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGa                                          <  2:2623800/100‑1 (MQ=255)
                                                        cAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGag                                         <  2:5437363/100‑1 (MQ=255)
                                                         aGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGagc                                        >  2:3691859/1‑100 (MQ=255)
                                                            tGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGagcaag                                     <  2:2181328/100‑1 (MQ=255)
                                                            tGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGagcaag                                     <  1:5275276/100‑1 (MQ=255)
                                                            tGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGagcaag                                     <  1:2996117/100‑1 (MQ=255)
                                                             gCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGagcaagc                                    <  2:2508348/100‑1 (MQ=255)
                                                             gCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGagcaagc                                    <  1:584166/100‑1 (MQ=255)
                                                              cTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGagcaagca                                   >  2:2003865/1‑100 (MQ=255)
                                                              cTGAAAGCTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCTAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTACTCCACAACCAGTGGagcaagta                                   <  2:164270/100‑3 (MQ=255)
                                                               tGAAAGGTGGTAACGCCCCTGCGGCGAAGAGAGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAgcaagcaa                                  >  1:4467962/1‑100 (MQ=255)
                                                                gAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAac                                 <  1:97847/100‑1 (MQ=255)
                                                                gAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAac                                 <  2:869782/100‑1 (MQ=255)
                                                                gAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAac                                 <  2:526848/100‑1 (MQ=255)
                                                                gAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAac                                 <  2:1408409/100‑1 (MQ=255)
                                                                gAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAac                                 <  1:2898203/100‑1 (MQ=255)
                                                                gAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAac                                 >  2:3221035/1‑100 (MQ=255)
                                                                 aaaGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAaca                                <  2:2284958/100‑1 (MQ=255)
                                                                 aaaGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAaca                                <  2:3113695/100‑1 (MQ=255)
                                                                  aaGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAacac                               >  1:4651323/1‑100 (MQ=255)
                                                                   aGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAacacn                              <  2:4178703/100‑2 (MQ=255)
                                                                   aGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACAcg                              <  1:937634/100‑1 (MQ=255)
                                                                    ggtggtAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACAcgt                             >  2:678750/1‑100 (MQ=255)
                                                                    ggtggtAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACAcgt                             <  1:4272125/100‑1 (MQ=255)
                                                                       ggtAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACAcgtcgt                          <  1:3182336/100‑1 (MQ=255)
                                                                        gtAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACgtcgtc                         <  2:3884126/100‑1 (MQ=255)
                                                                         tAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTcc                        <  2:1249720/100‑1 (MQ=255)
                                                                          aaCGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTcca                       <  1:5399770/100‑1 (MQ=255)
                                                                            cgcgCCTGCGGCGAAGATTGGTAGCTGTTCCAGCAATCTGCTGCTTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTccacc                     <  1:5563516/100‑1 (MQ=255)
                                                                                ccTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTc                 <  1:1899510/100‑1 (MQ=255)
                                                                                ccTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTc                 >  2:1933207/1‑100 (MQ=255)
                                                                                 cTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCa                >  2:3454029/1‑100 (MQ=255)
                                                                                   ggggcgAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTACTCCACAACCAGTGGAGCAAGCAACATGTCGTCCACGAGTCAgg              >  1:5115246/3‑100 (MQ=255)
                                                                                    cggcgAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAggg             <  1:4255132/100‑1 (MQ=255)
                                                                                        gAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGat         <  2:3552171/100‑1 (MQ=255)
                                                                                        gAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGat         >  1:4174583/1‑100 (MQ=255)
                                                                                           gagCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGAtatt      >  2:3614233/1‑100 (MQ=255)
                                                                                            agCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGAtatta     <  1:3346315/100‑1 (MQ=255)
                                                                                             gCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGAtattat    <  2:3335783/100‑1 (MQ=255)
                                                                                              cGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACCACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATg   >  2:1653384/1‑100 (MQ=255)
                                                                                              cGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATg   <  2:4852479/100‑1 (MQ=255)
                                                                                               gATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGCCGATATTATgg  >  1:4916785/1‑100 (MQ=255)
                                                                                                  |                                                                                                
CCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGG  >  NC_000913/4403050‑4403244

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