Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,983,543 (C)5→4 coding (609/990 nt) araF ← L‑arabinose‑binding periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,983,5390C.97.6% 914.4 / 9.5 208coding (613/990 nt)araFL‑arabinose‑binding periplasmic protein
Reads supporting (aligned to +/- strand):  ref base C (1/4);  new base . (100/103);  total (101/107)
Fisher's exact test for biased strand distribution p-value = 3.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.34e-01

Only 100 of 228 total aligned reads displayed.

GCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTA  >  NC_000913/1983447‑1983627
                                                                                            |                                                                                        
gccgccCAGCACGGTGCTGTCGTTCATACCGCCGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGAAGCGTCAAATGc                                                                                          >  2:1533179/1‑93 (MQ=255)
 ccgccCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGcc                                                                                         >  2:570833/1‑93 (MQ=255)
  cgccCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGACTTGGCAGCGTCAAATGccc                                                                                        <  2:6129183/93‑1 (MQ=255)
   gccCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGcccc                                                                                       <  2:2095934/93‑1 (MQ=255)
   gccCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGcccc                                                                                       <  1:1581297/93‑1 (MQ=255)
   gccCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGcccc                                                                                       <  2:3093610/93‑1 (MQ=255)
   gccCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGcccc                                                                                       <  2:4354360/93‑1 (MQ=255)
   gccCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGcccc                                                                                       <  2:4749766/93‑1 (MQ=255)
    cccAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGcgccg                                                                                      >  1:1111675/1‑92 (MQ=255)
     ccAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGccccgg                                                                                     <  2:5980589/93‑3 (MQ=255)
       aGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCggga                                                                                   <  1:3701874/93‑2 (MQ=255)
        gCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGAt                                                                                 <  2:2441413/93‑1 (MQ=255)
          ccGGTTCTGTCTTCCATGCCGTCGATCAGCCAATGTTTAACTGCCGGATGTTGAACGCGCATTGAGTCGGCAGCGTCATATGCGCCGGGATGt                                                                               <  2:4004760/92‑1 (MQ=255)
          aCGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGt                                                                               >  1:1298916/1‑93 (MQ=255)
            ggTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCg                                                                             >  2:2046018/1‑93 (MQ=255)
             gTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGt                                                                            >  1:1016596/1‑93 (MQ=255)
             gTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGATGGCAGCGTCAAATGCCCCGGGATGTCGt                                                                            <  2:4496491/93‑1 (MQ=255)
               gCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTa                                                                          <  2:133200/93‑1 (MQ=255)
               gCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTa                                                                          <  2:3422280/93‑1 (MQ=255)
               gCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTa                                                                          <  1:4931669/93‑1 (MQ=255)
                 tGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGTGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGa                                                                        >  2:3140076/1‑93 (MQ=255)
                     gatCATACCGACGATGAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTg                                                                    <  2:4194029/91‑1 (MQ=255)
                     gTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTg                                                                    <  1:5704330/93‑1 (MQ=255)
                      ttCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTgg                                                                   >  1:775763/1‑93 (MQ=255)
                        cATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTggta                                                                 <  2:1856909/93‑1 (MQ=255)
                         aTACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTggtag                                                                <  2:5551196/93‑1 (MQ=255)
                          tACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTggtagg                                                               >  2:3493167/1‑93 (MQ=255)
                            ccGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTggtaggta                                                             <  1:2390754/93‑1 (MQ=255)
                             cgacgaTCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTAc                                                            >  1:5055027/1‑93 (MQ=255)
                              gacgaTCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTAcc                                                           >  1:4373009/1‑93 (MQ=255)
                              gacgaTCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTAcc                                                           >  1:5534155/1‑93 (MQ=255)
                              gacgaTCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTAcc                                                           <  1:4761473/93‑1 (MQ=255)
                               acgaTCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCt                                                          <  2:3531807/93‑1 (MQ=255)
                               acgaTCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCt                                                          >  1:5930701/1‑93 (MQ=255)
                               acgaTCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCt                                                          >  1:2098645/1‑93 (MQ=255)
                                cgaTCAGCAAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTg                                                         <  1:5678432/93‑1 (MQ=255)
                                   tCAGCGAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATa                                                      <  2:5869906/93‑1 (MQ=255)
                                   tCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATa                                                      <  1:973967/93‑1 (MQ=255)
                                   tCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATa                                                      <  1:3212334/93‑1 (MQ=255)
                                   tCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATa                                                      >  1:4203017/1‑93 (MQ=255)
                                    cAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATaa                                                     >  2:5373069/1‑93 (MQ=255)
                                     aGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATaaa                                                    >  1:483805/1‑93 (MQ=255)
                                      gCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAAt                                                   >  2:3291912/1‑93 (MQ=255)
                                      gCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAAt                                                   >  1:3254880/1‑93 (MQ=255)
                                      gCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAAt                                                   <  2:1049057/93‑1 (MQ=255)
                                       ccAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAAtt                                                  <  1:5223273/93‑1 (MQ=255)
                                        cAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAAttt                                                 >  1:2127885/1‑93 (MQ=255)
                                        cAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAAttt                                                 <  2:1447891/93‑1 (MQ=255)
                                          aTGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGt                                               >  2:1058470/1‑93 (MQ=255)
                                           tGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGtt                                              >  1:3290039/1‑93 (MQ=255)
                                            gTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGttt                                             >  2:276635/1‑93 (MQ=255)
                                            gTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGttt                                             <  2:1405796/93‑1 (MQ=255)
                                             tttAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGtttt                                            <  2:5464158/93‑1 (MQ=255)
                                             tttAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGtttt                                            <  1:1730361/93‑1 (MQ=255)
                                             tttAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGtttt                                            >  1:4697134/1‑93 (MQ=255)
                                              ttAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGttttt                                           >  2:747600/1‑93 (MQ=255)
                                              ttAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGttttt                                           >  1:2170652/1‑93 (MQ=255)
                                               tAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTc                                          <  2:1576234/93‑1 (MQ=255)
                                               tAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTc                                          <  2:2261903/93‑1 (MQ=255)
                                                aaCTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTcc                                         <  1:4521610/93‑1 (MQ=255)
                                                  cTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCgg                                       <  2:6279028/93‑1 (MQ=255)
                                                  cTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCgg                                       <  2:939680/93‑1 (MQ=255)
                                                   ttCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCggg                                      >  2:13066/1‑93 (MQ=255)
                                                   ttCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCggg                                      >  1:4978833/1‑93 (MQ=255)
                                                   ttCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCggg                                      <  2:4315239/93‑1 (MQ=255)
                                                     ccGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGaa                                    >  1:2659552/1‑93 (MQ=255)
                                                      cGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAAt                                   >  2:4221570/1‑93 (MQ=255)
                                                      cGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAAt                                   <  2:1602027/93‑1 (MQ=255)
                                                       ggATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATc                                  >  1:3752141/1‑93 (MQ=255)
                                                        gATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGCTAAATTTGTTTTTCCGGGAATcg                                 >  2:1763869/1‑92 (MQ=255)
                                                        gATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATcc                                 <  2:2718554/93‑1 (MQ=255)
                                                        gATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATcc                                 >  1:5171631/1‑93 (MQ=255)
                                                          tGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGCTAAATTTGTTTTTCCGGGAATCCgc                               >  1:198816/1‑92 (MQ=255)
                                                           gTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGc                              >  2:2325832/1‑93 (MQ=255)
                                                                aCCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCtt                         >  1:1086916/1‑93 (MQ=255)
                                                                 ccAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCttt                        <  1:2496110/93‑1 (MQ=255)
                                                                  cAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTc                       >  1:5210191/1‑93 (MQ=255)
                                                                  cAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTc                       >  2:2620425/1‑93 (MQ=255)
                                                                   aGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCa                      <  1:3525802/93‑1 (MQ=255)
                                                                   aGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCa                      >  1:1073580/1‑93 (MQ=255)
                                                                   aGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCa                      >  1:692468/1‑93 (MQ=255)
                                                                   aGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCa                      <  2:2885181/93‑1 (MQ=255)
                                                                     cATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAgc                    <  1:4491868/93‑1 (MQ=255)
                                                                     cATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAgc                    <  2:5920695/93‑1 (MQ=255)
                                                                      aTTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAgcg                   >  2:5792588/1‑93 (MQ=255)
                                                                          aGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATc               <  1:3333858/93‑1 (MQ=255)
                                                                           gTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTGCGGGAATCCGGCCGCTTTCAGCGCATca              >  1:1460556/1‑92 (MQ=255)
                                                                            ttGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCa             <  1:5994779/93‑1 (MQ=255)
                                                                              ggCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATa           <  1:190526/93‑1 (MQ=255)
                                                                              ggCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATa           <  1:3039803/93‑1 (MQ=255)
                                                                               gCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAg          >  1:5197724/1‑93 (MQ=255)
                                                                                cAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGa         >  1:5521247/1‑93 (MQ=255)
                                                                                  gCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATc       <  1:5426719/93‑1 (MQ=255)
                                                                                  gCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATc       >  1:1959050/1‑93 (MQ=255)
                                                                                  gCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATc       <  2:2025039/93‑1 (MQ=255)
                                                                                   cGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATcc      >  2:4417107/1‑93 (MQ=255)
                                                                                   cGTAAAATGCCCCGGAATGTCGTTAGATTTAGTAGGTAACTGATAAATTGGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATcc      <  2:4931863/93‑1 (MQ=255)
                                                                                    gTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATccc     <  1:4585197/93‑1 (MQ=255)
                                                                                    gTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATccc     <  2:2037620/93‑1 (MQ=255)
                                                                                       aaaTGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCgta  <  2:4639600/93‑1 (MQ=255)
                                                                                            |                                                                                        
GCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTA  >  NC_000913/1983447‑1983627

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