Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 254,223 A→T 18.3% intergenic (+21/+36) 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,2230AT18.3% 283.5 / 54.2 120intergenic (+21/+36)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/50);  new base T (10/12);  total (58/62)
Fisher's exact test for biased strand distribution p-value = 8.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

Only 100 of 122 total aligned reads displayed.

GGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTC  >  NC_000913/254126‑254321
                                                                                                 |                                                                                                  
ggcgAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAgg                                                                                                  <  1:2568565/100‑1 (MQ=255)
ggcgAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAgg                                                                                                  >  1:757787/1‑100 (MQ=255)
  cgAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCg                                                                                                >  2:210103/1‑100 (MQ=255)
  cgAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCg                                                                                                >  2:710892/1‑100 (MQ=255)
    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)
           cATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGc                                                                                       >  2:244178/1‑100 (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                                                                              <  2:1780831/100‑1 (MQ=38)
                     ttGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCGgcag                                                                             <  1:987085/100‑1 (MQ=38)
                     ttGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCGgcag                                                                             <  1:783311/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)
                         aCGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaa                                                                          >  1:922254/1‑99 (MQ=255)
                           gTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaaatc                                                                       >  1:1104927/1‑100 (MQ=255)
                            tGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaaatc                                                                       >  1:1374309/1‑99 (MQ=255)
                             ggCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCGgcagtttcatcc                                                                     <  2:111384/100‑9 (MQ=21)
                             ggCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaaatcaa                                                                     <  2:543174/100‑1 (MQ=255)
                               cAACCCGCGACGGACATTTACAGGGATGGCTTTTACCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaaatcaaat                                                                   <  2:906685/100‑1 (MQ=255)
                                  cccGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAaataat                                                                >  1:585553/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)
                                            aCATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCg                                                      <  2:2185291/100‑1 (MQ=255)
                                                 tccAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCAGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtt                                                 >  1:570746/3‑100 (MQ=255)
                                                 tACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCgataaaagc                                                 <  1:209479/100‑9 (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)
                                                     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:646388/100‑1 (MQ=255)
                                                          ggCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagc                                        <  2:1107017/100‑1 (MQ=255)
                                                          ggCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATAAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagc                                        <  1:1393686/100‑1 (MQ=255)
                                                           gCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagca                                       >  2:1897736/1‑100 (MQ=255)
                                                           gCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagca                                       <  2:922257/100‑1 (MQ=255)
                                                            cTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagcag                                      >  2:2316028/1‑100 (MQ=255)
                                                            cTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagcag                                      <  2:224726/100‑1 (MQ=255)
                                                            cTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGAGTAAACGCCATCCGGCAGCCAATCAAATAATTACTTCGCCGGAATTTCTTTcagcag                                      >  1:2389137/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)
                                                                  tCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagt                                >  1:1118081/1‑100 (MQ=255)
                                                                   cGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagtc                               <  2:919645/100‑1 (MQ=255)
                                                                    gAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagtca                              <  1:1671248/100‑1 (MQ=255)
                                                                      aGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTcagc                            <  2:828909/100‑1 (MQ=255)
                                                                       ggATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTcagca                           >  2:1639125/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:679528/100‑1 (MQ=255)
                                                                               aaGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAg                   >  1:2075375/1‑100 (MQ=255)
                                                                               aaGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAg                   <  2:2055929/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)
                                                                                     gAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAAt               >  2:1606362/1‑98 (MQ=255)
                                                                                     gAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGa             <  2:1812584/100‑1 (MQ=255)
                                                                                       tgcTGCCGGATGGCGTTTAATCGCCTTCCGGCAGTAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGAcc           >  1:1831735/3‑100 (MQ=25)
                                                                                       aaCTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGAcc           >  2:1777222/1‑100 (MQ=255)
                                                                                       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:456844/100‑1 (MQ=255)
                                                                                         cTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTa         >  2:1476157/1‑100 (MQ=255)
                                                                                          tGCCGGATGGCGTTTAATCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTAc        >  1:1318311/1‑100 (MQ=38)
                                                                                          tGCCGGATGGCGTTTAATCGCCTTCCGGCAGCAAATCAAATAAGTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGAAGTGTCCAGTAAGGACCTAc        >  1:1909142/1‑100 (MQ=37)
                                                                                           gCCGGATGGCGTTTAATCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACg       >  2:2221583/1‑100 (MQ=38)
                                                                                           gCCGGATGGCGTTTAATCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACg       >  1:398553/1‑100 (MQ=38)
                                                                                           gCCGGATGGCGTTTAATCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACg       >  1:74593/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)
                                                                                             cGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCt     <  1:2574680/100‑1 (MQ=255)
                                                                                             cGGAAGGCGATTAAACGACATCCGGCAGCACATCAAATAATTACTTCGCCGGAATTTCTTTCGGCAGTTCAGTCAGCAGTGTCCAGTAATTACCTACGc      <  1:477207/99‑1 (MQ=255)
                                                                                                aaGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTc  <  2:1099850/100‑1 (MQ=255)
                                                                                                 |                                                                                                  
GGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTC  >  NC_000913/254126‑254321

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