Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 254,233 A→T 16.7% intergenic (+31/+26) 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,2330AT16.7% 323.1 / 52.5 133intergenic (+31/+26)prfH/pepDpseudogene, RF‑1 domain family; probable peptide chain release factor/aminoacyl‑histidine dipeptidase (peptidase D)
Reads supporting (aligned to +/- strand):  ref base A (54/55);  new base T (10/12);  total (65/68)
Fisher's exact test for biased strand distribution p-value = 8.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.31e-01

Only 100 of 141 total aligned reads displayed.

TTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACT  >  NC_000913/254134‑254331
                                                                                                   |                                                                                                  
ttCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTaaa                                                                                                    <  2:685881/100‑1 (MQ=255)
ttCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTaaa                                                                                                    >  1:1366951/1‑100 (MQ=255)
 tCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAAc                                                                                                   >  2:833328/1‑100 (MQ=255)
  ccATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACg                                                                                                  >  1:1011325/1‑100 (MQ=255)
  ccATCACCAGATTGAACGTGGCAAACCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACg                                                                                                  <  1:1864323/100‑1 (MQ=255)
   cATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGc                                                                                                 <  2:2172242/100‑1 (MQ=255)
    aTCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGc                                                                                                 <  2:876113/99‑1 (MQ=255)
       aCCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATc                                                                                             >  1:1921365/1‑100 (MQ=255)
       aCCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATc                                                                                             >  2:196015/1‑100 (MQ=255)
       aCCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATc                                                                                             >  2:455885/1‑100 (MQ=255)
       aCCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATc                                                                                             <  1:2249175/100‑1 (MQ=255)
       aCCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATc                                                                                             >  2:1594634/1‑100 (MQ=255)
          aGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCgg                                                                                          <  1:1852036/100‑1 (MQ=38)
          aGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCgg                                                                                          <  2:1326601/100‑1 (MQ=38)
           gATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCGgc                                                                                         <  2:848901/100‑1 (MQ=38)
           gATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgc                                                                                         >  2:772663/1‑100 (MQ=255)
           gATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTAGCGAAGGATAATGAAGGATGAATCTGCCGGAAGGCGATTAAACGCCATCCGgc                                                                                         <  2:1485769/100‑1 (MQ=255)
            aTTGAACGTGGCAACCCGCGACTGACATTTACAAGCATCGCTTTTATCGAAGGCTAATGAAGCATGAAACTGCCGAATGGCGTTTAATCGCCTTCCgg                                                                                          <  1:283875/98‑1 (MQ=25)
            aTTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCGgca                                                                                        <  1:2296677/100‑1 (MQ=38)
            aTTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCGgca                                                                                        <  2:440415/100‑1 (MQ=38)
            aTTGAAAGTGGCAACCCTCGACGGAGATTTACAGGGGTCGTTTTTATCGAAGGATAATTAAGGATGAAACTGCCGGATGTCGTTTAATCGCCTTCCggt                                                                                         <  1:2052090/99‑2 (MQ=21)
             ttGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGGAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCGgcag                                                                                       <  2:2480319/100‑1 (MQ=37)
              tGAACGTGGCAACCCGCGACGGACATTGACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCGgcagt                                                                                      <  1:930650/100‑2 (MQ=37)
                aaCGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaa                                                                                    >  1:75996/1‑100 (MQ=255)
                aaCGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaa                                                                                    >  2:2290015/1‑100 (MQ=255)
                aaCGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaa                                                                                    <  1:1418345/100‑1 (MQ=255)
                aaCGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaa                                                                                    <  1:1143878/100‑1 (MQ=255)
                   gTGGCAACCCGCGACGGACATTTACAGTGATGGCTGTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCTATTAAACGCCATCCGGCAGcaaatc                                                                                 <  1:1318396/100‑1 (MQ=255)
                    tGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaaatca                                                                                <  2:172805/100‑1 (MQ=255)
                     ggCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaaatcaa                                                                               <  1:2780470/100‑1 (MQ=255)
                      gCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAGCTGCCGGATGGCGTTTAATCGCCTTCCGgcagtttcatcct                                                                              <  2:1792267/100‑10 (MQ=21)
                          cccGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAaataat                                                                          >  2:1910011/1‑100 (MQ=255)
                             gcgACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTAc                                                                       <  1:2519759/100‑1 (MQ=255)
                                aCGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTc                                                                    <  2:1489860/100‑1 (MQ=255)
                                  ggACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGc                                                                  <  2:1402261/100‑1 (MQ=255)
                                    aCATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCg                                                                <  2:87995/100‑1 (MQ=255)
                                       tttACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGaa                                                             >  2:1146192/1‑100 (MQ=255)
                                        ttACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAt                                                            <  1:453882/100‑1 (MQ=255)
                                          aCAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAttt                                                          >  2:1754508/1‑100 (MQ=255)
                                            aGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtttct                                                        >  2:143667/1‑100 (MQ=255)
                                             gggATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCGGCAAATCAAATAATTACTTCGCCGGAAtttctt                                                       <  1:1737601/100‑1 (MQ=255)
                                             gggATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtttct                                                        >  1:1451159/1‑99 (MQ=255)
                                              ggATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtttcttt                                                      >  2:280710/1‑100 (MQ=255)
                                               gATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtttctttc                                                     <  1:2718031/100‑1 (MQ=255)
                                               gATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtttctttc                                                     <  1:574854/100‑1 (MQ=255)
                                                 tGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcag                                                   <  1:2675896/100‑1 (MQ=255)
                                                 tGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcag                                                   <  1:2639476/100‑1 (MQ=255)
                                                 tGGCTTTGATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcag                                                   <  1:2639534/100‑1 (MQ=255)
                                                  ggCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagc                                                  <  1:1348872/100‑1 (MQ=255)
                                                  ggCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagc                                                  <  1:1037364/100‑1 (MQ=255)
                                                   gCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagca                                                 >  2:758417/1‑100 (MQ=255)
                                                    cTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagcag                                                >  1:707150/1‑100 (MQ=255)
                                                    cTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagcag                                                >  1:1360141/1‑100 (MQ=255)
                                                       ttATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtc                                             >  1:1538756/1‑100 (MQ=255)
                                                       ttATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtc                                             <  2:2115647/100‑1 (MQ=255)
                                                       ttATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtc                                             <  1:531152/100‑1 (MQ=255)
                                                        tATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtca                                            <  2:1483101/100‑1 (MQ=255)
                                                         aTCGAAGGATAATGAAGGATGAGACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcag                                           >  1:1693938/1‑100 (MQ=255)
                                                         aTCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcag                                           >  2:1476563/1‑100 (MQ=255)
                                                          tCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagt                                          >  1:784089/1‑100 (MQ=255)
                                                          tCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagt                                          >  1:2752319/1‑100 (MQ=255)
                                                            gAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTGCTTCGCCGGAATTTCTTTCAGCAGTtcagtca                                        >  2:1688493/1‑100 (MQ=255)
                                                            gAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagtca                                        >  2:292809/1‑100 (MQ=255)
                                                            gAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagtca                                        >  1:2058486/1‑100 (MQ=255)
                                                              aGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTcagc                                      <  2:530763/100‑1 (MQ=255)
                                                                gATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTcagcag                                    <  2:2790608/100‑1 (MQ=255)
                                                                  tAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAgtg                                  >  1:2342463/1‑100 (MQ=255)
                                                                     tGAAGGATGAACCTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTcc                               <  1:1277846/100‑1 (MQ=255)
                                                                     tGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTcc                               <  1:60342/100‑1 (MQ=255)
                                                                      gAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCa                              >  1:2042158/1‑100 (MQ=255)
                                                                         ggATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTa                           >  2:1586227/1‑100 (MQ=255)
                                                                          gATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTaa                          >  2:674039/1‑100 (MQ=255)
                                                                            ttaaaCTGCCGGAATGCGATTAACCGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGGCCAGTAATg                        <  1:591283/98‑1 (MQ=255)
                                                                             gAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGa                       <  1:2552010/100‑1 (MQ=255)
                                                                             gAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGa                       >  1:1898420/1‑100 (MQ=255)
                                                                               aaCTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGAcc                     <  1:1123409/100‑1 (MQ=255)
                                                                               aaCTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGAcc                     <  1:651016/100‑1 (MQ=255)
                                                                                gcTGCCGGATGGCGTTTAATCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGAAGTTCAGTCAGCAGTGTCCAGTAATGACCt                    >  1:2902709/2‑100 (MQ=37)
                                                                                   nccGGATGGCGTTTAATCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACg                 >  1:1615197/2‑100 (MQ=38)
                                                                                   gCCGGATGGCGTTTAATCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACg                 >  2:2687841/1‑100 (MQ=38)
                                                                                   gCCGGATGGCGTTTAAGCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACg                 >  1:1050881/1‑100 (MQ=39)
                                                                                      ggAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCtt              >  2:1672468/1‑100 (MQ=255)
                                                                                         aGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCg           <  2:2472061/100‑1 (MQ=255)
                                                                                         aGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCg           >  2:2909165/1‑100 (MQ=255)
                                                                                          ggCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGa          <  1:405566/100‑1 (MQ=255)
                                                                                          ggCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGa          <  2:1884937/100‑1 (MQ=255)
                                                                                            cGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGAtg        <  1:1519136/100‑1 (MQ=255)
                                                                                            cGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGAtg        <  2:898352/100‑1 (MQ=255)
                                                                                            cGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGAtg        <  1:778586/100‑1 (MQ=255)
                                                                                            cGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGAtg        >  2:2102270/1‑100 (MQ=255)
                                                                                            cGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGAtg        >  2:1997661/1‑100 (MQ=255)
                                                                                            cGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGAtg        >  1:2500454/1‑100 (MQ=255)
                                                                                            cGATTAAACGCCATCCGGAAGCAAATCAAATAATTACTTCTCCGTAATGTCTTTCTGCAGTTCAGTCAGCAGTGTCCTTTAATTACCTACGCTTTCg           <  1:1719997/97‑1 (MQ=255)
                                                                                              aTTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGAtgtg      >  2:2045818/1‑100 (MQ=255)
                                                                                              aTTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGAtgtg      >  2:1755762/1‑100 (MQ=255)
                                                                                              aTTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGAtgtg      >  1:2794107/1‑100 (MQ=255)
                                                                                                 aaaCGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTCATGACCTACGCTTTCGATGTGAAc   >  1:1244484/1‑100 (MQ=255)
                                                                                                 aaaCGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAAc   >  1:1438632/1‑100 (MQ=255)
                                                                                                 aaaCGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAAc   >  1:686445/1‑100 (MQ=255)
                                                                                                  aaCGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACt  >  2:1358763/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
TTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACT  >  NC_000913/254134‑254331

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: