Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,851,006 A→G intergenic (‑391/‑537) ivbL ← / → emrD ilvB operon leader peptide/2‑module integral membrane pump

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,851,0060AG95.5% 1744.3 / 45.0 594intergenic (‑391/‑537)ivbL/emrDilvB operon leader peptide/2‑module integral membrane pump
Reads supporting (aligned to +/- strand):  ref base A (14/13);  new base G (251/316);  total (265/329)
Fisher's exact test for biased strand distribution p-value = 5.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 628 total aligned reads displayed.

GTCCAAATATTACTGTTTATTTATACAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCACTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGA  >  NC_000913/3850916‑3851098
                                                                                          |                                                                                            
 tCCAAATATTACTGTTTATTTATACAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTa                                                                                           >  1:4395611/1‑93 (MQ=255)
 tCCAAATATTACTGTTTATTTATACAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTa                                                                                           <  2:1700149/93‑1 (MQ=255)
  ccAAATATTACTGTTTATTAATACAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTaa                                                                                          <  2:241630/93‑1 (MQ=255)
   cAAATATTACTGTTTATTTATACAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAg                                                                                         <  1:1871692/93‑1 (MQ=255)
   cAAATATTACTGTTTATTTATACAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAg                                                                                         >  1:5980780/1‑93 (MQ=255)
     aaTATTACTGTTTATTTATACAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTc                                                                                       >  2:6113547/1‑93 (MQ=255)
     aaTATTACTGTTTATTTATACAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTc                                                                                       <  1:4271895/93‑1 (MQ=255)
     aaTATTACTGTTTATTTATACAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTc                                                                                       <  2:5512171/93‑1 (MQ=255)
       tatTACTGTTTATTTATACAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCt                                                                                     <  1:804767/93‑1 (MQ=255)
       tatTACTGTTTATTTATACAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCt                                                                                     >  1:7346497/1‑93 (MQ=255)
         ttACTGTTTATTTATACAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTgg                                                                                   <  1:75108/93‑1 (MQ=255)
          tACTGTTTATTTATACAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGc                                                                                  <  2:1999958/93‑1 (MQ=255)
           aCTGTTTATTTATACAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCt                                                                                 >  1:4667973/1‑93 (MQ=255)
           aCTGTTTATTTATACAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCACTAAGTCCTGGCt                                                                                 <  1:9257395/93‑1 (MQ=255)
            cTGTTTATTTATACAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTAGCCGGCGCGCTAAGTCCTGGCTg                                                                                >  1:5072627/1‑93 (MQ=255)
             tGTTTATTTATACAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGa                                                                               >  1:8530167/1‑93 (MQ=255)
              gTTTATTTATACAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGaa                                                                              >  1:1026424/1‑93 (MQ=255)
              gTTTATTTATACAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGaa                                                                              >  1:7283189/1‑93 (MQ=255)
              gTTTATTTATACAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGaa                                                                              <  2:777259/93‑1 (MQ=255)
               tttatttaTACAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGaaa                                                                             <  2:9118383/93‑1 (MQ=255)
                ttatttatACAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAAc                                                                            <  2:3468720/93‑1 (MQ=255)
                ttatttatACAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAAc                                                                            <  2:6501736/93‑1 (MQ=255)
                 tatttatACAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACg                                                                           <  1:4012341/93‑1 (MQ=255)
                   tttatACAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACggg                                                                         <  2:3906872/93‑1 (MQ=255)
                   tttatACAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACggg                                                                         >  2:4094906/1‑93 (MQ=255)
                   tttatACAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACggg                                                                         <  1:8295051/93‑1 (MQ=255)
                     tataCAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGggtg                                                                       <  2:1608784/93‑1 (MQ=255)
                     tataCAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGggtg                                                                       <  2:10045650/93‑1 (MQ=255)
                      ataCAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGggtgg                                                                      <  2:4311984/93‑1 (MQ=255)
                      ataCAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGggtgg                                                                      <  1:2498352/93‑1 (MQ=255)
                       taCAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGggtggt                                                                     <  1:8812315/93‑1 (MQ=255)
                       taCAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGggtggt                                                                     <  2:3100553/93‑1 (MQ=255)
                       taCAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGggtggt                                                                     >  1:1336449/1‑93 (MQ=255)
                       taCAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGggtggt                                                                     <  1:2880878/93‑1 (MQ=255)
                       taCAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGggtggg                                                                     >  1:8462122/1‑92 (MQ=255)
                       taCAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTAGCTGAAACGggcggc                                                                     >  1:6654607/1‑92 (MQ=255)
                        aCAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGgtggtg                                                                    >  1:6185232/1‑93 (MQ=255)
                        aCAGTAAACTTCTATAATATCACCGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGATGAAACGggtggtc                                                                    >  2:6719651/1‑92 (MQ=255)
                          aGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGcc                                                                  >  2:9309285/1‑93 (MQ=255)
                          aGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGCGcc                                                                  >  2:2656574/1‑93 (MQ=255)
                             aaaCTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCACTAAGTCCTGGCTGAAACGGGTGGTGCCGTc                                                               >  1:571945/1‑93 (MQ=255)
                              aaCTTCTATAATATCACTGTACGCATTGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCa                                                              >  2:1550912/1‑93 (MQ=255)
                              aaCTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCa                                                              >  1:10001309/1‑93 (MQ=255)
                               aCTTCTATAATATCCCTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAg                                                             <  1:1731055/93‑1 (MQ=255)
                               aCTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTTGCTGAAACGGGTGGTGCCGTCAg                                                             >  1:4810536/1‑93 (MQ=255)
                               aCTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAg                                                             <  2:4824230/93‑1 (MQ=255)
                                cTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAgc                                                            <  2:7593407/93‑1 (MQ=255)
                                cTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAgc                                                            <  2:5771992/93‑1 (MQ=255)
                                  tCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAgcgc                                                          >  1:3072890/1‑93 (MQ=255)
                                   cTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGcc                                                         >  2:4748251/1‑93 (MQ=255)
                                       aataTCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTaa                                                     >  2:8434539/1‑93 (MQ=255)
                                        atatCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAAc                                                    >  2:1311577/1‑93 (MQ=255)
                                        atatCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCACTAAGTCCTGGCTGAAATGGGTGGTGCCGTCAGCGCCTTAAc                                                    >  2:5426501/1‑93 (MQ=255)
                                          atCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAAccc                                                  <  1:9383506/93‑1 (MQ=255)
                                             aCTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCcgc                                               >  2:2721100/1‑93 (MQ=255)
                                               tGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCcgcgg                                             >  1:3761587/1‑92 (MQ=255)
                                                 tACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGa                                           <  2:6314296/93‑1 (MQ=255)
                                                 tACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCACTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGa                                           >  2:2896421/1‑93 (MQ=255)
                                                   cGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGGCCTGGCTGAAGCGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGc                                         <  1:9114127/93‑1 (MQ=255)
                                                    gCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGca                                        <  1:6550203/93‑1 (MQ=255)
                                                    gCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTAGTGCCGTCAGCGCCTTAACCCCGCGTGAGca                                        <  2:5980779/93‑1 (MQ=255)
                                                      aaTGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGcaca                                      >  1:7007346/1‑93 (MQ=255)
                                                      aaTGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCACTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGcaca                                      >  1:6884202/1‑93 (MQ=255)
                                                       aTGTGTTATGCTGGGGCCGCACCGTTCCCCGGCGCGCTAAGTCCTGCCTGAAACGCGTGGTGCCGTCAGCGCCTTACCCCCGCGTGAGCacac                                     <  2:8902836/93‑1 (MQ=255)
                                                       aTGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCacac                                     <  2:7948282/93‑1 (MQ=255)
                                                       aTGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCacac                                     >  1:9889406/1‑93 (MQ=255)
                                                        tgtgTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACt                                    >  2:366576/1‑93 (MQ=255)
                                                        tgtgTTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACt                                    >  2:6946779/1‑93 (MQ=255)
                                                         tggtTATGCGGGTGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACtg                                   <  1:4919911/91‑1 (MQ=255)
                                                         gtgtTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACtg                                   >  1:3603675/1‑93 (MQ=255)
                                                         gtgtTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACtg                                   <  2:6758871/93‑1 (MQ=255)
                                                           gtTATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACtgtg                                 >  2:3249360/1‑93 (MQ=255)
                                                            ttATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTgtgt                                >  1:9686958/1‑93 (MQ=255)
                                                            ttATGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTgtgt                                <  2:4948854/93‑1 (MQ=255)
                                                              aTGCGGGGGCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTa                              <  1:7724468/93‑1 (MQ=255)
                                                               tGCGGTGGCCGCATCGTTGCCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCCCACTGTGTTAt                             <  2:2160201/93‑1 (MQ=255)
                                                                   ggggCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTca                         >  1:1442064/1‑93 (MQ=255)
                                                                    gggCCGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTcaa                        <  1:8968352/93‑1 (MQ=255)
                                                                    gggCCGCATCGTTACCCGGCGCACTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTcaa                        <  2:7360273/93‑1 (MQ=255)
                                                                        cGCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAg                    >  1:2319358/1‑93 (MQ=255)
                                                                         gCATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGc                   <  2:1107808/93‑1 (MQ=255)
                                                                          cATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGca                  <  1:7882492/93‑1 (MQ=255)
                                                                          cATCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGca                  >  2:3329190/1‑93 (MQ=255)
                                                                           aTCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGcac                 >  2:7094475/1‑93 (MQ=255)
                                                                           aTCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGcac                 <  2:9588944/93‑1 (MQ=255)
                                                                           aTCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGcac                 <  1:8766009/93‑1 (MQ=255)
                                                                            tCGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGcaca                >  2:6412397/1‑93 (MQ=255)
                                                                             cGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACaa               <  2:5288117/93‑1 (MQ=255)
                                                                             cGTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGGGTTATGTCAACAAGCACaa               <  1:5366985/93‑1 (MQ=255)
                                                                              gTTACCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAAc              >  1:5853112/1‑93 (MQ=255)
                                                                                 aCCCGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGtt           <  1:3390702/93‑1 (MQ=255)
                                                                                  cccGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGttt          <  1:5531022/93‑1 (MQ=255)
                                                                                  cccGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGttt          >  2:8659966/1‑93 (MQ=255)
                                                                                   ccGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTtc         >  2:8214062/1‑93 (MQ=255)
                                                                                   ccGGCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTtc         <  2:1865110/93‑1 (MQ=255)
                                                                                     ggCGCGCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTtctc       >  1:1982169/1‑93 (MQ=255)
                                                                                      gcgcgcTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTcc      <  1:3254373/93‑1 (MQ=255)
                                                                                       cgcgCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCt     <  1:1146702/93‑1 (MQ=255)
                                                                                          aCTANGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTga  >  1:2637299/1‑93 (MQ=255)
                                                                                          aCTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTga  <  2:4552181/93‑1 (MQ=255)
                                                                                          |                                                                                            
GTCCAAATATTACTGTTTATTTATACAGTAAACTTCTATAATATCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCACTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGA  >  NC_000913/3850916‑3851098

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: