Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,776,109 (G)5→4 coding (191/459 nt) yfjX → CP4‑57 prophage; putative antirestriction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,776,1050G.97.6% 556.1 / 8.1 124coding (187/459 nt)yfjXCP4‑57 prophage; putative antirestriction protein
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base . (59/61);  total (61/63)
Fisher's exact test for biased strand distribution p-value = 6.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 141 total aligned reads displayed.

ATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATG  >  NC_000913/2776006‑2776200
                                                                                                   |                                                                                               
atCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCg                                                                                                 >  2:477064/1‑100 (MQ=255)
atCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCg                                                                                                 >  1:347037/1‑100 (MQ=255)
agccGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCg                                                                                                 >  1:2290444/3‑100 (MQ=255)
 tCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCgg                                                                                                <  1:3141639/100‑1 (MQ=255)
  ccGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCggg                                                                                               >  2:3033509/1‑100 (MQ=255)
  acGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGGGGATAACGCGGGAACCGCGCAGGGTCGCCTGGACGGACCGCGTCTGTGATGAGTACTGCgg                                                                                                >  2:1987258/2‑99 (MQ=255)
   cGCTTCTGGCCGCAGCACTTTGGCACCATCCCGGAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGACGGACCGCTTCTGTGATGAGTACTGCgggg                                                                                              >  2:2047850/1‑100 (MQ=255)
     cTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCggggta                                                                                            <  2:1803414/100‑3 (MQ=255)
     cTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCggggta                                                                                            <  1:3261349/100‑3 (MQ=255)
     cTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCggggta                                                                                            >  2:1314787/1‑98 (MQ=255)
     cTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTCTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCggggta                                                                                            <  1:352555/100‑3 (MQ=255)
       tCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCggggtatc                                                                                          <  1:2148623/100‑5 (MQ=255)
       tCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCggggtatc                                                                                          >  2:2465558/1‑96 (MQ=255)
        cTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCt                                                                                        >  1:2396923/1‑100 (MQ=255)
          ggCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTgg                                                                                      >  1:1836620/1‑100 (MQ=255)
            ccGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTc                                                                                    >  1:2858574/1‑100 (MQ=255)
            ccGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTc                                                                                    >  1:2857129/1‑100 (MQ=255)
             cGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTcc                                                                                   <  1:65763/100‑1 (MQ=255)
             cGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGCACTGAGGGGTATCTGGTct                                                                                   <  2:281884/100‑2 (MQ=37)
              gcagcaCTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCt                                                                                  >  1:2368535/1‑100 (MQ=255)
                 gcaCTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCtttt                                                                               <  1:1715901/100‑1 (MQ=255)
                 gcaCTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCtttt                                                                               <  2:2145840/100‑1 (MQ=255)
                  caCTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTa                                                                              <  1:3380383/100‑1 (MQ=255)
                  caCTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTa                                                                              <  1:743990/100‑1 (MQ=255)
                  caCTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTa                                                                              >  1:1519669/1‑100 (MQ=255)
                  caCTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTa                                                                              <  2:347036/100‑1 (MQ=255)
                  caCTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTa                                                                              <  1:615310/100‑1 (MQ=255)
                  caCTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTAGATGGACCGCTTCTGAGATGAGTACTGCGGGGTATCTGGTCCTTTTa                                                                              >  2:2292355/1‑100 (MQ=39)
                   aCTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTac                                                                             >  2:826039/1‑100 (MQ=255)
                    cTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTaca                                                                            <  2:2857118/100‑1 (MQ=255)
                    cTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTaca                                                                            <  2:2858563/100‑1 (MQ=255)
                       tGGCACCATCCCGCAGTGTATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGc                                                                         <  2:2643760/100‑1 (MQ=255)
                        ggCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCt                                                                        >  2:3194561/1‑100 (MQ=255)
                         gCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTc                                                                       <  1:2872262/100‑1 (MQ=255)
                         gCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTc                                                                       <  2:1603094/100‑1 (MQ=255)
                         gCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTc                                                                       >  1:2588338/1‑100 (MQ=255)
                         gCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTc                                                                       >  1:1913410/1‑100 (MQ=255)
                         gCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTc                                                                       <  1:511127/100‑1 (MQ=255)
                          cACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCa                                                                      >  1:1531807/1‑100 (MQ=255)
                          cACCATCCCGCAGAGGATAACGCTGGAACCGCGCATTTGCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTc                                                                       >  2:2106155/1‑100 (MQ=255)
                             cATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCa                                                                   >  2:1783822/1‑100 (MQ=255)
                             cATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCa                                                                   >  2:1594764/1‑100 (MQ=255)
                             cATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCa                                                                   <  1:1783826/100‑1 (MQ=255)
                              aTCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCaa                                                                  >  1:96005/1‑100 (MQ=255)
                               tCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCGGGTCCTTTTACACGCTCAGCAAt                                                                 <  2:3010658/100‑1 (MQ=255)
                                 ccGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATgg                                                               <  2:1519666/100‑1 (MQ=255)
                                  cGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGc                                                              >  1:727762/1‑100 (MQ=255)
                                  cGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGc                                                              <  1:1139345/100‑1 (MQ=255)
                                  cGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGc                                                              <  1:532943/100‑1 (MQ=255)
                                   gCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCg                                                             <  1:799318/100‑1 (MQ=255)
                                   gCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCg                                                             >  2:1227870/1‑100 (MQ=255)
                                      gTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGa                                                           >  2:2346278/1‑99 (MQ=255)
                                      gTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAg                                                          <  2:727762/100‑1 (MQ=255)
                                      gTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAg                                                          >  1:1565732/1‑100 (MQ=255)
                                       tGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGc                                                         >  1:148770/1‑100 (MQ=255)
                                        ggATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCg                                                        >  1:3039414/1‑100 (MQ=255)
                                           tAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGttt                                                     >  2:3360827/1‑100 (MQ=255)
                                            aaCGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGttt                                                     <  1:997936/100‑1 (MQ=255)
                                             aagCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTAt                                                   >  1:2291098/3‑100 (MQ=255)
                                             aCGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTAt                                                   <  2:96006/100‑1 (MQ=255)
                                             aCGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTAt                                                   >  2:2120762/1‑100 (MQ=255)
                                             aCGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTAt                                                   >  1:3140218/1‑100 (MQ=255)
                                              cGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCGATGGCGGAGCGTTTAt                                                   >  2:1948419/1‑99 (MQ=255)
                                              cGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATg                                                  <  2:3140206/100‑1 (MQ=255)
                                              cGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATg                                                  <  2:1913406/100‑1 (MQ=255)
                                              cGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATg                                                  <  2:1149891/100‑1 (MQ=255)
                                               gCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATgg                                                 >  1:3253577/1‑100 (MQ=255)
                                                 tGGAACCGCGCATTTTCGCCTGGATAGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCt                                               <  1:1314790/100‑1 (MQ=255)
                                                      ccGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGa                                          <  1:2013187/100‑1 (MQ=255)
                                                        gcgcATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATg                                        <  1:2446787/100‑1 (MQ=255)
                                                         cgcATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGc                                       <  2:2290438/100‑1 (MQ=255)
                                                          ccgatctCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCt                                      <  2:1622/94‑1 (MQ=37)
                                                                tCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGt                                <  2:2842386/100‑1 (MQ=255)
                                                                 cGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCtagatcgg                               >  1:1622/1‑93 (MQ=37)
                                                                 cGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTg                               >  2:2022274/1‑100 (MQ=255)
                                                                  gCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGa                              >  1:3022263/1‑100 (MQ=255)
                                                                   ccTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGAc                             <  1:1973510/100‑1 (MQ=255)
                                                                      ggatggaCCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGAt                          >  1:293565/1‑100 (MQ=255)
                                                                        atggaCCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGACGGTTTATGGCTCCTGATGCTGACGGTGACGATa                         >  1:2703970/1‑99 (MQ=37)
                                                                          ggaCCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAAt                      <  2:1565729/100‑1 (MQ=255)
                                                                          ggaCCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAAt                      <  1:2022278/100‑1 (MQ=255)
                                                                           gaCCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATg                     <  1:881714/100‑1 (MQ=255)
                                                                           gaCCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATg                     >  2:1589829/1‑100 (MQ=255)
                                                                           gaCCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATg                     >  2:557558/1‑100 (MQ=255)
                                                                            aCCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATgg                    <  1:3192104/100‑1 (MQ=255)
                                                                            aCCGCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATgg                    >  1:1386660/1‑100 (MQ=255)
                                                                               gCTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGGTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCAt                 >  2:2997162/1‑100 (MQ=255)
                                                                                cTTCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCAtt                >  2:580755/1‑100 (MQ=255)
                                                                                  tgaaggatgaGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTg              <  1:447305/95‑1 (MQ=255)
                                                                                  tCTGTGATGAGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTg              >  1:2036021/1‑100 (MQ=255)
                                                                                       gatgaGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCaa         >  1:1115469/1‑100 (MQ=255)
                                                                                         tgaGTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACg       <  2:247892/100‑1 (MQ=255)
                                                                                            gTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGc     >  2:2416646/1‑100 (MQ=255)
                                                                                            gTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCa    >  1:1170529/1‑100 (MQ=255)
                                                                                            gTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCa    >  1:2829832/1‑100 (MQ=255)
                                                                                            gTACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCa    <  1:302022/100‑1 (MQ=255)
                                                                                             tACTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGTCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCAt   <  2:2396916/100‑1 (MQ=255)
                                                                                              aCTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCAt   <  1:557559/99‑1 (MQ=255)
                                                                                              aCTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATg  <  1:221658/100‑1 (MQ=255)
                                                                                              aCTGCGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATg  <  2:87263/100‑1 (MQ=255)
                                                                                                   |                                                                                               
ATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATG  >  NC_000913/2776006‑2776200

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