Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 254,228 A→T 17.5% intergenic (+26/+31) prfH → / ← pepD pseudogene, RF‑1 domain family; probable peptide chain release factor/aminoacyl‑histidine dipeptidase (peptidase D)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913254,2280AT17.5% 302.8 / 49.8 125intergenic (+26/+31)prfH/pepDpseudogene, RF‑1 domain family; probable peptide chain release factor/aminoacyl‑histidine dipeptidase (peptidase D)
Reads supporting (aligned to +/- strand):  ref base A (48/55);  new base T (10/12);  total (58/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 129 total aligned reads displayed.

AATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTG  >  NC_000913/254130‑254327
                                                                                                  |                                                                                                   
aaTGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGa                                                                                                     <  1:178581/99‑1 (MQ=255)
aaTGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGAt                                                                                                    <  2:1991792/100‑1 (MQ=255)
 aTGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGAtt                                                                                                   <  2:1501839/100‑1 (MQ=255)
  tGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTa                                                                                                  <  1:2399076/100‑1 (MQ=255)
    ttCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTaaa                                                                                                <  1:2480445/100‑1 (MQ=255)
    ttCCAACACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTaaa                                                                                                <  1:874392/100‑1 (MQ=255)
      ccATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACg                                                                                              <  1:864662/100‑1 (MQ=255)
       cATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGc                                                                                             >  2:2246102/1‑100 (MQ=255)
        aTCACGAGATTGTATGTTGTATCCCGCGACGTACATGTACATGGATGGCTTTTATCGAAGGATAATGAAGGTTGAAACTGCCGGAAGGCGATTAAACg                                                                                              <  1:493926/98‑1 (MQ=255)
           aCCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATc                                                                                         >  1:2082336/1‑100 (MQ=255)
              aGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCgg                                                                                      <  2:1804371/100‑1 (MQ=38)
               gATTGAGCGTGGCAACGCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCGgc                                                                                     <  2:1715251/100‑1 (MQ=37)
               gATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCGgc                                                                                     <  2:2182775/100‑1 (MQ=38)
               gATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCGgc                                                                                     <  1:433610/100‑1 (MQ=38)
               gATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgc                                                                                     >  1:420611/1‑100 (MQ=255)
                aTTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCGgca                                                                                    <  1:1766395/100‑1 (MQ=38)
                aTTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCGgca                                                                                    <  2:1780831/100‑1 (MQ=38)
                 ttGAACGTGGCAACCCGCGTCGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCGgcag                                                                                   <  1:1114391/100‑1 (MQ=38)
                 ttGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCGgcag                                                                                   <  1:783311/100‑1 (MQ=38)
                 ttGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCGgcag                                                                                   <  1:987085/100‑1 (MQ=38)
                 ttGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgcag                                                                                   >  2:411268/1‑100 (MQ=255)
                 ttGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgcag                                                                                   <  2:2082350/100‑1 (MQ=255)
                  tGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCGgcagt                                                                                  <  1:1003568/100‑2 (MQ=37)
                   gAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgcagca                                                                                 <  2:2455964/100‑1 (MQ=255)
                    aaCGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaa                                                                                >  2:1944729/1‑100 (MQ=255)
                        tGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaaatc                                                                             >  1:1374309/1‑99 (MQ=255)
                         ggCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCGgcagtttcatcc                                                                           <  2:111384/100‑9 (MQ=21)
                              cccGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAaataat                                                                      >  1:2434560/1‑100 (MQ=255)
                                  cgACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACt                                                                  >  1:242731/1‑100 (MQ=255)
                                    aCGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACtac                                                                >  1:265336/1‑98 (MQ=255)
                                     cGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCg                                                               >  2:2275062/1‑100 (MQ=255)
                                     cGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCg                                                               <  2:1641414/100‑1 (MQ=255)
                                      ggACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGc                                                              >  1:2285889/1‑100 (MQ=255)
                                      ggACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGc                                                              >  1:2185276/1‑100 (MQ=255)
                                      ggACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGc                                                              >  1:164322/1‑100 (MQ=255)
                                        ttccgatctgggATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCg                                                            <  2:1462248/91‑1 (MQ=255)
                                            ttACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAt                                                        <  1:2007411/100‑1 (MQ=255)
                                             tccAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCAGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtt                                                       >  1:570746/3‑100 (MQ=255)
                                             tACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtt                                                       >  1:550969/1‑100 (MQ=255)
                                             tACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtt                                                       >  1:326024/1‑100 (MQ=255)
                                              aCAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAttt                                                      >  2:2013523/1‑100 (MQ=255)
                                              aCAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAttt                                                      >  2:1556980/1‑100 (MQ=255)
                                               cAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtttc                                                     >  2:2524124/1‑100 (MQ=255)
                                               cAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAttt                                                      >  1:597243/1‑99 (MQ=255)
                                                 gggATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAGATCAAATAATTACTTCGCCGGAAtttctt                                                   <  1:2140294/100‑1 (MQ=255)
                                                 gggATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtttctt                                                   <  1:191221/100‑1 (MQ=255)
                                                 gggATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCAGGCAGCAAATCAAATAATTACTTCGCCGAAAtttct                                                    >  1:1045493/1‑99 (MQ=255)
                                                  ggATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtttcttt                                                  <  2:2334363/100‑1 (MQ=255)
                                                   gATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtttctttc                                                 <  1:1097346/100‑1 (MQ=255)
                                                      ggCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagc                                              <  2:1107017/100‑1 (MQ=255)
                                                      ggCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATAAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagc                                              <  1:1393686/100‑1 (MQ=255)
                                                       gCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagca                                             <  2:922257/100‑1 (MQ=255)
                                                       gCTTTTATCGAAGAATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagca                                             >  2:1740259/1‑100 (MQ=255)
                                                        cTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagcag                                            <  2:224726/100‑1 (MQ=255)
                                                        cTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagcag                                            >  2:2073575/1‑100 (MQ=255)
                                                         ttttATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagca                                             >  1:105108/1‑98 (MQ=255)
                                                           ttATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtc                                         <  2:2566526/100‑1 (MQ=255)
                                                             aTCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcag                                       <  2:59530/100‑1 (MQ=255)
                                                             aTCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcag                                       <  2:880310/100‑1 (MQ=255)
                                                              tCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagt                                      >  1:1118081/1‑100 (MQ=255)
                                                               cGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagtc                                     <  2:919645/100‑1 (MQ=255)
                                                                  aGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTcagc                                  <  2:1016553/100‑1 (MQ=255)
                                                                  aGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTcagc                                  <  2:828909/100‑1 (MQ=255)
                                                                   ggATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTcagca                                 <  2:567681/100‑1 (MQ=255)
                                                                        aTGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTc                            >  1:274824/1‑100 (MQ=255)
                                                                        aTGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTc                            <  1:1441342/100‑1 (MQ=255)
                                                                         tGAAGGATGAAACTGCCGGAAGGCGATTAAACACCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTcc                           <  2:1554630/100‑1 (MQ=255)
                                                                          gAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGACAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCa                          >  2:530347/1‑100 (MQ=255)
                                                                           ctggATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAg                         <  1:90854/98‑1 (MQ=255)
                                                                           aaGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAg                         >  1:2075375/1‑100 (MQ=255)
                                                                           aaGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAg                         <  1:679528/100‑1 (MQ=255)
                                                                           aaGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAg                         <  1:681664/100‑1 (MQ=255)
                                                                           aaGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGACAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAg                         <  1:530345/100‑1 (MQ=255)
                                                                             ggATGAAACTGCCTGAAGGCGATTAAACGCCATCCGGCGGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTa                       >  1:1595390/1‑100 (MQ=255)
                                                                             ggATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAgag                       >  2:90855/1‑98 (MQ=255)
                                                                             ggATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTa                       <  2:1595399/100‑1 (MQ=255)
                                                                                 gAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAAt                     >  2:1606362/1‑98 (MQ=255)
                                                                                 gAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGa                   <  2:1812584/100‑1 (MQ=255)
                                                                                   tgcTGCCGGATGGCGTTTAATCGCCTTCCGGCAGTAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGAcc                 >  1:1831735/3‑100 (MQ=25)
                                                                                   aaCTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGAcc                 <  2:1077034/100‑1 (MQ=255)
                                                                                    gcTGCCGGATGGCGTTTAATCGCCTTCCGGCAGTAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCt                >  2:2534319/2‑100 (MQ=37)
                                                                                    gcTGCCGGATGGCGTTTAATCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCt                >  2:1116883/2‑100 (MQ=37)
                                                                                    gcTGCCGGATGGCGTTTAATCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCt                >  2:2025667/2‑100 (MQ=37)
                                                                                     cTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTa               >  2:1476157/1‑100 (MQ=255)
                                                                                     cTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTa               <  2:456844/100‑1 (MQ=255)
                                                                                       gCCGGATGGCGTTTAATCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACg             >  2:2221583/1‑100 (MQ=38)
                                                                                       gCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCACATCATTACTTCGGCGGAATTTCTCTACGCACGTAAGTCAGCAGTGTCCAGCaaccgc                  >  1:2037007/1‑91 (MQ=255)
                                                                                       gCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACg             >  2:684875/1‑100 (MQ=255)
                                                                                        ccGGATGGCGTTTAATCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGc            >  1:1609539/1‑100 (MQ=38)
                                                                                        ccGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGc            >  1:1839914/1‑100 (MQ=255)
                                                                                            aaGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTc        >  2:2373714/1‑100 (MQ=255)
                                                                                              ggCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGa      <  1:467220/100‑1 (MQ=255)
                                                                                                cGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGAtg    >  2:548678/1‑100 (MQ=255)
                                                                                                cGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGAtg    <  2:1104931/100‑1 (MQ=255)
                                                                                                cGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGAtg    <  1:2073561/100‑1 (MQ=255)
                                                                                                cGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGAtg    <  1:2057134/100‑1 (MQ=255)
                                                                                                cGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGAGTACAGTAATGACATACgaat          >  1:1125211/1‑91 (MQ=255)
                                                                                                  aTTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCATTCAGCAGTGTCCAGTAATGACCTACGCTTTCGAtgtg  <  1:1379785/100‑1 (MQ=255)
                                                                                                  aTTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGAtgtg  <  2:597246/100‑1 (MQ=255)
                                                                                                  aTTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGAtgtg  <  2:2353390/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
AATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTG  >  NC_000913/254130‑254327

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: