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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,240,9771.G95.0% 212.8 / 3.9 60coding (1034/1245 nt)sstTsodium:serine/threonine symporter
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base G (28/29);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATCAATATGGCGGGCGCAGCAATCACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGCGTATTGCAGGACTCTTGCGAAACCGCGCTGAACTCTTCAACTGACG  >  NC_000913/3240833‑3241120
                                                                                                                                                 |                                                                                                                                                  
ccaTCAATATGGCGGGCGCAGCAATCACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGTCgggggg                                                                                                                                              <  2:205623/151‑1 (MQ=255)
    cAATATGGCGGGAGCGCCAAAAACTATTACCGTTGTGGCGCTGGCTGCGGTTGATATGATGGGTATTTCTGTCGACTTGCCCACCGCGCTGCCGTTTGGAGGGGGGGGTTCTCTGTGTGCCTGTGGCGCATCCCGCCTggggggggggggg                                                                                                                                          >  2:179993/1‑148 (MQ=255)
        atGGCGGGCGCAGCAATCACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTctg                                                                                                                                      <  1:191196/151‑1 (MQ=255)
          ggcgggcgCAGCAATCACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTctgct                                                                                                                                    <  1:580639/151‑1 (MQ=255)
          ggcgggcgCAGCAATCACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTctgc                                                                                                                                     <  2:119682/150‑1 (MQ=255)
                   agcaATCACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGCTTTTCTGGTGATccc                                                                                                                          >  2:462994/1‑151 (MQ=255)
                   agcaATCACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGCCCTGCTGCTGATccc                                                                                                                          >  1:378216/1‑151 (MQ=255)
                     caATCACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCAc                                                                                                                        <  2:637828/151‑1 (MQ=255)
                     caATCACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCAc                                                                                                                        <  1:747189/151‑1 (MQ=255)
                       aTCACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTg                                                                                                                      <  2:536818/151‑1 (MQ=255)
                               tACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGCTCTGCTGGTGATCCCCCTGGCCTGt                                                                                                                >  2:805106/1‑149 (MQ=255)
                                    tgtTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGCTCTGCTGCTGATACCACCTGGCTGTaa                                                                                                              >  2:908585/1‑146 (MQ=255)
                                    tgtTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGCCCTGCTGCTGAACCCACCGGCCGGTTatat                                                                                                           >  2:955447/1‑149 (MQ=255)
                                      ttGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGCTCTGCTGCTGGTCCCACTGGCCTGTAATAtttt                                                                                                        >  1:95336/1‑150 (MQ=255)
                                       tGACGCTGGCTGCGGTTGATACGCTGGGTATTCCGGTCGATCTGCCCTCGGCGCTGCTGTTGAGCGTGGTGGCTTCGCTGTGTGCCTGTGGCGCATACGGCTTGGCGGGGGGGCTCTCTGCTGCTGATCCCACTAGCCTGTAatat                                                                                                           <  2:187499/146‑1 (MQ=255)
                                        gACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCgg                                                                                                     <  1:277641/151‑1 (MQ=255)
                                        gACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCgg                                                                                                     <  2:111234/151‑1 (MQ=255)
                                        gACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCgg                                                                                                     <  2:559287/150‑1 (MQ=255)
                                           gctggctgCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCCGtct                                                                                                  >  1:603364/1‑149 (MQ=255)
                                             tggctgCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGCTCTGCTGCTGATCCCACCGGCCTGTAATATGTTCGGTAttt                                                                                                >  2:1069256/1‑151 (MQ=255)
                                                ctgCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTTGGTATTTCg                                                                                              >  1:9474/1‑150 (MQ=255)
                                                  gCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGCTCTGCTGCTGATCCCACCGGCCTGTAATATGTTTGGGATTTCgga                                                                                            >  1:1044293/1‑148 (MQ=255)
                                                    ggTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGCTCTGCTGCTGATTCCACTGGCCTGTAATATGTTTGGTATTTCGAAc                                                                                           >  1:948620/1‑149 (MQ=255)
                                                    ggTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGGATTTCGAACg                                                                                          >  1:423060/1‑150 (MQ=255)
                                                        aaTACGCTGGGTATTCCGTTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGatat                                                                                      <  2:173657/150‑1 (MQ=255)
                                                        aaTACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGatat                                                                                      <  1:379212/150‑1 (MQ=255)
                                                         aTACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGCTCTGCTGCTGATCCCCCTGGCCTGTAATATGTTCGGGATTTCGAACGatat                                                                                      >  2:775264/1‑149 (MQ=255)
                                                            cGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCa                                                                                 <  2:972643/151‑1 (MQ=255)
                                                            cGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGCTCTTCTTCTGATTCCCCTTGCCTGTCAAATGTTTGGGATTTTGCACGGTATccgc                                                                                  >  2:824951/1‑147 (MQ=255)
                                                              cTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGCTCTCCTGCTGATCACACTGGCCTTTAATATTTTCGGGATTTTGAACGCTATCGCCa                                                                                 >  2:640364/1‑149 (MQ=255)
                                                               tGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTTGCCTGTAATATGTTCGGGATTTCGAACGATATCGcc                                                                                  >  1:564720/1‑147 (MQ=255)
                                                                  gTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGCTCTGCCGGTGGTCCCCCTTGCCTCTCATATGTTCGGGATTTCGAACGATTTCGCCCTg                                                                               >  2:229222/1‑147 (MQ=255)
                                                                    aTTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAggt                                                                          <  1:787748/150‑1 (MQ=255)
                                                                    aTTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAggt                                                                          <  1:779993/150‑1 (MQ=255)
                                                                     ttCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCCCTGGCCTGTAATATGTTCGGGATTTCGAACGATATCGCCCTGCAggtgg                                                                        >  2:856610/1‑151 (MQ=255)
                                                                       ccGGTGGGTCTGCCCACGGCGCTGCTGTGGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCGTCCGGCGTGGCGGGGGGGGTGTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGtt                                                                      <  2:216057/151‑1 (MQ=255)
                                                                         ggTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAata                                                                                                            <  2:432100/111‑1 (MQ=255)
                                                                           tCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGCTCTGCTGCTGGTCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGGGGGTgtc                                                                   >  2:1064959/1‑148 (MQ=255)
                                                                            cGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACCGGCCTGTAATATGTTCGGGATTTCGAACGATATCGCCATGCAGGTGGTTGCCg                                                                  >  1:806451/1‑150 (MQ=255)
                                                                               tCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGCTCTGCTGCTGATTCCACTGGCCTGTAATATGTTCGGGATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCgg                                                             >  1:195353/1‑151 (MQ=255)
                                                                                 tGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGCTCTGCTGCTGATCCCACTGGCCTGTTATATGTTCGGGATTTCGAACGATATCGCCCTGCAGGTGGGTGCCGTCgg                                                             >  2:308092/1‑149 (MQ=255)
                                                                                  gCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGGTTCCACTGGGCTGTAATATGTTCGGTATTTCGAACGGTATCGCCCTGCAGGTGGTTGCCGTCGGCtt                                                          >  1:421767/1‑151 (MQ=255)
                                                                                  gCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTc                                                               <  1:377487/146‑1 (MQ=255)
                                                                                  gCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTc                                                               <  2:489625/146‑1 (MQ=255)
                                                                                  gCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCt                                                           <  2:942338/150‑1 (MQ=255)
                                                                                       tggCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTTGGGATTTTGAACGATATCGCCCTGCAGGTGGGTGCCGTCGGCttttt                                                       >  2:46799/2‑147 (MQ=255)
                                                                                          cgcTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCg                                                  <  1:775306/151‑1 (MQ=255)
                                                                                           gctgctGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCgg                                                 <  2:466966/151‑1 (MQ=255)
                                                                                           gctgctGTTGAGCGTGGTGGCTTCGCTGTGGGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGCTCTGGTGCTGATACCAACTGGCTGTTATTTTTTTGGGTTTTTGGACGATATCGGCACGCAGGGGGGTGGCGTTCGGCttatca                                                      >  2:463037/1‑143 (MQ=255)
                                                                                             tgctgTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGCTCTGCTGCTGATCCCACTGGGCTGTAATATGTTCGGGATTTTGAAAGATATCGCCATGCAGGTGGGTGGCGTCGGCTTTttcatctggg                                               >  2:73434/1‑147 (MQ=255)
                                                                                              gctgTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGc                                                <  1:908632/149‑1 (MQ=255)
                                                                                              gctgTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGCGt                                              <  1:313852/151‑1 (MQ=255)
                                                                                                     aGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGCGTATTGCAg                                       <  1:638109/151‑1 (MQ=255)
                                                                                                     aGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGCGTATTGCAg                                       <  1:446484/151‑1 (MQ=255)
                                                                                                                 tctcTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCCCTGGCCTGTAATATGTTCGGTATTTCGAGCGATATCGCCATGCAGGTGGTTGCCGTCGGCGTTATCATCGGCGTATTGCAGGACTCTTg                               <  2:86336/147‑1 (MQ=255)
                                                                                                                            cTGTGGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTTGAAAGATATCGCCATGCAGGGGGGTGCCGTCGGCTTTTTCATCCGCGTATTGCAGGAATCTTGCGaaa                          >  1:683679/1‑141 (MQ=255)
                                                                                                                              gtgGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGCGTATTGCAGGACTCTTGCGAAACCGCGCTGaa                <  1:803814/149‑1 (MQ=255)
                                                                                                                               tgGCGCATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGGATTTCGAACGATATCGCCATGCAGGTGGGTGCCGTCGGCTTTTTTATCGGCGTATTGCAGGACCCTTGCGAAAACGCGCTGaa                >  1:951703/1‑148 (MQ=255)
                                                                                                                                  cgcATCCGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGCGTATTGCAGGACTCTTGCGAAACCGCGCTGAACTCtt           <  2:195350/150‑1 (MQ=255)
                                                                                                                                       ccGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGCGTATTGCAGGACTCTTGCGAAACCGCGCTGAACTCTTCAACt      <  2:475445/150‑1 (MQ=255)
                                                                                                                                        cGGCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGCGTATTGCAGGACTCTTGCGAAACCGCGCTGAACTCTTCAACTGa    <  2:344669/151‑1 (MQ=255)
                                                                                                                                          gCGTGGCGGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGCGTATTGCAGGACTCTTGCGAAACCGCGCTGAACTCTTCAACTGACg  <  1:229247/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CCATCAATATGGCGGGCGCAGCAATCACTATTACCGTGTTGACGCTGGCTGCGGTTAATACGCTGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGCGTATTGCAGGACTCTTGCGAAACCGCGCTGAACTCTTCAACTGACG  >  NC_000913/3240833‑3241120

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