Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 254,223 A→T 19.0% 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,2230AT19.0% 281.8 / 58.9 121intergenic (+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 (47/51);  new base T (11/12);  total (58/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.19e-01

Only 100 of 123 total aligned reads displayed.

CGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCG  >  NC_000913/254125‑254322
                                                                                                  |                                                                                                   
 ggcgAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAgg                                                                                                   >  2:1854297/1‑100 (MQ=255)
  gcgAATGTTCCATCACCAGATTGAACGTGGCAAACCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGc                                                                                                  >  1:1924311/1‑100 (MQ=255)
       tGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTa                                                                                             <  2:365946/100‑1 (MQ=255)
         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)
           ccATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACg                                                                                         <  2:1424926/100‑1 (MQ=255)
           ccATCACCAGATTGAACGTGGCAAACCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACg                                                                                         <  1:1864323/100‑1 (MQ=255)
            cATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAAc                                                                                          >  2:876933/1‑98 (MQ=255)
            cATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGc                                                                                        <  2:2172242/100‑1 (MQ=255)
             aTCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGc                                                                                        <  2:876113/99‑1 (MQ=255)
                aCCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATc                                                                                    >  2:1594634/1‑100 (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)
                   aGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCgg                                                                                 <  2:1326601/100‑1 (MQ=38)
                   aGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCgg                                                                                 <  1:1852036/100‑1 (MQ=38)
                   aGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGATTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCgg                                                                                 <  1:2049581/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)
                      ttGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGGAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCGgcag                                                                              <  2:2480319/100‑1 (MQ=37)
                      ttGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCGgcag                                                                              <  2:2400867/100‑1 (MQ=38)
                       tGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgcagc                                                                             <  2:1910617/100‑1 (MQ=255)
                       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)
                             tGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaaatca                                                                       <  1:1330179/100‑1 (MQ=255)
                              ggCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaaatcaa                                                                      <  1:2780470/100‑1 (MQ=255)
                                cAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaaatcaaat                                                                    >  1:1531716/1‑100 (MQ=255)
                                cAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaaatcaaat                                                                    <  2:1761197/100‑1 (MQ=255)
                                   cccGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAaataat                                                                 >  2:1910011/1‑100 (MQ=255)
                                         aCGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTc                                                           <  2:1489860/100‑1 (MQ=255)
                                             aCATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCg                                                       <  1:2532567/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)
                                                        gATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtttctttc                                            <  1:386870/100‑1 (MQ=255)
                                                        gATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtttctttc                                            <  1:574854/100‑1 (MQ=255)
                                                        gATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtttctttc                                            <  1:2718031/100‑1 (MQ=255)
                                                        gATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtttctttc                                            <  2:415485/100‑1 (MQ=255)
                                                         aTGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTca                                           <  1:2740616/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:1037364/100‑1 (MQ=255)
                                                           ggCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagc                                         <  1:1348872/100‑1 (MQ=255)
                                                            gCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagca                                        >  2:758417/1‑100 (MQ=255)
                                                             cTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagcag                                       >  1:1360327/1‑100 (MQ=255)
                                                             cTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagcag                                       >  1:707150/1‑100 (MQ=255)
                                                                ttATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtc                                    >  1:1538756/1‑100 (MQ=255)
                                                                ttATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtc                                    <  2:2115647/100‑1 (MQ=255)
                                                                 tATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtca                                   >  2:1481259/1‑100 (MQ=255)
                                                                  aTCGAAGGATAATGAAGGATGAGACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcag                                  >  1:1693938/1‑100 (MQ=255)
                                                                  aTCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcag                                  >  2:1476563/1‑100 (MQ=255)
                                                                   tCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagt                                 >  1:784739/1‑100 (MQ=255)
                                                                   tCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagt                                 >  1:784089/1‑100 (MQ=255)
                                                                   tCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagt                                 >  1:2752319/1‑100 (MQ=255)
                                                                     gAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagtca                               >  1:2058486/1‑100 (MQ=255)
                                                                     gAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagtca                               >  2:292809/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                      >  2:613600/1‑100 (MQ=255)
                                                                              tGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTcc                      <  1:60342/100‑1 (MQ=255)
                                                                              tGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTcc                      <  1:1443535/100‑1 (MQ=255)
                                                                               gAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCa                     >  1:2042158/1‑100 (MQ=255)
                                                                                 aGGATGAAACTGCCGCAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGt                   <  2:2868836/100‑1 (MQ=255)
                                                                                  ggATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTa                  >  2:1586227/1‑100 (MQ=255)
                                                                                  ggATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTa                  >  2:1728414/1‑100 (MQ=255)
                                                                                    gatttgcTGCCGGATGGCGTTTAATCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAAt                >  1:77820/7‑100 (MQ=25)
                                                                                     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)
                                                                                          cTGCCGGATGGCGTTTAATCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTa          >  2:2578554/1‑100 (MQ=38)
                                                                                          cTGCCGGATGGCGTTTAATCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTa          >  2:2116354/1‑100 (MQ=38)
                                                                                            nccGGATGGCGTTTAATCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACg        >  1:1615197/2‑100 (MQ=38)
                                                                                            gCCGGATGGCGTTTAATCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACg        >  2:1349046/1‑100 (MQ=38)
                                                                                            gCCGGATGGCGTTTAATCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACg        >  2:2687841/1‑100 (MQ=38)
                                                                                            gCCGGATGGCGTTTAAGCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACg        >  1:1050881/1‑100 (MQ=39)
                                                                                              cGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGc       >  1:1838984/1‑99 (MQ=255)
                                                                                               ggAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCtt     >  2:1672468/1‑100 (MQ=255)
                                                                                                  aGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCg  >  2:2909165/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
CGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCG  >  NC_000913/254125‑254322

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: