Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,479,020 A→T 16.3% intergenic (+131/+14) 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,0200AT16.3% 407.0 / 64.5 165intergenic (+131/+14)rraB/yjgMprotein inhibitor of RNase E/GNAT family putative N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (72/65);  new base T (12/15);  total (85/80)
Fisher's exact test for biased strand distribution p-value = 5.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.43e-01

Only 100 of 173 total aligned reads displayed.

CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCT  >  NC_000913/4478921‑4479119
                                                                                                   |                                                                                                   
cTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAaca                                                                                                     <  1:251344/100‑1 (MQ=255)
  gCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACAcg                                                                                                   <  2:4142756/100‑1 (MQ=255)
  gCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGAGAGTTGAacagg                                                                                                   >  2:3139067/1‑98 (MQ=255)
   cAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACAcgc                                                                                                  >  2:3936845/1‑100 (MQ=255)
      ccTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCggg                                                                                               >  2:181223/1‑100 (MQ=255)
      ccTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCggg                                                                                               <  1:3876190/100‑1 (MQ=255)
      ccTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCggg                                                                                               >  1:3904549/1‑100 (MQ=255)
      ccTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCggg                                                                                               <  2:2090807/100‑1 (MQ=255)
      ccTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCggg                                                                                               <  1:1648431/100‑1 (MQ=255)
       cTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCggg                                                                                               <  1:3393879/99‑1 (MQ=39)
       cTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGc                                                                                              >  1:3071565/1‑100 (MQ=39)
       cTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGc                                                                                              <  1:2739711/100‑1 (MQ=39)
       cTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGc                                                                                              >  2:1940402/1‑100 (MQ=255)
        tGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGCa                                                                                             <  1:3309403/100‑1 (MQ=39)
         ggTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAt                                                                                            >  1:496834/1‑100 (MQ=255)
          gTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGCAtt                                                                                           <  2:2209528/100‑1 (MQ=39)
          gTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGCAtt                                                                                           <  1:1378408/100‑1 (MQ=38)
           tGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGCAttc                                                                                          <  2:3614044/100‑2 (MQ=38)
           tGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGCAttc                                                                                          <  1:2022446/100‑2 (MQ=38)
           tGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAttt                                                                                          >  2:212345/1‑100 (MQ=255)
            gAATTTGTAGCACGGATAAGGCATTCATGTCACATCCGGCATAGACAGAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGCAtt                                                                                           <  2:275535/98‑1 (MQ=37)
             aaTTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGCAttcgn                                                                                        <  1:1239266/100‑4 (MQ=38)
              aTTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGCAttcgaa                                                                                       <  1:1283906/100‑5 (MQ=37)
                ttGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTa                                                                                     <  1:3193161/100‑1 (MQ=255)
                ttGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTa                                                                                     >  2:1905921/1‑100 (MQ=255)
                  gTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGGTCAAATGCCCGCGTGTTCAACTCGCGAGCAtacaaatgat                                                                                   >  1:951202/1‑91 (MQ=25)
                   tAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTagag                                                                                  <  1:3604422/100‑1 (MQ=255)
                     gCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCACGAGTTGAACACGCGGGCATTTTTTTAGAGtt                                                                                >  2:61496/1‑100 (MQ=255)
                     gCCAGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGtt                                                                                >  1:1956065/1‑100 (MQ=255)
                       ccGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTcg                                                                              >  2:2961577/1‑100 (MQ=255)
                           aTAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCa                                                                          <  1:838620/100‑1 (MQ=255)
                              aGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACa                                                                       <  2:1352654/100‑1 (MQ=255)
                              aGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACa                                                                       <  1:4030451/100‑1 (MQ=255)
                              aGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACa                                                                       >  1:2716256/1‑100 (MQ=255)
                               ggCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACAt                                                                      >  1:3047672/1‑100 (MQ=255)
                                 cATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATcc                                                                    >  2:1096831/1‑100 (MQ=255)
                                  aTTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCg                                                                   <  1:414010/100‑1 (MQ=255)
                                    tCATGTCACAGCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCa                                                                 <  2:1818591/100‑1 (MQ=255)
                                     cATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCAc                                                                <  2:3072096/100‑1 (MQ=255)
                                        gTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTc                                                             >  1:353062/1‑100 (MQ=255)
                                        gTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTc                                                             <  1:1438227/100‑1 (MQ=255)
                                        gTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTc                                                             >  1:664682/1‑100 (MQ=255)
                                         tCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTca                                                            >  2:3069388/1‑100 (MQ=255)
                                         tCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTca                                                            >  2:1677224/1‑100 (MQ=255)
                                         tCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCAGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTca                                                            >  2:743163/1‑100 (MQ=255)
                                         gatcttCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTca                                                            <  2:1699531/95‑1 (MQ=255)
                                          cacaTCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTcac                                                           >  2:1726387/1‑100 (MQ=255)
                                          cacaTCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTcac                                                           <  2:2502494/100‑1 (MQ=255)
                                            caTCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACaa                                                         >  1:3753006/1‑100 (MQ=255)
                                            caTCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACaa                                                         >  1:4026844/1‑100 (MQ=255)
                                             aTCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAAt                                                        <  1:3598575/100‑1 (MQ=255)
                                              tCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTcaagatc                                                       >  1:1699535/1‑95 (MQ=255)
                                              tCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCAGAACTTCACAATc                                                       >  1:1130577/1‑100 (MQ=255)
                                               ccGGCATAAACAAGACGTACGTTGTTAATGTTCGGATGCCCGCGAGTTGAACACGCGGGAGTTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCg                                                      <  2:700809/100‑1 (MQ=255)
                                                 ggCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGAc                                                    <  2:1130577/100‑1 (MQ=255)
                                                  gCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACa                                                   <  2:526815/100‑1 (MQ=255)
                                                  gCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACa                                                   >  1:1920355/1‑100 (MQ=255)
                                                    aTAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATg                                                 <  2:3511175/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)
                                                         cAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCg                                            >  1:1179052/1‑100 (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)
                                                            aaCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGc                                         >  1:2462192/1‑100 (MQ=255)
                                                             acgtacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCa                                        <  2:3898468/100‑1 (MQ=255)
                                                              cgtacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAg                                       <  2:505190/100‑1 (MQ=255)
                                                              cgtacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAg                                       <  1:4073314/100‑1 (MQ=255)
                                                              cgtacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAg                                       >  1:778289/1‑100 (MQ=255)
                                                               gtacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGc                                      <  1:2887582/100‑1 (MQ=255)
                                                                  cgtTGTTAATGTTCGAATGCACGCGAGTTGAACACGCGGGCATTTTTTTAGAGGTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCCGCCaa                                   >  2:3607754/1‑100 (MQ=255)
                                                                     tgttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAgcg                                <  2:3028615/100‑1 (MQ=255)
                                                                      gttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAgcgc                               >  1:3356229/1‑100 (MQ=255)
                                                                        tAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCAt                             <  2:3356213/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)
                                                                               cGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCg                        >  2:2218849/1‑98 (MQ=255)
                                                                               aaaaTGCCCGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGa                       >  2:1294263/3‑99 (MQ=38)
                                                                                gAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGata                     >  1:886121/1‑100 (MQ=255)
                                                                                gAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGata                     <  1:2357265/100‑1 (MQ=255)
                                                                                   tGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGc                  >  2:2850055/1‑100 (MQ=255)
                                                                                   tGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGc                  <  2:2679497/100‑1 (MQ=255)
                                                                                    gCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCt                 <  1:3690289/100‑1 (MQ=255)
                                                                                     cccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTc                >  2:1043499/1‑100 (MQ=255)
                                                                                     cccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTc                >  1:831366/1‑100 (MQ=255)
                                                                                       cgcgTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCaa              >  1:239966/1‑100 (MQ=255)
                                                                                       cgcgAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACCTTGCCCGTGCAGCCAAGCGCATAGTCGATATg                   >  2:2002819/1‑95 (MQ=255)
                                                                                       cgcgAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGatatattca               >  2:857932/1‑94 (MQ=255)
                                                                                       cgcgAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCaa              >  1:727151/1‑100 (MQ=255)
                                                                                         cgAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAg            <  1:3836962/100‑1 (MQ=255)
                                                                                           aGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGcc          >  1:1022637/1‑100 (MQ=255)
                                                                                              tGAACACGCGGGCATTTTTTTTTTTTTTTTTCAACATCCGCACTGCACAATCTACATGGCCCTTACAGCCAAGCGCACCGCCGCTCTTCTCAAATcccca       >  2:3146476/1‑98 (MQ=255)
                                                                                              tGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCaa       >  2:2046183/1‑100 (MQ=255)
                                                                                                aaCACGCGGGCATTTTTTTAGAGTTCGCGCAACAACCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAAt     <  2:742741/100‑1 (MQ=255)
                                                                                                 acacGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATg    >  1:2215185/1‑100 (MQ=255)
                                                                                                  cacGCGGGCATTTTTTTAGAGTTAGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGc   >  2:1149612/1‑100 (MQ=255)
                                                                                                   acGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCt  <  1:1905924/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCT  >  NC_000913/4478921‑4479119

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: