Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 867,718 G→C 7.5% E67Q (GAA→CAA)  gsiA → glutathione ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913867,7180GC 7.5% 313.4 / 13.1 92E67Q (GAA→CAA) gsiAglutathione ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (44/41);  new base C (4/3);  total (48/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.65e-01

TGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACAT  >  NC_000913/867572‑867865
                                                                                                                                                  |                                                                                                                                                   
tGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGaa                                                                                                                                                   <  4:224481/149‑1 (MQ=255)
    tatTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAgg                                                                                                                                               >  3:714087/1‑149 (MQ=255)
     atTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAgg                                                                                                                                               <  1:784793/148‑1 (MQ=255)
     atTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTCCAACAggc                                                                                                                                              <  4:131779/149‑1 (MQ=255)
      ttGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCGATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAggcg                                                                                                                                             >  3:536547/1‑149 (MQ=255)
          cTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAggcgggcg                                                                                                                                         >  4:317378/1‑149 (MQ=255)
           tttATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAggcgggcg                                                                                                                                         >  2:123216/1‑148 (MQ=255)
             tATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGtt                                                                                                                                      <  4:529220/149‑1 (MQ=255)
               tGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGttt                                                                                                                                     <  1:392423/148‑1 (MQ=255)
                gCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAg                                                                                                                                   >  4:132651/1‑149 (MQ=255)
                 cAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAg                                                                                                                                   >  2:691551/1‑148 (MQ=255)
                 cAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAg                                                                                                                                   >  2:310816/1‑148 (MQ=255)
                   ggACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTAc                                                                                                                                >  4:690605/1‑149 (MQ=255)
                     aCCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACa                                                                                                                               >  1:341265/1‑148 (MQ=255)
                     aCCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACa                                                                                                                               >  1:600269/1‑148 (MQ=255)
                      ccAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAg                                                                                                                              <  1:691551/148‑1 (MQ=255)
                         gcagAAAATCGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCg                                                                                                                          <  3:769806/149‑1 (MQ=255)
                            gAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGAt                                                                                                                        >  1:907016/1‑148 (MQ=255)
                                 tAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAAt                                                                                                                   >  2:346232/1‑148 (MQ=255)
                                     tGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGt                                                                                                              <  3:875669/149‑1 (MQ=255)
                                      gCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGt                                                                                                              <  2:803413/148‑1 (MQ=255)
                                           ccGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAgc                                                                                                        <  3:267783/149‑1 (MQ=255)
                                             gCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGt                                                                                                                             >  2:233220/1‑126 (MQ=255)
                                             gCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGt                                                                                                                             <  1:233220/126‑1 (MQ=255)
                                             gCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAgcg                                                                                                       >  1:618987/1‑148 (MQ=255)
                                              cAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAgcggc                                                                                                     <  4:806943/149‑1 (MQ=255)
                                              cAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAgcgg                                                                                                      <  1:85333/148‑1 (MQ=255)
                                                 tctctcTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCa                                                                                                  >  3:324324/1‑149 (MQ=255)
                                                   tctctTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGt                                                                                                >  3:379780/1‑149 (MQ=255)
                                                    ctctTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGt                                                                                                <  2:258638/148‑1 (MQ=255)
                                                        tttAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTcgcg                                                                                            >  1:399145/1‑148 (MQ=255)
                                                         ttAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGaa                                                                                          <  3:132651/149‑1 (MQ=255)
                                                         ttAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGa                                                                                           <  1:508630/148‑1 (MQ=255)
                                                              cTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGa                                                                                      >  1:705119/1‑148 (MQ=255)
                                                              cTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGAt                                                                                     >  2:123140/1‑148 (MQ=255)
                                                                  aaCGCGGTGAGATGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGa                                                                                  <  1:721137/148‑1 (MQ=255)
                                                                  aaCGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGc                                                                                                                 <  1:48915/117‑1 (MQ=255)
                                                                  aaCGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGc                                                                                                                 >  2:48915/1‑117 (MQ=255)
                                                                  aaCGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGaa                                                                                 <  3:235184/149‑1 (MQ=255)
                                                                  aaCGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGaa                                                                                 <  4:556759/149‑1 (MQ=255)
                                                                  aaCGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGa                                                                                  <  1:110457/148‑1 (MQ=255)
                                                                       ggTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAg                                                                            <  3:690605/149‑1 (MQ=255)
                                                                        gTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAg                                                                            <  2:99147/148‑1 (MQ=255)
                                                                          gagaCGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCg                                                                          <  1:748433/148‑1 (MQ=255)
                                                                            gaCGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAg                                                                        <  1:310816/148‑1 (MQ=255)
                                                                             aCGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGc                                                                       <  1:210048/148‑1 (MQ=255)
                                                                                 tGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGaa                                                                   <  2:468771/148‑1 (MQ=255)
                                                                                   gCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGc                                                                >  4:185180/1‑149 (MQ=255)
                                                                                    cAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGc                                                                <  2:792843/148‑1 (MQ=255)
                                                                                     aaTTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTg                                                              <  4:379780/149‑1 (MQ=255)
                                                                                         gttgGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGcac                                                           >  2:311612/1‑148 (MQ=255)
                                                                                          ttgGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACaa                                                         >  4:790961/1‑149 (MQ=255)
                                                                                           tgGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACaaa                                                        <  3:548217/149‑1 (MQ=255)
                                                                                            ggCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACaaa                                                        >  1:155853/1‑148 (MQ=255)
                                                                                            ggCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACaaa                                                        >  1:858229/1‑148 (MQ=255)
                                                                                              cGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATg                                                      >  2:898807/1‑148 (MQ=255)
                                                                                               gAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATgcg                                                    >  3:835636/1‑149 (MQ=255)
                                                                                                     ggCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATg                                               >  2:476175/1‑148 (MQ=255)
                                                                                                      gCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGtt                                             >  4:359507/1‑149 (MQ=255)
                                                                                                         ccGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTcg                                           <  2:318330/148‑1 (MQ=255)
                                                                                                            gTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCgg                                        >  1:114776/1‑148 (MQ=255)
                                                                                                              aaGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGc                                     <  4:744269/149‑1 (MQ=255)
                                                                                                                    gTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGatat                               >  4:523848/1‑149 (MQ=255)
                                                                                                                    gTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGata                                >  2:216748/1‑148 (MQ=255)
                                                                                                                        cTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGc                            >  1:69101/1‑148 (MQ=255)
                                                                                                                         tGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCg                           <  2:618987/148‑1 (MQ=255)
                                                                                                                          gCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCga                          >  1:709786/1‑148 (MQ=255)
                                                                                                                            gTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCgatga                       >  4:88613/1‑149 (MQ=255)
                                                                                                                             ttGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGa                                                                         <  1:499977/98‑1 (MQ=255)
                                                                                                                             ttGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGa                                                                         >  2:499977/1‑98 (MQ=255)
                                                                                                                             ttGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCgatga                       >  1:834051/1‑148 (MQ=255)
                                                                                                                                gCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATAtt                   >  3:78197/1‑149 (MQ=255)
                                                                                                                                gCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATAtt                   >  4:323315/1‑149 (MQ=255)
                                                                                                                                  aTTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTc                 >  3:14485/1‑149 (MQ=255)
                                                                                                                                   ttGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCgatg                        <  1:93915/141‑1 (MQ=255)
                                                                                                                                   ttGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGCCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCa                <  4:536547/149‑1 (MQ=255)
                                                                                                                                   ttGATGCGCCTGTTGGAACAGGCGGGCGGATTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCgatg                        >  2:93915/1‑141 (MQ=255)
                                                                                                                                    tGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCa                <  2:907016/148‑1 (MQ=255)
                                                                                                                                      aTGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGa             <  4:324324/149‑1 (MQ=255)
                                                                                                                                        ccgcCTGTTCCAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGc           >  3:835023/2‑149 (MQ=255)
                                                                                                                                         cgcCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGcc          <  3:185180/149‑1 (MQ=255)
                                                                                                                                         cgcCTGTTCCAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACTCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGc           >  2:125886/1‑148 (MQ=255)
                                                                                                                                         cgcCTGTTCCAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGc           >  2:424935/1‑148 (MQ=255)
                                                                                                                                          gcCTGTTCCAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGc                                     <  4:421744/121‑1 (MQ=255)
                                                                                                                                          gcCTGTTCCAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGc                                     >  3:421744/1‑121 (MQ=255)
                                                                                                                                           ccTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATAtt                   <  1:245568/138‑1 (MQ=255)
                                                                                                                                           ccTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATAtt                   >  2:245568/1‑138 (MQ=255)
                                                                                                                                            cTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGAt       >  4:466950/1‑149 (MQ=255)
                                                                                                                                            cTGTTCCAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGa        <  1:424935/148‑1 (MQ=255)
                                                                                                                                             tGTTGGAACAGGCGGGCGGTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATg      <  1:123140/148‑1 (MQ=255)
                                                                                                                                               ttGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGa     <  2:705119/148‑1 (MQ=255)
                                                                                                                                                tGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACa   >  3:842900/1‑149 (MQ=255)
                                                                                                                                                 ggAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACAt  >  3:826841/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACAT  >  NC_000913/867572‑867865

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