Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,767,667 T→G 10.1% intergenic (+312/‑43) rnlB → / → yfjP CP4‑57 prophage; antitoxin RnlB/CP4‑57 prophage; putative GTP‑binding protein YfjP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,767,6670TG10.1% 132.1 / 9.7 62intergenic (+312/‑43)rnlB/yfjPCP4‑57 prophage; antitoxin RnlB/CP4‑57 prophage; putative GTP‑binding protein YfjP
Reads supporting (aligned to +/- strand):  ref base T (27/27);  new base G (3/3);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-02

TGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGATAAATCAGTTAACCC  >  NC_000913/2767525‑2767809
                                                                                                                                              |                                                                                                                                               
tGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGACCGGCAGAACCt                                                                                                                                                 >  2:737248/1‑143 (MQ=255)
    aCCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTcgcg                                                                                                                                             >  1:522795/1‑143 (MQ=255)
     ccAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGACAAAATCTTCAGATAAATATCCTTCTCGTGAAGAAACCAATGACCGGCAGAACCTCGCgg                                                                                                                                            >  2:81338/1‑143 (MQ=255)
            gcgggAAATCCCATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGcc                                                                                                                                     <  1:158493/140‑1 (MQ=255)
            cttGGAAGTCCTATGAGGCCTTGCCCGGGAAGTGGTGCTGTGAGTGCAGAAACATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATCGATATCATTCTCCTGAACAATCAGATTAACGGCAGAACCTCGCGGTTCTGcc                                                                                                                                     <  1:1013991/141‑1 (MQ=255)
             cTGGAAATACTATGATGCCTTTCCCGGGAAATGGTGCTGTGAGCTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCg                                                                                                                                    <  1:964733/143‑1 (MQ=255)
                gAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGACCAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGttt                                                                                                                                 >  1:45886/1‑143 (MQ=255)
                gAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGttt                                                                                                                                 >  1:529946/1‑143 (MQ=255)
                     ccTATGATGCCTTTTCCGGTAAATGGTGCTGTGAGTTCATTAATATTCTGCTATCTATGTGCATTGAAAAAATCGTCAGATAGAGATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGtt                                                                                                                           <  1:373187/143‑1 (MQ=255)
                      cTAGGATGCCTTTCCCTGTAAATGGTGCTGTAAGTTAAGAGAGATTCTGCAATGTATGTGCATTGAAAGAACCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGttcatcgga                                                                                                                           <  1:1019351/143‑8 (MQ=37)
                        atgatgCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTa                                                                                                                        <  1:240914/143‑1 (MQ=255)
                         tgatgCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGCTTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTAc                                                                                                                       >  1:749407/1‑143 (MQ=255)
                         tgatgCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTAc                                                                                                                       <  1:1075270/143‑1 (MQ=255)
                         tgatgCCTTTCCCGGTAAATGGTGCTGTCAGTTAAGAACGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTAc                                                                                                                       <  1:444663/143‑1 (MQ=255)
                           atgCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACtt                                                                                                                     >  1:850485/1‑143 (MQ=255)
                             gCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACtttt                                                                                                                   >  1:806560/1‑143 (MQ=255)
                               cTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGGAATATATGTGCATTGAAAAAATCTTTAGATAGATATCATTCTCGTGAACAATCCGATGAACGGCAGAACCCCGCGGTTCTGCCGTTGTTGTTTTACTTTTCa                                                                                                                 <  1:962092/143‑1 (MQ=255)
                                  tCCCGGGAAATTGGGCTGTGAGTTCGGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCagag                                                                                                              <  1:162850/143‑1 (MQ=255)
                                      ggTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGCTCTGACGTTTTTGTTTTACTTCTCAGAGGAgt                                                                                                          >  2:910325/1‑143 (MQ=255)
                                      ggTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAAACTATGTGCATTGAAAAAATATTAAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACAACGCGGTTCTGCCGTTATTGTTTTACTTTTCAGAAGAgt                                                                                                          >  2:99717/1‑143 (MQ=255)
                                           aTGGTGCGGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATGCTCCTCATCAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtga                                                                                                     <  2:649026/143‑1 (MQ=255)
                                            tGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgat                                                                                                    >  2:389144/1‑143 (MQ=255)
                                              gTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgatga                                                                                                  <  1:737248/143‑1 (MQ=255)
                                               tGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGGGATGaa                                                                                                 >  1:471614/1‑143 (MQ=255)
                                                gCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCATGAAGAATCCGATGAACGGAAGAACCTCGCGGTTCTGCCGTTTTCGTATTACGTGTCAGAGGAGTGGTGAGGaaa                                                                                                >  2:537241/1‑143 (MQ=255)
                                                  tgtgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGAGGTGATgaagaa                                                                                              >  2:896832/1‑143 (MQ=255)
                                                    tgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGGTATGAATAAtt                                                                                            >  1:329653/1‑143 (MQ=255)
                                                      aGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGTCAGAACCTAGCGGTTCTGCCGTTTTTGTTTTACTATTCAGAGGAGAGGTCACGAAACAgttt                                                                                          >  2:11020/1‑143 (MQ=255)
                                                      aGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAtttt                                                                                          >  2:191368/1‑143 (MQ=255)
                                                         tCATAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGaa                                                                                       <  2:1067408/143‑1 (MQ=255)
                                                         tCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGaa                                                                                       >  2:226805/1‑143 (MQ=255)
                                                               tgagtCTGTAATCTGTGTGCATTGGAATGATCTTCAGGTAGATATCATTCTCCTGATCAATCCGGTGAACGGCAGAGCCTCGCGGTTCTGCCGTTTTTGTGTTACTTTTCGGAGGAGTGGTGATGAACAATTTTGAAGTTCTg                                                                                 <  2:544624/139‑1 (MQ=255)
                                                                  ttCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGCGTGGTGATGAAAAATTTTGAAGTTCTGCAg                                                                              >  1:204757/1‑143 (MQ=255)
                                                                   tCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGc                                                                             <  1:615414/143‑1 (MQ=255)
                                                                    cTGCAATCTATGTGCATTGAAAAAATCTTCCGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGcc                                                                            <  2:63001/143‑1 (MQ=255)
                                                                       cAATCTATGTACACTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACt                                                                         <  1:761541/143‑1 (MQ=255)
                                                                          tCTTTGTGCATTGATAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACa                                                                      <  1:513798/143‑1 (MQ=255)
                                                                             aggtgCATGGTAAAACTCCTCGGATGGAAGTCACACTCCTGAGGAATCCGATGTACCGCATAACCTCGCGGGCCTGCCGTTCTTGTTTTACTCTTCAGAGGAGTGGTGATGAAAAATTTTGACGTTCTGCAGCCACTACAGaa                                                                   <  2:787888/141‑1 (MQ=255)
                                                                                tcaagtGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTc                                                                <  1:825651/138‑1 (MQ=255)
                                                                                 gCATTGAAAAAATCTTCAGATAGACATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTcc                                                               <  1:522332/143‑1 (MQ=255)
                                                                                  cATTGAAAAAATCTTCAGACAGACATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTccc                                                              <  1:62556/143‑1 (MQ=255)
                                                                                        aaaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCg                                                        <  1:818285/143‑1 (MQ=255)
                                                                                                    atagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTAt                                            <  1:239858/143‑1 (MQ=255)
                                                                                                         tatCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTCTCCGGTCTACCGCTATGGGTa                                       >  1:461280/1‑143 (MQ=255)
                                                                                                                ctcCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTCAAGTTCTGCAGCCACTACAGAACTACCTTTCCGTTCTACCGCTATGGGTATCTGAAc                                >  1:641094/1‑143 (MQ=255)
                                                                                                                 tcCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTACGGTCTACCGCACTGGGTATCTGAACg                               >  2:602047/1‑143 (MQ=255)
                                                                                                                   cTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTTAGTGTGCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACgc                              >  2:227539/1‑142 (MQ=255)
                                                                                                                   cTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGAAGTGGGGATGAAACACTCAGAAGATTTGCAGCCCCTACAGAACTCCCTTTCCGGTCTAGCGCTATGCGCATCAGAACgca                             >  2:641441/1‑141 (MQ=255)
                                                                                                                    tGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGtat                            >  2:836478/1‑143 (MQ=255)
                                                                                                                     gAACAGTCCGATGAACGGCAGGACCTCGCGGTTCTGACGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGATATATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGtata                           <  2:376570/143‑1 (MQ=255)
                                                                                                                     gAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGtata                           >  1:355930/1‑143 (MQ=255)
                                                                                                                     gAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGtaga                           >  1:286487/1‑141 (MQ=255)
                                                                                                                          aTCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACGACAGAACTCCCTTTCCGGCCTACCGCTATGAGTATCTGAACGTATATTgca                      <  2:558978/143‑1 (MQ=255)
                                                                                                                             cGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTgcagca                   <  1:511371/143‑1 (MQ=255)
                                                                                                                               aTGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGa                 >  1:766829/1‑143 (MQ=255)
                                                                                                                                  accGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGATaa              <  1:389144/141‑1 (MQ=255)
                                                                                                                                   aCGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATCGCAGCAGATaaa             >  1:352880/1‑143 (MQ=255)
                                                                                                                                      gCAGAACCTCGCGGTTCTGCGGTGTGTGTGTGACTGGTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGATAAATCa          <  1:413850/143‑1 (MQ=255)
                                                                                                                                       cAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGGTGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTCGGGGGAGCTGAGCGTATCTTGCAGCGGGTaaagcgg         >  1:880281/1‑139 (MQ=255)
                                                                                                                                       cAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCCGCAGATAACTCac         >  2:667889/1‑142 (MQ=255)
                                                                                                                                       cAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGATAAATCAg         >  2:911096/1‑143 (MQ=255)
                                                                                                                                        agcaCCTCGCGGTTCTGCCGTTTTTGTTTTTATTTTCAGCGGAGTGGTGATGACATTTTTTGCAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATTGGTAGCTGAACGTATATTGCAGCAGATAAATCAGt        <  2:522795/140‑1 (MQ=255)
                                                                                                                                             cTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGATAAATCAGTTAAcc   <  1:401978/143‑1 (MQ=255)
                                                                                                                                              tCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTCTGGGTATCTGAACGTATATTGCAGCAGATAAATCAGTTAAccc  >  1:515442/1‑143 (MQ=255)
                                                                                                                                              tCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGCTAAATCAGCTAAccc  >  1:249193/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                               
TGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGATAAATCAGTTAACCC  >  NC_000913/2767525‑2767809

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