Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 281,581 G→A 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,5810GA97.9% 457.8 / ‑0.9 144R135W (CGG→TGG) yagACP4‑6 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (67/73);  total (69/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

Only 100 of 144 total aligned reads displayed.

TCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGG  >  NC_000913/281484‑281678
                                                                                                 |                                                                                                 
tcCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGGCGTGTTCGAACCCGc                                                                                                 >  1:123069/1‑100 (MQ=255)
  cagcagGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGccc                                                                                               >  2:2967148/1‑100 (MQ=255)
  cagcagGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGccc                                                                                               >  1:3310975/1‑100 (MQ=255)
   agcagGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCg                                                                                              <  2:1690807/100‑1 (MQ=255)
     cagGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTg                                                                                            >  2:1669994/1‑100 (MQ=255)
       gggTGGGCGCATGGCAGCGTCCACCGCCAAAAGGAGAGTGGCCCTGAAAGTCCATCTGCCAGAGGCGGGTCGGCGCGGCGTGGTCGAAGCAGCCCgggg                                                                                           >  1:2026139/1‑99 (MQ=255)
       gggTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCGTACAATCCATCTGCCAGAGGCGGTGCGGCGCGTCGGGTTCTAAACAGccc                                                                                               >  1:1081259/1‑95 (MQ=255)
         gTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCgg                                                                                        >  2:3106923/1‑100 (MQ=255)
         gTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCgg                                                                                        >  1:1597409/1‑100 (MQ=255)
         gTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCgg                                                                                        >  1:587257/1‑100 (MQ=255)
         gTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCgg                                                                                        >  1:1272797/1‑100 (MQ=255)
          tGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCggg                                                                                       <  1:2413820/100‑1 (MQ=255)
            aGCGGATGGCAGCGTCCACCGCCAAACGGAAAGTGGCCCTTAAAATCCATCGGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGaa                                                                                     <  1:1942794/100‑1 (MQ=255)
             gCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAAt                                                                                    >  2:558817/1‑100 (MQ=255)
             gCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGGTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAAt                                                                                    <  1:1441332/100‑1 (MQ=255)
                  tGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCg                                                                               <  2:3192924/100‑1 (MQ=255)
                   ggCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTTAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCg                                                                               <  1:881656/99‑1 (MQ=255)
                   ggCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCTAACCAGCCCGTGGCGGGAATGCCCgg                                                                              <  2:2121937/100‑1 (MQ=255)
                     cAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTg                                                                            <  1:1402513/100‑1 (MQ=255)
                     cAGCGTCCACCGCCAAAAGGAAAGTGGCCCCTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTg                                                                            >  1:2281481/1‑100 (MQ=255)
                      aGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGa                                                                           <  2:1544894/100‑1 (MQ=255)
                      aGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGGGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTg                                                                            <  1:2294003/99‑1 (MQ=255)
                       gCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGaa                                                                          <  1:2085078/100‑1 (MQ=255)
                        cGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAg                                                                         <  2:458217/100‑1 (MQ=255)
                         gTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAgc                                                                        >  2:594456/1‑100 (MQ=255)
                         gTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAgc                                                                        >  1:1093989/1‑100 (MQ=255)
                         gTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCACCCCGTGGCGGGAATGCCCGGTGAAgc                                                                        <  1:3139721/100‑1 (MQ=255)
                          tCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAgcg                                                                       >  1:1739679/1‑100 (MQ=255)
                           ccACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAgcgc                                                                      >  1:166214/1‑100 (MQ=255)
                            cACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGcc                                                                     <  1:3171648/100‑1 (MQ=255)
                            cACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGcc                                                                     <  2:1429373/100‑1 (MQ=255)
                             aCCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGccc                                                                    <  2:1128444/100‑1 (MQ=255)
                              ccGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCg                                                                   >  1:1528354/1‑100 (MQ=255)
                                gCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGgc                                                                 >  1:586583/1‑100 (MQ=255)
                                gCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGgc                                                                 >  2:2779744/1‑100 (MQ=255)
                                gCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGgc                                                                 >  2:3466798/1‑100 (MQ=255)
                                  cAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGgcag                                                               >  2:2137223/1‑100 (MQ=255)
                                   aaaaGGCAAGTGGCCCTTAAAATCCATCTGCCAGAAGCGGTTCGGCGCGTCGTGTTCGAAACAGCCCGTGGCGGGAATGCCCCGTGAAGCGCCCGgcagc                                                              >  1:1307024/1‑100 (MQ=255)
                                   aaaaGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGgcagc                                                              >  2:1548511/1‑100 (MQ=255)
                                    aaaGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTAGGCGCGTCGTGTTCGAACCAGCACGTGGCGGGAATGCCCGGTGAAGCGCCCGggagca                                                             >  1:2099938/1‑100 (MQ=255)
                                      aGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAgg                                                           >  2:1875819/1‑100 (MQ=255)
                                      aGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGAAGCAgg                                                           >  2:708666/1‑100 (MQ=255)
                                       ggAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGc                                                          >  2:1065412/1‑100 (MQ=255)
                                       ggAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGc                                                          <  1:628358/100‑1 (MQ=255)
                                       ggAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGc                                                          >  1:87580/1‑100 (MQ=255)
                                          aaGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCAt                                                       <  2:3041053/100‑1 (MQ=255)
                                            gTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATgg                                                     <  2:217195/100‑1 (MQ=255)
                                             tGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGc                                                    <  2:3315363/100‑1 (MQ=255)
                                             tGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGc                                                    <  1:1800699/100‑1 (MQ=255)
                                              ggCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCACCCGGCAGCAGGCCATGGCg                                                   >  1:1079345/1‑100 (MQ=255)
                                               gCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCgg                                                  <  2:438878/100‑1 (MQ=255)
                                               gCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCgg                                                  >  2:238623/1‑100 (MQ=255)
                                                 ccTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGc                                                <  2:3011715/100‑1 (MQ=255)
                                                  cTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGcc                                               <  1:1644302/100‑1 (MQ=255)
                                                      aaaTCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCa                                           >  2:3073698/1‑100 (MQ=255)
                                                      aaaTCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCa                                           <  1:674829/100‑1 (MQ=255)
                                                       aaTCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAg                                          >  1:2140903/1‑100 (MQ=255)
                                                         tCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGt                                        <  2:2379419/100‑1 (MQ=255)
                                                          ccATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGtt                                       <  1:1875792/100‑1 (MQ=255)
                                                           cATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCCGGCCATGGCGGGCCATCAGGtt                                       >  1:1327356/1‑99 (MQ=255)
                                                            aTCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTAt                                     >  1:898787/1‑100 (MQ=255)
                                                            aTCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTAt                                     <  2:166216/100‑1 (MQ=255)
                                                             tCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCCTGGCGGGCCATCAGGTTATg                                    <  1:2124551/100‑1 (MQ=255)
                                                             tCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATg                                    <  2:2194800/100‑1 (MQ=255)
                                                             tCTGCCAGAGGCGGGTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATg                                    <  1:3073648/100‑1 (MQ=255)
                                                              cTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATgg                                   <  2:1333111/100‑1 (MQ=255)
                                                               tGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGa                                  >  1:1678872/1‑100 (MQ=255)
                                                                gCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGAc                                 >  2:3439135/1‑100 (MQ=255)
                                                                  cAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACgg                               >  2:1803211/1‑100 (MQ=255)
                                                                  cAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACgg                               >  1:821359/1‑100 (MQ=255)
                                                                  cAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACgg                               <  2:1936371/100‑1 (MQ=255)
                                                                   agagGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGt                              <  1:2487893/100‑1 (MQ=255)
                                                                   agagGCGGTTCGGCGCGTCGTGTGCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACgg                               <  1:2859765/99‑1 (MQ=255)
                                                                      ggCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCt                           >  2:1251165/1‑100 (MQ=255)
                                                                      ggCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCt                           >  2:1459412/1‑100 (MQ=255)
                                                                      ggCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCt                           >  1:584017/1‑100 (MQ=255)
                                                                      ggCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCt                           >  1:2153914/1‑100 (MQ=255)
                                                                          gTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAg                       <  2:240323/100‑1 (MQ=255)
                                                                          gTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAg                       <  1:3371845/100‑1 (MQ=255)
                                                                             cGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGGAGGCCAGGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCg                    <  1:1495243/100‑1 (MQ=255)
                                                                             cGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCg                    <  1:558810/100‑1 (MQ=255)
                                                                              ggCGCGTCGTGTTCGAACCAGCGCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCgg                   >  1:2021193/1‑100 (MQ=255)
                                                                              ggCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCgg                   >  2:44709/1‑100 (MQ=255)
                                                                               gcgcGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCggg                  <  2:898799/100‑1 (MQ=255)
                                                                                 gcgTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCa                <  2:2845053/100‑1 (MQ=255)
                                                                                 gcgTCGTATTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCa                >  2:1830947/1‑100 (MQ=255)
                                                                                  cgtcgtGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCAt               >  1:1954/1‑100 (MQ=255)
                                                                                  cgtcgtGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCAt               <  1:426326/100‑1 (MQ=255)
                                                                                    tcgtGTTCGAACCGGCCCGTGTCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATgg             <  2:1155993/100‑1 (MQ=255)
                                                                                    tcgtGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATgg             <  2:1054568/100‑1 (MQ=255)
                                                                                     cgtGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGgt            <  2:1253445/100‑1 (MQ=255)
                                                                                      gtgtTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGgtg           <  2:3048265/100‑1 (MQ=255)
                                                                                      gtgtTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCAAGgtg           >  2:276051/1‑100 (MQ=255)
                                                                                       tgtTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGgtgt          >  2:939483/1‑100 (MQ=255)
                                                                                        gtTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGtgtg         >  1:2447640/1‑100 (MQ=255)
                                                                                        gtTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGtgtg         >  1:821710/1‑100 (MQ=255)
                                                                                         ttCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGc        <  1:1394641/100‑1 (MQ=255)
                                                                                          tCGAACCAGCCCGTGGCGGGAATGCCCGGGGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGc        <  1:1668429/99‑1 (MQ=255)
                                                                                            gAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGcccc     <  1:1188870/100‑1 (MQ=255)
                                                                                            gAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGcccc     <  2:976491/100‑1 (MQ=255)
                                                                                                 |                                                                                                 
TCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGG  >  NC_000913/281484‑281678

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