Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 867,718 G→C 9.1% 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 9.1% 303.9 / 18.4 88E67Q (GAA→CAA) gsiAglutathione ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (40/40);  new base C (4/4);  total (44/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 = 7.60e-01

GAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCG  >  NC_000913/867573‑867831
                                                                                                                                                 |                                                                                                                 
gAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGaa                                                                                                                 >  4:599118/1‑148 (MQ=255)
  atatTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGa                                                                                                                  <  1:97633/145‑1 (MQ=255)
  atatTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGa                                                                                                                  >  2:97633/1‑145 (MQ=255)
          tttATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAggcgggcg                                                                                                       >  4:330080/1‑148 (MQ=255)
              tGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGttt                                                                                                   <  3:330080/148‑1 (MQ=255)
                  ggACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTCCAACAgg                                                                                                             >  2:390603/1‑134 (MQ=255)
                  ggACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTCCAACAgg                                                                                                             <  1:390603/134‑1 (MQ=255)
                       agcagAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGtt                                                                                                    >  1:428031/1‑138 (MQ=255)
                       agcagAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGtt                                                                                                    <  2:428031/138‑1 (MQ=255)
                         cagAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAg                                                                                            >  1:49443/1‑144 (MQ=255)
                         cagAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAg                                                                                            <  2:49443/144‑1 (MQ=255)
                         cagAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCg                                                                                        <  3:599118/148‑1 (MQ=255)
                                tAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGttt                                                                                                   >  4:158754/1‑130 (MQ=255)
                                tAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGttt                                                                                                   <  3:158754/130‑1 (MQ=255)
                                            gCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAAt                                                                                 >  1:1091887/1‑136 (MQ=255)
                                            gCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAAt                                                                                 <  2:1091887/136‑1 (MQ=255)
                                                tctctcTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGttt                                                                                                   >  4:1121422/1‑114 (MQ=255)
                                                tctctcTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGttt                                                                                                   <  3:1121422/114‑1 (MQ=255)
                                                   ctctTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAg                                                                                            <  3:1019990/118‑1 (MQ=255)
                                                   ctctTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAg                                                                                            >  4:1019990/1‑118 (MQ=255)
                                                     ctTTTAGTCTGCAAGGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTcg                                                            >  2:774553/1‑148 (MQ=255)
                                                     ctTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGc                                                                               <  4:214826/129‑1 (MQ=255)
                                                     ctTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGc                                                                               >  3:214826/1‑129 (MQ=255)
                                                     ctTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTcg                                                            >  4:734684/1‑148 (MQ=255)
                                                       tttAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCg                                                                                        <  1:553466/118‑1 (MQ=255)
                                                       tttAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCg                                                                                        >  2:553466/1‑118 (MQ=255)
                                                            tCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGt                                                                                                     >  2:195160/1‑100 (MQ=255)
                                                            tCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGt                                                                                                     <  1:195160/100‑1 (MQ=255)
                                                               gCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGAtt                                                  >  3:786644/1‑148 (MQ=255)
                                                                cAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTg                                                                             >  3:752329/1‑120 (MQ=255)
                                                                cAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTg                                                                             <  4:752329/120‑1 (MQ=255)
                                                                 aaCGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGa                                                <  3:734684/148‑1 (MQ=255)
                                                                  aCGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAg                                                                                                 <  2:1086587/98‑1 (MQ=255)
                                                                  aCGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAg                                                                                                 <  4:841595/98‑1 (MQ=255)
                                                                  aCGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAg                                                                                                 >  1:1086587/1‑98 (MQ=255)
                                                                  aCGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAg                                                                                                 >  3:841595/1‑98 (MQ=255)
                                                                       gTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAg                                          <  4:786644/148‑1 (MQ=255)
                                                                         gagaCGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGttt                                                                                                   >  4:36289/1‑89 (MQ=255)
                                                                         gagaCGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGttt                                                                                                   <  3:36289/89‑1 (MQ=255)
                                                                          agaCGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGa                                                         <  2:729217/130‑1 (MQ=255)
                                                                          agaCGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGa                                                         >  1:729217/1‑130 (MQ=255)
                                                                              gCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAg                                                                                            <  4:152687/91‑1 (MQ=255)
                                                                              gCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAg                                                                                            >  3:152687/1‑91 (MQ=255)
                                                                                   cAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATaa                                                                                    >  4:829938/1‑94 (MQ=255)
                                                                                   cAATCGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATaa                                                                                    <  3:829938/94‑1 (MQ=255)
                                                                                     aTTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCCTTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGa                                                         <  4:743923/119‑1 (MQ=255)
                                                                                     aTTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGa                                                         >  3:743923/1‑119 (MQ=255)
                                                                                     aTTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTg                            >  3:231825/1‑148 (MQ=255)
                                                                                          tgGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGAt                                                   >  1:100348/1‑120 (MQ=255)
                                                                                          tgGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGAt                                                   <  2:100348/120‑1 (MQ=255)
                                                                                          tgGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACaa                       <  1:774553/148‑1 (MQ=255)
                                                                                             cGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGc                                                                               >  2:457798/1‑89 (MQ=255)
                                                                                             cGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGc                                                                               <  1:457798/89‑1 (MQ=255)
                                                                                                aTCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGt                                                                                                     <  3:978086/64‑1 (MQ=255)
                                                                                                aTCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGt                                                                                                     >  4:978086/1‑64 (MQ=255)
                                                                                                    ggCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATaaa                                                                                   >  4:510284/1‑78 (MQ=255)
                                                                                                    ggCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATaaa                                                                                   <  3:510284/78‑1 (MQ=255)
                                                                                                    ggCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGc                                         >  2:231296/1‑120 (MQ=255)
                                                                                                    ggCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGc                                         <  1:231296/120‑1 (MQ=255)
                                                                                                      cTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGa                                                <  3:573519/111‑1 (MQ=255)
                                                                                                      cTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGa                                                >  4:573519/1‑111 (MQ=255)
                                                                                                       tCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAg                                                                                            <  2:267940/66‑1 (MQ=255)
                                                                                                       tCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAg                                                                                            >  1:267940/1‑66 (MQ=255)
                                                                                                        ccGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGttt                                                                                                   <  2:352762/58‑1 (MQ=255)
                                                                                                        ccGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGttt                                                                                                   >  1:352762/1‑58 (MQ=255)
                                                                                                            tAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGc                                                                               <  4:1128411/74‑1 (MQ=255)
                                                                                                            tAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGc                                                                               >  3:1128411/1‑74 (MQ=255)
                                                                                                              aGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGc   <  4:231825/148‑1 (MQ=255)
                                                                                                                 cAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAg                                                                                            <  1:743372/56‑1 (MQ=255)
                                                                                                                 cAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAg                                                                                            >  2:743372/1‑56 (MQ=255)
                                                                                                                       cTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCg                                        >  3:165721/1‑102 (MQ=255)
                                                                                                                       cTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCg                                        <  4:165721/102‑1 (MQ=255)
                                                                                                                                  ttGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGa                                                >  3:907265/1‑83 (MQ=255)
                                                                                                                                  ttGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGa                                                <  4:907265/83‑1 (MQ=255)
                                                                                                                                        cgcCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGa                                                <  3:1004380/77‑1 (MQ=255)
                                                                                                                                        cgcCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGa                                                >  4:1004380/1‑77 (MQ=255)
                                                                                                                                        cgcCTGTTCCAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGa                                  <  3:127750/91‑1 (MQ=255)
                                                                                                                                        cgcCTGTTCCAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGa                                  >  2:50853/1‑91 (MQ=255)
                                                                                                                                        cgcCTGTTCCAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGa                                  >  4:127750/1‑91 (MQ=255)
                                                                                                                                        cgcCTGTTCCAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGa                                  <  1:50853/91‑1 (MQ=255)
                                                                                                                                          ccTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCg  >  3:150690/1‑121 (MQ=255)
                                                                                                                                          ccTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCg  <  4:150690/121‑1 (MQ=255)
                                                                                                                                          ccTGTTCCAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGtt                                                                           <  1:689243/48‑1 (MQ=38)
                                                                                                                                          ccTGTTCCAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGtt                                                                           >  2:689243/1‑48 (MQ=38)
                                                                                                                                            tGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTa                                           <  2:1085580/78‑1 (MQ=255)
                                                                                                                                            tGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTa                                           >  1:1085580/1‑78 (MQ=255)
                                                                                                                                               tGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGca                          >  2:1033705/1‑92 (MQ=255)
                                                                                                                                               tGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGca                          <  1:1033705/92‑1 (MQ=255)
                                                                                                                                                 gAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATgcg                  >  1:124792/1‑98 (MQ=255)
                                                                                                                                                 gAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATgcg                  <  2:124792/98‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                 
GAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCG  >  NC_000913/867573‑867831

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