Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,479,020 A→T 15.2% 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,0200AT15.2% 400.4 / 58.8 158intergenic (+131/+14)rraB/yjgMprotein inhibitor of RNase E/GNAT family putative N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (69/65);  new base T (11/13);  total (80/78)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-02

Only 100 of 168 total aligned reads displayed.

CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCT  >  NC_000913/4478921‑4479119
                                                                                                   |                                                                                                   
cTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAACGTTCGAATGCCCGCGAGTTGAaca                                                                                                     >  2:2042629/1‑100 (MQ=255)
cTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGATGTTAATGTTCGAATGCCCGCGAGTTGAaca                                                                                                     >  1:1062927/1‑100 (MQ=255)
  gCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACAcg                                                                                                   >  2:1114778/1‑100 (MQ=255)
   cAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACAcgc                                                                                                  <  1:3289151/100‑1 (MQ=255)
      ccTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCggg                                                                                               <  1:2104332/100‑1 (MQ=38)
       cTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGc                                                                                              <  2:1463825/100‑1 (MQ=39)
       cTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGc                                                                                              <  1:145154/100‑1 (MQ=39)
        tGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCa                                                                                             >  1:3420631/1‑100 (MQ=255)
         ggTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAt                                                                                            >  1:2121954/1‑100 (MQ=255)
          gTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAtt                                                                                           >  1:1664897/1‑100 (MQ=255)
          gTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAtt                                                                                           >  2:2761232/1‑100 (MQ=255)
           tGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGCAttc                                                                                          <  2:3093795/100‑2 (MQ=38)
              aTTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAtttttt                                                                                       <  2:1445764/100‑1 (MQ=255)
              aTTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAtttttt                                                                                       >  1:3237353/1‑100 (MQ=255)
              aTTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAtttttt                                                                                       >  2:1813454/1‑100 (MQ=255)
               tttGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGCAttcgaac                                                                                      <  2:1147900/100‑6 (MQ=37)
               tttGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAttttttt                                                                                      <  1:2761256/100‑1 (MQ=255)
               tttGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAttttttt                                                                                      <  2:853560/100‑1 (MQ=255)
                ttGTAGCCCGGATACGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTa                                                                                     <  2:568651/100‑1 (MQ=255)
                 tGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTag                                                                                    <  1:2711285/100‑1 (MQ=255)
                 cgTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTGCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTag                                                                                    <  1:2920890/99‑1 (MQ=255)
                   tAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTagag                                                                                  >  2:2245228/1‑100 (MQ=255)
                    aGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGt                                                                                 >  1:2343940/1‑100 (MQ=255)
                    aGCCCGGATAAGGCATTCATGACACATCCAGCATAAACAAAACGTACGTTGTTAATGTTAGAATGCACGCGAGTTGAAAACGCGGGCATTTTTTTAGAGt                                                                                 >  2:1669504/1‑100 (MQ=255)
                     gCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAttttttttttttt                                                                                >  2:904936/1‑94 (MQ=255)
                     gCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGtt                                                                                >  1:1078235/1‑100 (MQ=255)
                     gCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGtt                                                                                >  1:2987307/1‑100 (MQ=255)
                     gCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCACGCGAGTTGAACACACGGGAATTTTTTTagagat                                                                                >  1:674254/1‑98 (MQ=39)
                         ggATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTcgcg                                                                            >  1:1413419/1‑100 (MQ=255)
                         ggATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTAGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTcgcg                                                                            <  2:3527872/100‑1 (MQ=255)
                           aTAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCgcgc                                                                           <  2:780596/99‑1 (MQ=255)
                             aaGGCATTCATGTCACATCCGGAATAAACAAAACGTACGATGTTAATATTCGAATGCCAGCGAGATGAAAACGCGGGCATTTTTTTAGAGTTCGCGCaa                                                                         >  2:3429363/1‑99 (MQ=255)
                              aGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACa                                                                       <  2:2121935/100‑1 (MQ=255)
                              aGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACa                                                                       <  1:3306038/100‑1 (MQ=255)
                              aGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACa                                                                       >  2:3540341/1‑100 (MQ=255)
                               ggCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACAt                                                                      >  1:1030240/1‑100 (MQ=255)
                                 tccgatcttCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATcc                                                                    <  2:2118348/92‑1 (MQ=255)
                                  aTTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCg                                                                   <  1:2362848/100‑1 (MQ=255)
                                   ttCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGc                                                                  <  1:3092205/100‑1 (MQ=255)
                                   ttCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGc                                                                  <  1:2827633/100‑1 (MQ=255)
                                     cATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCAc                                                                <  1:2078378/100‑1 (MQ=255)
                                     cATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAAGGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCAc                                                                >  2:2172066/1‑100 (MQ=255)
                                      aTGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACt                                                               >  1:3225536/1‑100 (MQ=255)
                                      aTGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCAAGTTAAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACt                                                               >  1:172498/1‑100 (MQ=255)
                                        gTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTc                                                             <  1:717744/100‑1 (MQ=255)
                                        gTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTAAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTc                                                             >  1:2335/1‑100 (MQ=255)
                                        gTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTAAACAAGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACtt                                                              >  1:64349/1‑99 (MQ=255)
                                        gTCACATCCGGCATAAAAAAAACGTAAGTGGATAATGTTCGAATGCACGCGAGTTGAAAACGCGGGCATTTTTTTAGAGTTAGCGCAACATCCGCACtt                                                              >  2:2249594/1‑99 (MQ=255)
                                         tCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTca                                                            >  1:101432/1‑100 (MQ=255)
                                          cacaTCCGGCATAAACAAAACGTACGTTGTGAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTcac                                                           <  2:1271041/100‑1 (MQ=255)
                                           acaTCAGGCATAAACAAAACGTACGATGTTAATATTCGAATGCCCGCGAGTTGAACACGCGGTCATTTTTTTAGAGTTCGCGCAACATCCGCACTTca                                                            <  2:2326734/98‑1 (MQ=255)
                                              tCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTAGAGTTCGCGCAACATCCGCAATTCACAATCg                                                      >  1:3638009/1‑100 (MQ=255)
                                               ccGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCg                                                      <  1:1181118/100‑1 (MQ=255)
                                                cGGCATAAACTATACTTACTTTGTTATTGTTCGAATGGCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTGCGCGCAACATCCGCACTTCACAATCGa                                                     <  1:582245/100‑1 (MQ=255)
                                                  gCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACa                                                   >  1:1112043/1‑100 (MQ=255)
                                                   cATAAACAAAACGTACGTTGTTAATGTTCAAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGAc                                                    <  1:3850209/98‑1 (MQ=255)
                                                       aaCAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGcc                                              >  1:1247805/1‑100 (MQ=255)
                                                       aaCAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAAAGATCGCGCAACATCCGCACTTAACAATCGACATGGcc                                              >  1:3525893/1‑100 (MQ=255)
                                                       aaCAAAACGTACGTTGTTAATGTTCGAAGGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGcc                                              <  1:3214201/100‑1 (MQ=255)
                                                            aaCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGc                                         <  1:3161780/100‑1 (MQ=255)
                                                             acgtacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCa                                        <  1:2564614/100‑1 (MQ=255)
                                                               gtacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGc                                      >  1:2614217/1‑100 (MQ=255)
                                                               gtacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGc                                      <  2:2063582/100‑1 (MQ=255)
                                                               gtacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGc                                      <  1:522497/100‑1 (MQ=255)
                                                                 acgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCa                                    <  1:1899958/100‑1 (MQ=255)
                                                                   gttgttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAg                                  >  1:2210267/1‑100 (MQ=255)
                                                                   gttgttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAg                                  <  1:3899543/100‑1 (MQ=255)
                                                                     tgttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAgcg                                >  2:2481207/1‑100 (MQ=255)
                                                                     tgtgAATGTTCGAATGCCCGCGAGTTGTACACGCGGGCATTTTTTTAGAGATCGCGCAACATCAGAACTTCAAAAACGACATGGCCCGTGCAGCCAAg                                  >  1:2152234/1‑98 (MQ=255)
                                                                      gttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAgcgc                               >  2:3136974/1‑100 (MQ=255)
                                                                      gttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAgcgc                               >  2:1003234/1‑100 (MQ=255)
                                                                      gttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAgcgc                               >  1:3768907/1‑100 (MQ=255)
                                                                        tAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCAt                             >  1:2527216/1‑100 (MQ=255)
                                                                        naaTGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACAACCACACTTCACAATCGACATGGCCCGTGCAACCAAGCGCa                              >  1:1329182/2‑99 (MQ=255)
                                                                          aTGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAg                           <  1:3651860/100‑1 (MQ=255)
                                                                           tGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAg                           >  1:3697337/1‑99 (MQ=255)
                                                                            gTTCGAATGCCCGTGAGTTGAACACGCGGGCATTTTTTTGGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTc                         <  1:474010/100‑1 (MQ=255)
                                                                            aaaaaaaTGCCCGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGt                          >  1:1237961/6‑99 (MQ=37)
                                                                            aaaaaaaTGCCCGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTc                         >  1:3482527/6‑100 (MQ=37)
                                                                             ttCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCg                        <  1:1665246/100‑1 (MQ=255)
                                                                               cGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGat                      <  1:708257/100‑1 (MQ=255)
                                                                                 aaTGCCCGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGatat                    >  2:3164202/1‑100 (MQ=38)
                                                                                  aTGCCCGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATg                   >  2:769467/1‑100 (MQ=39)
                                                                                  aTGCCCGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAAACGACATGGCACGTGCAGCCAAACGCAAAGACGaaa                     >  1:3672927/1‑96 (MQ=25)
                                                                                    gCCCGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGCTATGCt                 >  2:3400465/1‑100 (MQ=39)
                                                                                    gCCCGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCt                 >  1:1971522/1‑100 (MQ=39)
                                                                                    gCCCGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTAGATATGCt                 >  2:3235284/1‑100 (MQ=39)
                                                                                     cccGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTc                >  1:3078451/1‑100 (MQ=255)
                                                                                      ccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATATTCGATACGCTc                >  2:3717343/1‑99 (MQ=255)
                                                                                       cgcgAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCaa              >  1:2503076/1‑100 (MQ=255)
                                                                                        gcgAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCaaa             >  1:2943304/1‑100 (MQ=255)
                                                                                         cgtgTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAg            <  2:3482488/97‑1 (MQ=255)
                                                                                          gAGTTGAACACTCGGGCATTTTTTTAGAGTTCTCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCTGTCAAGTTCATAGTCGATATGCTCAAAGc           <  2:1852231/100‑1 (MQ=255)
                                                                                            gTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGccc         >  2:1591320/1‑100 (MQ=255)
                                                                                             ttGAACACGCTGTCGTTTTTTTAGAGTTTGCGCAACATCCGAACTTCACAATCGACATGGCACGTGCAGCCAATCGCATAGTCGATATGCTCAAAGccc         <  2:3900404/99‑1 (MQ=255)
                                                                                                  cacGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGc   <  1:1591335/100‑1 (MQ=255)
                                                                                                  cacGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGc   <  1:1813469/100‑1 (MQ=255)
                                                                                                  cacGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGc   >  2:1494160/1‑100 (MQ=255)
                                                                                                   acGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCt  <  2:3768868/100‑1 (MQ=255)
                                                                                                   acGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCt  <  1:2087120/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 < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: