Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,984,089 A→G A21A (GCT→GCC araF ← L‑arabinose‑binding periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,984,0890AG99.7% 1060.1 / ‑3.9 303coding (63/990 nt)araFL‑arabinose‑binding periplasmic protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (156/146);  minor base . (0/1);  total (156/147)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.93e-01

Only 100 of 313 total aligned reads displayed.

CCCGGCTTTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATAGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCTCCAGCTTTAGTGTCGTTTT  >  NC_000913/1983999‑1984178
                                                                                          |                                                                                         
 ccGGCTTTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCg                                                                                        >  1:5121510/1‑93 (MQ=255)
  cGGCTTTATCGGCAAACTTCCATTCTGTCTGGAACCACGGCCCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGGTCTCCGCCATGGCgg                                                                                       >  2:3794865/1‑93 (MQ=255)
   ggCTTTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGa                                                                                      >  2:1267727/1‑93 (MQ=255)
   ggCTTTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGa                                                                                      <  2:6722585/93‑1 (MQ=255)
   ggCTTTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGa                                                                                      >  2:2338879/1‑93 (MQ=255)
   ggCTTTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGa                                                                                      <  1:5321513/93‑1 (MQ=255)
    gCTTTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGAt                                                                                     <  2:5041185/93‑1 (MQ=255)
    gCTTTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGAt                                                                                     <  2:8367285/93‑1 (MQ=255)
     cTTTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGAtt                                                                                    >  1:7583672/1‑93 (MQ=255)
     cTTTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGAtt                                                                                    >  1:1755064/1‑93 (MQ=255)
      tttATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATtg                                                                                   <  2:103327/93‑1 (MQ=255)
       ttATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATtgt                                                                                  <  1:1462109/93‑1 (MQ=255)
        tATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATtgtg                                                                                 <  1:8053394/93‑1 (MQ=255)
        tATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATtgtg                                                                                 >  1:5465028/1‑93 (MQ=255)
        tATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATtgtg                                                                                 <  1:8150004/93‑1 (MQ=255)
            ggCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACAt                                                                             <  2:597704/93‑1 (MQ=255)
            ggCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACAt                                                                             <  1:4118775/93‑1 (MQ=255)
            ggCAAACTTCCATTCGGTCTGGAACAACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACAt                                                                             <  1:476999/93‑1 (MQ=255)
             gCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATa                                                                            <  2:4325285/93‑1 (MQ=255)
              cAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATaa                                                                           <  1:4280180/93‑1 (MQ=255)
               aaaCTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAAc                                                                          >  2:7576228/1‑93 (MQ=255)
               aaaCTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAAc                                                                          >  2:5830461/1‑93 (MQ=255)
                aaCTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACg                                                                         <  2:1859039/93‑1 (MQ=255)
                 aCTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACgg                                                                        >  1:4953222/1‑93 (MQ=255)
                  cTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGc                                                                       <  1:1930314/93‑1 (MQ=255)
                        ttCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAg                                                                 >  1:4241682/1‑93 (MQ=255)
                        ttCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAg                                                                 <  2:6471683/93‑1 (MQ=255)
                        ttCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAg                                                                 >  2:3877137/1‑93 (MQ=255)
                          cGGTCTGGAACCACGGCTCTGCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGGTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGAc                                                               <  2:2960635/93‑1 (MQ=255)
                            gTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAGCCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCa                                                             <  2:4836434/93‑1 (MQ=255)
                            gTCTGGAACCACGGCTCTTCCGGATGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCa                                                             >  1:6253011/1‑93 (MQ=255)
                             tCTGGAGCCACGGCTCTTCCGGTTGCTTCACCAGAACACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCaa                                                            <  2:7776801/93‑1 (MQ=255)
                             tCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCaa                                                            <  1:259882/93‑1 (MQ=255)
                             tCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCaa                                                            >  1:6145507/1‑93 (MQ=255)
                              cTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAAt                                                           <  1:1035823/93‑1 (MQ=255)
                               tGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCCATg                                                          >  1:8451437/1‑93 (MQ=255)
                               tGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATg                                                          <  2:1201698/93‑1 (MQ=255)
                               tGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATg                                                          >  2:6099643/1‑93 (MQ=255)
                               tGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATg                                                          <  1:4495154/93‑1 (MQ=255)
                               tGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATg                                                          >  1:6056229/1‑93 (MQ=255)
                               tGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATg                                                          >  2:1585391/1‑93 (MQ=255)
                                 gAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGc                                                        <  1:3338241/93‑1 (MQ=255)
                                 gAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGc                                                        <  2:3411888/93‑1 (MQ=255)
                                 gAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACCGCTGCCCGACCAATGGc                                                        >  2:5666398/1‑93 (MQ=255)
                                   aCCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTg                                                      <  2:5273843/93‑1 (MQ=255)
                                    ccACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGc                                                     >  2:2851981/1‑93 (MQ=255)
                                     cACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGcc                                                    >  2:6208863/1‑93 (MQ=255)
                                     cACGGATCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGGCCAATGCATGcc                                                    >  1:5195868/1‑93 (MQ=255)
                                      aCGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGccg                                                   >  1:2860790/1‑92 (MQ=255)
                                      aCGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGAGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCa                                                   >  2:8179447/1‑93 (MQ=255)
                                       cGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAg                                                  <  2:8005342/93‑1 (MQ=255)
                                       cGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAg                                                  <  1:4333000/93‑1 (MQ=255)
                                        ggCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAgg                                                 <  1:3872180/93‑1 (MQ=255)
                                         gCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAggg                                                <  1:2759211/93‑1 (MQ=255)
                                         gCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAggg                                                >  2:5095134/1‑93 (MQ=255)
                                          ctctTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGc                                               <  2:5646207/93‑1 (MQ=255)
                                           tctTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCCGGGCt                                              >  2:7868102/1‑93 (MQ=255)
                                           tctTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCt                                              >  1:8280965/1‑93 (MQ=255)
                                            ctTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCCGACCAATAGCTGCCAGGGCtt                                             >  1:5688452/1‑93 (MQ=255)
                                             ttCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCttt                                            >  2:2175411/1‑93 (MQ=255)
                                             ttCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCttt                                            >  1:203563/1‑93 (MQ=255)
                                               ccGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAg                                          <  2:5467192/93‑1 (MQ=255)
                                                cGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGt                                         >  1:8751984/1‑93 (MQ=255)
                                                 ggTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTa                                        >  1:3142876/1‑93 (MQ=255)
                                                 ggTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTa                                        >  1:7105473/1‑93 (MQ=255)
                                                 ggTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTa                                        >  1:4389539/1‑93 (MQ=255)
                                                  gTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTaa                                       >  2:3660059/1‑93 (MQ=255)
                                                  gTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTaa                                       >  1:7320158/1‑93 (MQ=255)
                                                    tGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAAt                                     <  1:4041562/93‑1 (MQ=255)
                                                    tGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAAt                                     >  2:7115232/1‑93 (MQ=255)
                                                     gCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAAtt                                    >  1:8408823/1‑93 (MQ=255)
                                                        tCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTtgt                                 >  2:5776689/1‑93 (MQ=255)
                                                        tCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTtgg                                 >  1:7889384/1‑92 (MQ=255)
                                                          aCCAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGc                               <  2:8476533/93‑1 (MQ=255)
                                                            cAGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCAt                             <  2:8091633/93‑1 (MQ=255)
                                                             aGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATg                            >  1:2335270/1‑93 (MQ=255)
                                                             aGAAAACCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATg                            >  2:1520120/1‑93 (MQ=255)
                                                                aaaCCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGtt                         >  2:456306/1‑93 (MQ=255)
                                                                aaaCCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCAAGGtt                         >  2:7770990/1‑93 (MQ=255)
                                                                 aaCCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGttt                        >  1:6339844/1‑92 (MQ=255)
                                                                 aaCCGAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTtc                        <  1:1071540/93‑1 (MQ=255)
                                                                     gAGCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCCGGGCTTTAGTAAATTTGTGCATGGTtctctc                    >  2:7406330/1‑93 (MQ=255)
                                                                       gCTTCAGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCTCCa                  >  1:6900162/1‑93 (MQ=255)
                                                                       cagtCAGGTTCGCAGACATGGCGGATTGTGACATATCGGCTGCCAGACCAATTGCTGCCAGGTCTTTAGTAAATTTGTGCATGGTTCTCGCCa                  <  2:263051/90‑1 (MQ=255)
                                                                          tCAGGTTCACCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCTCCAGCt               <  2:4688114/93‑1 (MQ=255)
                                                                            aGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCTCCAGCttt             >  2:3612033/1‑93 (MQ=255)
                                                                            aGGTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCTCCAGCttt             >  1:6581951/1‑93 (MQ=255)
                                                                              gTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCTCCAGCTTTAg           <  2:6018065/93‑1 (MQ=255)
                                                                              gTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCTCCAGCTTTAg           >  2:693764/1‑93 (MQ=255)
                                                                              gTTCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCTCCAGCTTTAg           <  1:78033/93‑1 (MQ=255)
                                                                               ttCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCTCCAGCTTTAgt          >  1:6166680/1‑93 (MQ=255)
                                                                               ttCTCCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCTCCAGCTTTAgt          <  1:7010669/93‑1 (MQ=255)
                                                                                 ctcCGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCTCCAGCTTTAgtgt        <  2:1338086/93‑1 (MQ=255)
                                                                                   ccGCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCTCCAGCTTTAGTGTCg      >  2:8676050/1‑93 (MQ=255)
                                                                                     gCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCTCCAGCTTTAGTGTCGtt    >  1:3516248/1‑93 (MQ=255)
                                                                                     gCCATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCTCCAGCTTTAGTGTCGtt    <  2:3125978/93‑1 (MQ=255)
                                                                                      ccATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCTCCAGCTTTAGTGTCGttt   <  1:770577/93‑1 (MQ=255)
                                                                                      ccATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCTCCAGCTTTAGTGTCGttt   >  1:3008912/1‑93 (MQ=255)
                                                                                      ccATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCTCCAGCTTTAGTGTCGttt   <  1:1221319/93‑1 (MQ=255)
                                                                                       cATGGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCTCCAGCTTTAGTGTCGtttt  <  1:1497078/93‑1 (MQ=255)
                                                                                          |                                                                                         
CCCGGCTTTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATAGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCTCCAGCTTTAGTGTCGTTTT  >  NC_000913/1983999‑1984178

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