Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 282,841 G→A P188P (CCG→CCA)  yagE → CP4‑6 prophage; putative 2‑keto‑3‑deoxygluconate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913282,8410GA99.7% 1181.8 / NA 291P188P (CCG→CCA) yagECP4‑6 prophage; putative 2‑keto‑3‑deoxygluconate aldolase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (144/146);  total (145/146)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 296 total aligned reads displayed.

CTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCG  >  NC_000913/282743‑282937
                                                                                                  |                                                                                                
    ccGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTc                                                                                            >  1:521167/1‑101 (MQ=255)
    ccGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTc                                                                                            <  1:520702/101‑1 (MQ=255)
      gACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCAc                                                                                          >  2:521171/1‑101 (MQ=255)
       aCTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCAcc                                                                                         >  2:521174/1‑101 (MQ=255)
        cTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCg                                                                                        <  2:520604/101‑1 (MQ=255)
         tCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGt                                                                                       >  1:521175/1‑101 (MQ=255)
           gcgcAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGc                                                                                     >  2:521181/1‑101 (MQ=255)
           gcgcAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGc                                                                                     >  1:521180/1‑101 (MQ=255)
            cgcAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGct                                                                                    >  1:521186/1‑101 (MQ=255)
             gcagcaATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGctc                                                                                   <  1:520522/101‑1 (MQ=255)
              cagcaATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGctct                                                                                  <  2:520694/101‑1 (MQ=255)
              cagcaATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGctct                                                                                  >  2:521189/1‑101 (MQ=255)
              cagcaATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGctct                                                                                  >  1:521190/1‑101 (MQ=255)
              cagcaATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGctct                                                                                  <  2:520733/101‑1 (MQ=255)
                gcaATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGc                                                                                <  2:520621/101‑1 (MQ=255)
                gcaATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGc                                                                                >  1:521194/1‑101 (MQ=255)
                 caATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCg                                                                               <  2:520838/101‑1 (MQ=255)
                 caATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCg                                                                               <  2:520743/101‑1 (MQ=255)
                 caATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCg                                                                               >  1:521197/1‑101 (MQ=255)
                 caATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCg                                                                               <  1:520828/101‑1 (MQ=255)
                  aaTATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCgg                                                                              <  2:520867/101‑1 (MQ=255)
                  aaTATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCgg                                                                              <  1:520831/101‑1 (MQ=255)
                   atatTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGc                                                                             >  1:521201/1‑101 (MQ=255)
                    tattatCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCt                                                                            <  1:520730/101‑1 (MQ=255)
                     attatCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTa                                                                           >  2:521204/1‑101 (MQ=255)
                      ttatCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTac                                                                          >  2:521209/1‑101 (MQ=255)
                      ttatCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTac                                                                          >  1:521207/1‑101 (MQ=255)
                       tatCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTacg                                                                         <  2:520824/101‑1 (MQ=255)
                        atCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGCTCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCt                                                                            >  2:521212/1‑97 (MQ=255)
                           ggCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTAcgacga                                                                     >  1:521215/1‑101 (MQ=255)
                            gCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGat                                                                    <  2:520581/101‑1 (MQ=255)
                             cATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGatc                                                                   >  2:521224/1‑101 (MQ=255)
                             cATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGatc                                                                   >  1:521221/1‑101 (MQ=255)
                              aTCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCAcc                                                                                         >  1:521228/1‑78 (MQ=255)
                               tCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGatcat                                                                 >  2:521229/1‑101 (MQ=255)
                                cAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGatcatc                                                                >  1:521233/1‑101 (MQ=255)
                                 aaaGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCt                                                               <  1:521145/101‑1 (MQ=255)
                                  aaGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTg                                                              <  2:520722/101‑1 (MQ=255)
                                   aGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGctct                                                                                  <  1:520737/80‑1 (MQ=255)
                                   aGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCg                                                                               >  2:521238/1‑83 (MQ=255)
                                   aGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGt                                                             <  1:520671/101‑1 (MQ=255)
                                   aGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGt                                                             >  1:521237/1‑101 (MQ=255)
                                        ccATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAAt                                                        >  2:521242/1‑101 (MQ=255)
                                         cATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATa                                                       >  1:521244/1‑101 (MQ=255)
                                         cATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATa                                                       <  1:520972/101‑1 (MQ=255)
                                         cATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATa                                                       <  2:520921/101‑1 (MQ=255)
                                          aTCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATAc                                                      <  1:520665/101‑1 (MQ=255)
                                          aTCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATAc                                                      <  2:520406/101‑1 (MQ=255)
                                          aTCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATAc                                                      >  1:521248/1‑101 (MQ=255)
                                           tCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATAcc                                                     >  2:521249/1‑101 (MQ=255)
                                            cGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATAccc                                                    >  1:521252/1‑101 (MQ=255)
                                             gACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGct                                                                                    <  2:520170/68‑1 (MQ=255)
                                              aCTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCctg                                                  <  2:520807/101‑1 (MQ=255)
                                               cTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGctct                                                                                  >  2:521257/1‑68 (MQ=255)
                                                tCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCctgct                                                <  1:520766/101‑1 (MQ=255)
                                                 ccGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCctgct                                                <  2:520796/100‑1 (MQ=255)
                                                  cGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCtgctgc                                              >  2:521260/1‑101 (MQ=255)
                                                  cGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCtgctgc                                              <  1:520565/101‑1 (MQ=255)
                                                   gTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTgctgct                                             >  2:521263/1‑101 (MQ=255)
                                                     cGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTAcgacga                                                                     <  1:520792/75‑1 (MQ=255)
                                                     cGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTcg                                           >  1:521264/1‑101 (MQ=255)
                                                      gCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTcgg                                          <  1:520847/101‑1 (MQ=255)
                                                      gCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTcgg                                          <  1:520498/101‑1 (MQ=255)
                                                       cccACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTcggc                                         <  1:520354/101‑1 (MQ=255)
                                                         cACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTcggcgg                                       >  2:521271/1‑101 (MQ=255)
                                                         cACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTcggcgg                                       <  1:520876/101‑1 (MQ=255)
                                                         cACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTcggcgg                                       >  1:521269/1‑101 (MQ=255)
                                                          aCCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCggcggc                                      <  2:520899/101‑1 (MQ=255)
                                                           ccTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGgcggcg                                     <  2:520636/101‑1 (MQ=255)
                                                           ccTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGgcggcg                                     <  1:520771/101‑1 (MQ=255)
                                                           ccTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGgcggcg                                     >  1:521273/1‑101 (MQ=255)
                                                            cTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTg                                                                                 <  1:520916/56‑1 (MQ=255)
                                                            cTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGa                                    <  1:520801/101‑1 (MQ=255)
                                                             tGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGAc                                   <  1:520754/101‑1 (MQ=255)
                                                              gcgcAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACg                                  <  2:521072/101‑1 (MQ=255)
                                                               cgcAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACgg                                 >  1:521279/1‑101 (MQ=255)
                                                               cgcAGCATGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACgg                                 >  2:521278/1‑101 (MQ=255)
                                                                gcagcaTGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATAcc                                                     <  2:521078/80‑1 (MQ=255)
                                                                  agcaTGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGc                              >  2:521284/1‑101 (MQ=255)
                                                                   gcaTGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCggcggc                                      >  1:521286/1‑92 (MQ=255)
                                                                     aTGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGat                           >  2:521289/1‑101 (MQ=255)
                                                                      tGATCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGata                          <  1:520894/101‑1 (MQ=255)
                                                                        aTCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATc                        >  1:521292/1‑101 (MQ=255)
                                                                        aTCCATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATc                        <  1:520684/101‑1 (MQ=255)
                                                                          ccATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCgg                      <  2:520770/101‑1 (MQ=255)
                                                                          ccATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCgg                      >  2:521294/1‑101 (MQ=255)
                                                                          ccATACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCgg                      <  2:520759/101‑1 (MQ=255)
                                                                            aTACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCg                    <  1:520862/101‑1 (MQ=255)
                                                                            aTACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCg                    >  2:521297/1‑101 (MQ=255)
                                                                             tACCGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGa                   <  2:520543/101‑1 (MQ=255)
                                                                               ccGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGc                 >  2:521301/1‑101 (MQ=255)
                                                                                cGTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCg                <  2:520644/101‑1 (MQ=255)
                                                                                 gTCAAAGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCgg               >  1:521303/1‑101 (MQ=255)
                                                                                    aaaGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGc                 <  1:521049/96‑1 (MQ=255)
                                                                                    aaaGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCaa            <  2:520998/101‑1 (MQ=255)
                                                                                     aaGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAAc           >  1:521305/1‑101 (MQ=255)
                                                                                      aGGTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACt          >  1:521308/1‑101 (MQ=255)
                                                                                       ggTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACtt         <  1:520707/101‑1 (MQ=255)
                                                                                        gTGCCCATCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACttt        <  1:520781/101‑1 (MQ=255)
                                                                                              aTCCACACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCg  >  2:521313/1‑101 (MQ=255)
                                                                                                  |                                                                                                
CTCGCCGACTCGCGCAGCAATATTATCGGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCG  >  NC_000913/282743‑282937

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