Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,479,012 A→T 15.7% intergenic (+123/+22) rraB → / ← yjgM protein inhibitor of RNase E/GNAT family putative N‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,479,0120AT15.7% 454.9 / 67.4 178intergenic (+123/+22)rraB/yjgMprotein inhibitor of RNase E/GNAT family putative N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (73/77);  new base T (12/16);  total (85/93)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.37e-01

Only 100 of 179 total aligned reads displayed.

GTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCC  >  NC_000913/4478914‑4479111
                                                                                                  |                                                                                                   
gTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAg                                                                                                    <  2:2900319/100‑1 (MQ=255)
 tAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAATGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGt                                                                                                   <  2:2375747/100‑1 (MQ=255)
 tAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGt                                                                                                   <  1:2242956/100‑1 (MQ=255)
   gTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGtag                                                                                                 <  1:2834769/100‑3 (MQ=255)
   gTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTg                                                                                                 <  1:3069929/100‑1 (MQ=255)
   gTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGCTGTTAATGTTCGAATGCCCGCGAGTTg                                                                                                 >  1:285708/1‑100 (MQ=255)
      ccTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAac                                                                                              >  1:4161166/1‑100 (MQ=255)
      ccTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAac                                                                                              <  1:3409099/100‑1 (MQ=255)
       cTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAaca                                                                                             <  1:251344/100‑1 (MQ=255)
         gCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACAcg                                                                                           <  2:4142756/100‑1 (MQ=255)
         gCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGAGAGTTGAacagg                                                                                           >  2:3139067/1‑98 (MQ=255)
             ccTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCggg                                                                                       <  1:1648431/100‑1 (MQ=255)
              cTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGc                                                                                      <  1:2739711/100‑1 (MQ=39)
              cTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGc                                                                                      >  2:1940402/1‑100 (MQ=255)
              cTGGTGAATTTGTAGCCCGGATAAGGCATTCATGGCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGc                                                                                      <  1:3749270/100‑1 (MQ=38)
               tGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGCa                                                                                     <  1:3309403/100‑1 (MQ=39)
                ggTGAAGTTGTAGCCCGGAGAAGGCAGTCATGTCGCATCCGGCATAAACAAAACGGACGTTGTTAATGTTCGAATGCCCGCGTGTTCAGCTCGCGGGc                                                                                      <  2:27287/98‑1 (MQ=25)
                 gTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGCAtt                                                                                   <  1:3523872/100‑1 (MQ=38)
                  tGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGCAttc                                                                                  <  1:2022446/100‑2 (MQ=38)
                   gAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAtttt                                                                                 >  1:757657/1‑100 (MQ=255)
                    aaTTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGCAttcgn                                                                                <  1:1239266/100‑4 (MQ=38)
                       ttGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTa                                                                             >  2:1905921/1‑100 (MQ=255)
                       ttGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTa                                                                             <  1:3193161/100‑1 (MQ=255)
                         gTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGCAttcgaacatt                                                                           <  2:3071552/100‑9 (MQ=25)
                         gTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGGTCAAATGCCCGCGTGTTCAACTCGCGAGCAtacaaatgat                                                                           >  1:951202/1‑91 (MQ=25)
                          tagctcttaTAATGCATTAATGTCACATCCGACATAAACAAAACGTACGTTGTTAATTTTCGAATCCCCCCTAGTTGCACACTCGTGCATTTTTTTaaag                                                                          <  2:219450/92‑1 (MQ=255)
                            gCCAGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGtt                                                                        >  1:1956065/1‑100 (MQ=255)
                              ccGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTcg                                                                      >  2:2961577/1‑100 (MQ=255)
                               cGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTcgc                                                                     <  1:40265/100‑1 (MQ=255)
                                 gATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCgcgc                                                                   <  1:1267239/100‑1 (MQ=255)
                                  aTAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCa                                                                  <  1:838620/100‑1 (MQ=255)
                                     aGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACa                                                               <  1:4030451/100‑1 (MQ=255)
                                        cATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATcc                                                            >  2:1096831/1‑100 (MQ=255)
                                         aTTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCg                                                           <  2:2638663/100‑1 (MQ=255)
                                         aTTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCg                                                           <  1:1300560/100‑1 (MQ=255)
                                           tCATGTCACAGCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCa                                                         <  2:1818591/100‑1 (MQ=255)
                                            cATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCAc                                                        <  2:3072096/100‑1 (MQ=255)
                                               gTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTc                                                     >  1:664682/1‑100 (MQ=255)
                                                tCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCAGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTca                                                    >  2:743163/1‑100 (MQ=255)
                                                gatcttCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTca                                                    <  2:1699531/95‑1 (MQ=255)
                                                 cacaTCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTcac                                                   <  2:2502494/100‑1 (MQ=255)
                                                 cacaTCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTcac                                                   >  2:1726387/1‑100 (MQ=255)
                                                 cacaTCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTcac                                                   <  2:3386101/100‑1 (MQ=255)
                                                   caTCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACaa                                                 >  1:1095072/1‑100 (MQ=255)
                                                    aTCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAAt                                                <  1:3598575/100‑1 (MQ=255)
                                                     tCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATc                                               <  1:4032040/100‑1 (MQ=255)
                                                      ccGGCATAAACAAGACGTACGTTGTTAATGTTCGGATGCCCGCGAGTTGAACACGCGGGAGTTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCg                                              <  2:700809/100‑1 (MQ=255)
                                                      ccGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCg                                              >  2:591119/1‑100 (MQ=255)
                                                      ccGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCg                                              >  1:1676892/1‑100 (MQ=255)
                                                         gCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACa                                           >  1:1920355/1‑100 (MQ=255)
                                                          cATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACAt                                          >  2:3675394/1‑100 (MQ=255)
                                                           aTAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATg                                         <  1:697318/100‑1 (MQ=255)
                                                              aaCAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGCCATTTTTTTAGAGTTCACGCAACATCCGCACTTCACAATCGACATGGcc                                      >  1:2346531/1‑100 (MQ=255)
                                                              aaCAAAACGTACGTTGTTAACGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGcc                                      >  1:3489640/1‑100 (MQ=255)
                                                               aCAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGccc                                     >  2:2898624/1‑100 (MQ=255)
                                                               aCAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGccc                                     >  1:3190845/1‑100 (MQ=255)
                                                                cAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGccc                                     <  2:1309449/99‑1 (MQ=255)
                                                                cAAAAAGTAGGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACACTGGCATCTTTTTGGTGTTAGCGCAATATCTGCACTTCACAATCGACATGGCCCg                                    <  2:776931/100‑1 (MQ=255)
                                                                  aaaCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGcc                                      <  2:1031647/96‑1 (MQ=255)
                                                                  aaaCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTg                                  <  1:1023493/100‑1 (MQ=255)
                                                                  aaaCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTg                                  >  1:2009548/1‑100 (MQ=255)
                                                                  aaaCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTg                                  >  1:2420859/1‑100 (MQ=255)
                                                                   aaCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGc                                 >  1:2462192/1‑100 (MQ=255)
                                                                    acgtacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCa                                <  2:3898468/100‑1 (MQ=255)
                                                                    acgtacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCa                                <  1:1396598/100‑1 (MQ=255)
                                                                     cgtacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAg                               <  2:505190/100‑1 (MQ=255)
                                                                         cgtTGTTAATGTTCGAATGCACGCGAGTTGAACACGCGGGCATTTTTTTAGAGGTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCCGCCaa                           >  2:3607754/1‑100 (MQ=255)
                                                                          gttgttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAg                          <  2:303579/100‑1 (MQ=255)
                                                                            tgttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAgcg                        >  1:3882737/1‑100 (MQ=255)
                                                                            tgttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAgcg                        <  1:3832749/100‑1 (MQ=255)
                                                                            tgttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAgcg                        <  1:2603946/100‑1 (MQ=255)
                                                                             gttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAgcgc                       >  1:3356229/1‑100 (MQ=255)
                                                                             gttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAgcgc                       <  1:1177071/100‑1 (MQ=255)
                                                                               tAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCAt                     <  2:3356213/100‑1 (MQ=255)
                                                                               tAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCAt                     <  1:1726391/100‑1 (MQ=255)
                                                                                tctaaaaaaaTGCCCGCGTGTTCAACTCGCGGGCATTTTTATAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATa                    >  2:4035088/9‑100 (MQ=25)
                                                                                 ctaaaaaaaTGCCCGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAg                   >  1:1465455/8‑100 (MQ=37)
                                                                                  tGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGt                  >  2:3864921/1‑100 (MQ=255)
                                                                                   gTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTc                 <  2:3582735/100‑1 (MQ=255)
                                                                                    aaaaaaTGCCCGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCg                >  1:828562/5‑100 (MQ=37)
                                                                                      cGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCg                >  2:2218849/1‑98 (MQ=255)
                                                                                      cGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGat              <  1:3675413/100‑1 (MQ=255)
                                                                                      aaaaTGCCCGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGa               >  2:1294263/3‑99 (MQ=38)
                                                                                       gAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGata             <  1:2357265/100‑1 (MQ=255)
                                                                                       gAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGata             >  1:886121/1‑100 (MQ=255)
                                                                                       gAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGAGCAACATCCGCACGTTACCATAGACACGGCCCGTGCAGCCAAGCGCATAGTCGa               >  2:3970605/1‑98 (MQ=255)
                                                                                         aTGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATg           <  2:662164/100‑1 (MQ=255)
                                                                                          tGCCCGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGc          >  2:1773583/1‑100 (MQ=255)
                                                                                          tGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGc          >  2:2850055/1‑100 (MQ=255)
                                                                                          tGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGc          <  1:2984776/100‑1 (MQ=255)
                                                                                          tGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTACGCGCAACATTCGCACTTCACAATCGACATGGGCCGGGCAGCGAAGCGCATAGTCGATATGc          <  2:2081396/100‑1 (MQ=255)
                                                                                           gCCCGCGAGTTGAACACGCGGTCACTTTTTTAGAGTTCGCGCACCATCCGGACTTCACAATCGACAGGGCCCTTGCAGCCAAGCGCATAGTCGatat            <  2:1540165/97‑1 (MQ=255)
                                                                                            cccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTc        >  2:1043499/1‑100 (MQ=255)
                                                                                            cccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTc        >  1:411544/1‑100 (MQ=255)
                                                                                              cgcgAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACCTTGCCCGTGCAGCCAAGCGCATAGTCGATATg           >  2:2002819/1‑95 (MQ=255)
                                                                                              cgcgAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGatatattca       >  2:857932/1‑94 (MQ=255)
                                                                                              cgcgAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCaa      >  2:209655/1‑100 (MQ=255)
                                                                                                cgAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAATATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAg    <  1:454571/100‑1 (MQ=255)
                                                                                                cgAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAg    <  1:3836962/100‑1 (MQ=255)
                                                                                                  aGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGcc  >  1:1022637/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
GTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCC  >  NC_000913/4478914‑4479111

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: