Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,240,984 (G)7→6 coding (1041/1245 nt) sstT → serine/threonine:Na(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,240,9780G.92.5% 936.5 / 70.7 212coding (1035/1245 nt)sstTserine/threonine:Na(+) symporter
Reads supporting (aligned to +/- strand):  ref base G (16/0);  new base . (36/160);  total (52/160)
Fisher's exact test for biased strand distribution p-value = 2.33e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

Only 100 of 267 total aligned reads displayed.

TGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGC  >  NC_000913/3240878‑3241074
                                                                                                    |                                                                                                
 ggctgCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCgg                                                                                                 >  2:5891678/1‑101 (MQ=255)
  gctgCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCggg                                                                                                <  1:5891587/101‑1 (MQ=255)
  gctgCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCggg                                                                                                <  2:5891220/101‑1 (MQ=255)
  gctgCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCggg                                                                                                >  1:5891679/1‑101 (MQ=255)
  gctgCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCggg                                                                                                <  2:5891047/101‑1 (MQ=255)
   ctgCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCgggg                                                                                               >  1:5891684/1‑101 (MQ=255)
    tgCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCggggg                                                                                              <  1:5891530/101‑1 (MQ=255)
    tgCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCggggg                                                                                              >  1:5891686/1‑101 (MQ=255)
       ggTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCggggggtc                                                                                           <  2:5891326/101‑3 (MQ=255)
       ggTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCggggggtc                                                                                           <  1:5891170/101‑3 (MQ=255)
        gTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCggggggtct                                                                                          <  1:5891393/101‑4 (MQ=255)
        gTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCggggggtct                                                                                          <  1:5891083/101‑4 (MQ=255)
        gTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCgggggggct                                                                                          >  2:5891691/1‑99 (MQ=255)
         ttAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCggggggtctc                                                                                         <  2:5891369/101‑5 (MQ=255)
          tAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTctct                                                                                       >  1:5891695/1‑101 (MQ=255)
             tACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCg                                                                                                  >  2:5891701/1‑88 (MQ=255)
             tACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTctgct                                                                                    <  1:5891375/101‑1 (MQ=255)
               cGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTctgctgc                                                                                  <  2:5891342/101‑1 (MQ=255)
                gCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTctgctgct                                                                                 <  1:5891242/101‑1 (MQ=255)
                gCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTctgctg                                                                                   <  2:5891364/99‑1 (MQ=255)
                 cTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTctgctgctg                                                                                <  2:5891184/101‑1 (MQ=255)
                  tGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGa                                                                               <  1:5891196/101‑1 (MQ=255)
                   gggTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGAt                                                                              <  2:5891061/101‑1 (MQ=255)
                    ggTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATc                                                                             <  1:5890979/101‑1 (MQ=255)
                    ggTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATc                                                                             <  2:5891125/101‑1 (MQ=255)
                    ggTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATc                                                                             <  2:5890974/101‑1 (MQ=255)
                     gTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATcc                                                                            >  2:5891708/1‑101 (MQ=255)
                        ttCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGCTCTGCTGCTGATCCCAc                                                                         >  1:5891711/1‑101 (MQ=255)
                        ttCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCAc                                                                         <  1:5891268/101‑1 (MQ=255)
                            ggTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCgg                                                                                                 >  1:5891716/1‑74 (MQ=255)
                            ggTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGCTCTGCTGCTGATccc                                                                           >  2:5891714/1‑95 (MQ=255)
                            ggTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGc                                                                     <  1:5891019/101‑1 (MQ=255)
                              tCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCt                                                                   <  1:5891440/101‑1 (MQ=255)
                              tCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCt                                                                   <  2:5891310/101‑1 (MQ=255)
                               cGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGCTCTGCTGCTGATCCCACTGGcc                                                                    >  2:5891718/1‑99 (MQ=255)
                               cGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTg                                                                  <  2:5891320/101‑1 (MQ=255)
                               cGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTg                                                                  <  1:5891288/101‑1 (MQ=255)
                                gATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGCTCTGCTGCTGATCCCACTGGCCTGt                                                                 >  2:5891721/1‑101 (MQ=255)
                                 aTCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTg                                                                       >  1:5891727/1‑94 (MQ=255)
                                 aTCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTa                                                                <  1:5891323/101‑1 (MQ=255)
                                 aTCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTa                                                                <  2:5891205/101‑1 (MQ=255)
                                  tCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTaa                                                               >  1:5891730/1‑101 (MQ=255)
                                    tGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAata                                                             <  1:5891345/101‑1 (MQ=255)
                                     gCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCggggg                                                                                              <  2:5891106/68‑1 (MQ=255)
                                     gCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAatat                                                            <  1:5891500/101‑1 (MQ=255)
                                     gCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAatat                                                            <  2:5891169/101‑1 (MQ=255)
                                       ccACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGt                                                          <  2:5891540/101‑1 (MQ=255)
                                        cACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGC‑‑GGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTc                                                        <  1:5891298/101‑1 (MQ=255)
                                         aCGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTc                                                        >  2:5891739/1‑101 (MQ=255)
                                           ggCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCgg                                                      >  1:5891742/1‑101 (MQ=255)
                                           ggCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCgg                                                      <  1:5891176/101‑1 (MQ=255)
                                             cgcTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCgggggggct                                                                                          >  2:5891745/1‑62 (MQ=255)
                                             cgcTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTa                                                    <  2:5891279/101‑1 (MQ=255)
                                              gctgctGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTgg                                                                      <  2:5891252/82‑1 (MQ=255)
                                               ctgctgTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTAtt                                                  <  2:5891386/101‑1 (MQ=255)
                                               ctgctgTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTAtt                                                  <  1:5891221/101‑1 (MQ=255)
                                               ctgctgTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTAtt                                                  <  1:5891232/101‑1 (MQ=255)
                                                 gctgTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTc                                                >  1:5891748/1‑101 (MQ=255)
                                                  ctgTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGCTCTGCTGCTGATCCCACTGGCCTGt                                                                 >  1:5891751/1‑83 (MQ=255)
                                                  ctgTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCg                                               <  1:5891366/101‑1 (MQ=255)
                                                  ctgTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCg                                               <  1:5891384/101‑1 (MQ=255)
                                                   tgTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGa                                              <  1:5891450/101‑1 (MQ=255)
                                                    gTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGaa                                             <  1:5891135/101‑1 (MQ=255)
                                                      tGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACg                                           >  1:5891755/1‑101 (MQ=255)
                                                       gAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCgggggggc                                                                                           >  2:5891758/1‑52 (MQ=255)
                                                        aGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGat                                         <  1:5891107/101‑1 (MQ=255)
                                                           gtggtgGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATc                                      <  1:5891311/101‑1 (MQ=255)
                                                             ggtggCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGc                                    <  1:5891258/101‑1 (MQ=255)
                                                             ggtggCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGc                                    <  1:5891359/101‑1 (MQ=255)
                                                              gtggCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGcc                                   <  1:5891278/101‑1 (MQ=255)
                                                               tggCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCa                                  >  2:5891762/1‑101 (MQ=255)
                                                                 gCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATg                                >  1:5891764/1‑101 (MQ=255)
                                                                  cTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGc                                    <  2:5891351/96‑1 (MQ=255)
                                                                   ttCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCa                              >  1:5891767/1‑101 (MQ=255)
                                                                     ctctGTGTGCCTGTGGCGCATCCGGCGTGGCggggg                                                                                              >  2:5891771/1‑36 (MQ=255)
                                                                     ctctGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAgg                            <  2:5891733/101‑1 (MQ=255)
                                                                       ctGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAggtg                          >  1:5891772/1‑101 (MQ=255)
                                                                       ctGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAggtg                          <  2:5891495/101‑1 (MQ=255)
                                                                             gCCTGTGGCGCATCCGGCGTGGCGGGGGGGCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGcc                    >  1:5891775/1‑101 (MQ=255)
                                                                              ccTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCg                   <  2:5891116/101‑1 (MQ=255)
                                                                               cTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGt                  <  2:5891233/101‑1 (MQ=255)
                                                                               cTGTGGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGt                  <  1:5891560/101‑1 (MQ=255)
                                                                                tgtgGCGCATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTc                 <  2:5891391/101‑1 (MQ=255)
                                                                                    gcgcATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCt             >  2:5891779/1‑101 (MQ=255)
                                                                                     cgcATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCtt            <  1:5891398/101‑1 (MQ=255)
                                                                                     cgcATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCtt            <  1:5891155/101‑1 (MQ=255)
                                                                                      gcATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCttt           <  2:5891556/101‑1 (MQ=255)
                                                                                      gcATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCttt           <  1:5891306/101‑1 (MQ=255)
                                                                                       cATCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTa          <  2:5891446/101‑1 (MQ=255)
                                                                                        aTCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTat         >  1:5891782/1‑101 (MQ=255)
                                                                                        aTCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTat         <  1:5891487/101‑1 (MQ=255)
                                                                                        aTCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTat         >  1:5891784/1‑101 (MQ=255)
                                                                                        aTCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTat         <  1:5891337/101‑1 (MQ=255)
                                                                                         tCCGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTc                 <  1:5891293/92‑1 (MQ=255)
                                                                                          ccGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTatca       <  2:5891430/101‑1 (MQ=255)
                                                                                          ccGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTatca       <  1:5891461/101‑1 (MQ=255)
                                                                                          ccGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTatca       <  1:5891415/101‑1 (MQ=255)
                                                                                          ccGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTatca       <  2:5891192/101‑1 (MQ=255)
                                                                                           cGGCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTatcat      <  2:5891331/101‑1 (MQ=255)
                                                                                            ggCGTGGCGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTatcatc     <  2:5891416/101‑1 (MQ=255)
                                                                                                    |                                                                                                
TGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGC  >  NC_000913/3240878‑3241074

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