Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 281,581 G→A 24.2% R135W (CGG→TGG)  yagA ← CP4‑6 prophage; putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913281,5810GA24.2% 180.4 / 51.5 100R135W (CGG→TGG) yagACP4‑6 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (50/25);  new base A (10/14);  total (61/39)
Fisher's exact test for biased strand distribution p-value = 3.39e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

Only 100 of 101 total aligned reads displayed.

TCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTC  >  NC_000913/281484‑281680
                                                                                                 |                                                                                                   
tcCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGc                                                                                                   >  2:768890/1‑100 (MQ=255)
  ttccgatctGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGccc                                                                                                 <  1:139730/92‑1 (MQ=255)
  cagcagGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGccc                                                                                                 >  1:705898/1‑100 (MQ=255)
    ggagGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAAGCCATCTGCCAGAGGCGGGTCGGCGCGTCGTGTTCGAACCGGc                                                                                                   <  2:2124493/94‑1 (MQ=255)
    gcagGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGt                                                                                               >  1:12283/1‑100 (MQ=255)
       gggTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGc                                                                                            >  2:1134587/1‑100 (MQ=255)
       gggTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGc                                                                                            <  1:1626977/100‑1 (MQ=255)
         gTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCg                                                                                           >  1:1332856/1‑99 (MQ=255)
          tgtgCGGATGGCAGCGTCCACCGTCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCgg                                                                                          <  1:172155/96‑1 (MQ=255)
          tGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCggg                                                                                         >  2:23954/1‑100 (MQ=255)
           gAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGcccagatcggaa                                                                                        >  2:139727/1‑91 (MQ=255)
            aGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGaa                                                                                       >  1:424968/1‑100 (MQ=255)
              cGGATGGCCGCGTCCCCCGCCAAACGGACAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATg                                                                                     <  1:946269/100‑1 (MQ=255)
                gATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAAt                                                                                      >  2:1851978/1‑97 (MQ=255)
                   ggCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCg                                                                                 >  1:783392/1‑99 (MQ=255)
                    gCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGt                                                                               <  1:190213/100‑1 (MQ=255)
                    gCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGt                                                                               >  2:644295/1‑100 (MQ=255)
                    gCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAAGCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGGGGCGGGAATGCCCGGt                                                                               <  2:1080496/100‑1 (MQ=255)
                      aGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGa                                                                             <  1:1116283/100‑1 (MQ=255)
                        cGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAg                                                                           <  1:1250673/100‑1 (MQ=255)
                         gTCCACCGCCAAAAGGAAAGTGGCCATTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTTGCGGGAATGCCCGGCGAAg                                                                           >  2:175476/1‑99 (MQ=255)
                          tacACCGTCAAATGGAAAGTGGCCCTTAAAGTGCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGaa                                                                            <  1:2060870/95‑1 (MQ=255)
                          tCCACCGCCAAAATGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAgcg                                                                         <  2:2351883/100‑1 (MQ=255)
                          tCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAgcg                                                                         <  1:1134597/100‑1 (MQ=255)
                           ccACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAgcg                                                                         >  1:838945/1‑99 (MQ=255)
                           ccACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAgcgc                                                                        <  1:817955/100‑1 (MQ=255)
                            cACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAGTACCCGGTGAAgcgc                                                                        >  2:213433/1‑99 (MQ=255)
                            cACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCAGGTGAAGCGcc                                                                       >  2:786082/1‑100 (MQ=255)
                             aCCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGccc                                                                      <  1:2383793/100‑1 (MQ=255)
                               cGCCAAAAGGAAAGTTGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAcgcg                                                                         >  2:2047271/1‑95 (MQ=255)
                                gCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGgc                                                                   >  2:913325/1‑100 (MQ=255)
                                gCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGgc                                                                   >  2:1713488/1‑100 (MQ=255)
                                  cgaaaGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGgca                                                                  >  1:1614356/3‑99 (MQ=255)
                                  cAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCg                                                                     >  2:78481/1‑96 (MQ=255)
                                    aaaGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGgcagca                                                               >  2:2212882/1‑100 (MQ=255)
                                    aaaGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGgcagca                                                               <  2:1035187/100‑1 (MQ=255)
                                      aGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGcagcag                                                              <  1:741816/99‑1 (MQ=255)
                                          aaGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCAt                                                         >  1:166942/1‑100 (MQ=255)
                                          aaGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCAt                                                         <  2:746502/100‑1 (MQ=255)
                                          aaGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTAGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCAt                                                         >  2:2373859/1‑100 (MQ=255)
                                          aaGTGGCCCGTAAAATCCATCTGCCAGCGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCTTAGAGTGAATGCCCTGTGAAGCTCCCGGCAGCAGGcc                                                           <  2:728587/98‑1 (MQ=255)
                                           aGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGGGGCGGGAATGGCCGGTGAAGCGCCCGGCAGCAGGCCATg                                                        >  2:874036/1‑100 (MQ=255)
                                              ggCCCTTAAAATCCGTCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGc                                                      <  2:143823/99‑1 (MQ=255)
                                               gCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCgg                                                    <  2:2295741/100‑1 (MQ=255)
                                                cccTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCggg                                                   >  1:1027627/1‑100 (MQ=255)
                                                 ccTTAAACTCCATCTGCCAGAGGCGATTCGGCGCGTCGTGTTCGACCCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGc                                                  <  2:920374/100‑1 (MQ=255)
                                                 ccTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGc                                                  <  1:1535249/100‑1 (MQ=255)
                                                 ccTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGc                                                  >  2:789059/1‑100 (MQ=255)
                                                  cTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGc                                                  <  1:822059/99‑1 (MQ=255)
                                                  cTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAACCGCCCGGCAACAGGCCaaaacagac                                                  >  1:837404/1‑91 (MQ=255)
                                                    tAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCa                                                >  1:2275136/1‑99 (MQ=255)
                                                      aaaTCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATc                                              >  2:494531/1‑99 (MQ=255)
                                                      aaaTCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCa                                             >  1:1479868/1‑100 (MQ=255)
                                                      aaaTCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTAGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCa                                             >  2:1282814/1‑100 (MQ=255)
                                                       aaTCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAg                                            <  2:1336174/100‑1 (MQ=255)
                                                         tctaTCTGCCAGAGGCGGTTCAGCGCGTCGTTTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGcc                                                 <  1:2162351/90‑1 (MQ=255)
                                                         tCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGt                                          >  1:155057/1‑100 (MQ=255)
                                                          ccATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAAGCGGCCCGTGGCGGGAATGCGCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGt                                          <  2:1305239/99‑1 (MQ=255)
                                                          ccATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAgg                                           >  2:123976/1‑98 (MQ=255)
                                                           cATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACGGGCTCGTGGCGGGACTGCCCGGTGAAGCGCCCGGCAGCCGGCCAGGGCGGGCCATCAGGtt                                         >  1:972269/1‑99 (MQ=255)
                                                           cATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGTCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTa                                        <  1:1769587/100‑1 (MQ=255)
                                                           cATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGtt                                         >  2:1114889/1‑99 (MQ=255)
                                                           cATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTa                                        >  2:160865/1‑100 (MQ=255)
                                                              cTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATgg                                     >  2:504285/1‑100 (MQ=255)
                                                              cTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATgg                                     >  2:1040658/1‑100 (MQ=255)
                                                              cTGCCAGAGGCGGGTCGGCGCGTCGTGTTCGACATGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCTATCAGGTTATgg                                     <  2:2086951/100‑1 (MQ=255)
                                                                gCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGAc                                   >  1:493989/1‑100 (MQ=255)
                                                                  cAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACgg                                 >  1:1504666/1‑100 (MQ=255)
                                                                  cAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGACCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACgg                                 >  2:1680451/1‑100 (MQ=255)
                                                                     agGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCTGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGc                              >  1:1683037/1‑100 (MQ=255)
                                                                      ggCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCt                             >  2:2057453/1‑100 (MQ=255)
                                                                      ggCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCt                             >  1:1544571/1‑100 (MQ=255)
                                                                      ggCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCt                             >  2:1466268/1‑100 (MQ=255)
                                                                      ggCGGTTCGGCGCGTCGTGTTCGAAAAGGCTGGTGGAGGGAATGGCTGGTGAAGGGCGCGGCAGCAGGCCATGGAGGGGGATAAGGGTATGGAGGGTGc                              >  2:2044742/1‑99 (MQ=255)
                                                                         ggTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGaa                          >  1:1442326/1‑100 (MQ=255)
                                                                              ggCGCGTCGTGTTCGAACCCGCCCGTGGCCGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGCACGGTGCTGAAGGCg                      >  2:838672/1‑99 (MQ=255)
                                                                               gcgcGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCggn                    <  1:66534/100‑2 (MQ=255)
                                                                               gcgcGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCggg                    <  1:1108755/100‑1 (MQ=255)
                                                                                cgcgTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGc                   >  2:414239/1‑100 (MQ=255)
                                                                                cgcgTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGc                   <  1:1850482/100‑1 (MQ=255)
                                                                                    tcgtGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATgg               <  1:1710976/100‑1 (MQ=255)
                                                                                    tcgtGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATgg               <  2:790921/100‑1 (MQ=255)
                                                                                     cgtGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGgt              >  2:1769250/1‑100 (MQ=255)
                                                                                     cgtGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCACGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGgt              >  2:598525/1‑100 (MQ=255)
                                                                                      gtgtTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGgtg             >  2:403102/1‑100 (MQ=255)
                                                                                      gtgtTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGgtg             >  1:236493/1‑100 (MQ=255)
                                                                                      gtgtTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGgtg             <  1:1477412/100‑1 (MQ=255)
                                                                                       tgtTCGAACCGGCCCGTGGCGGGAATGCCCAGTAAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCAGGCATGgtg             >  2:1575567/1‑99 (MQ=255)
                                                                                          tCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGcc         <  1:902421/100‑1 (MQ=255)
                                                                                           cGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGccc        >  2:1890766/1‑100 (MQ=255)
                                                                                           cGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGGTATGGACGGTGCTGAAGACGGGCATGGGGTGccc        >  1:2036880/1‑100 (MQ=255)
                                                                                            gAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGcccc       <  2:54457/100‑1 (MQ=255)
                                                                                            gAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGcc         >  1:753842/1‑98 (MQ=255)
                                                                                             aaCCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCt      >  2:416319/1‑100 (MQ=255)
                                                                                              aCCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCt      >  1:2203158/1‑99 (MQ=255)
                                                                                               ccGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTgg    >  1:1223755/1‑100 (MQ=255)
                                                                                               ccGGCCCCTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTgg    >  1:1252911/1‑100 (MQ=255)
                                                                                                cGGGCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGt   <  2:1008888/100‑1 (MQ=255)
                                                                                                cGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTg     <  1:1513231/98‑1 (MQ=255)
                                                                                                 ggCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTc  >  1:1409967/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
TCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTC  >  NC_000913/281484‑281680

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