Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,715,020 G→T 8.9% intergenic (+38/+6) gstA → / ← pdxY glutathione S‑transferase GstA/pyridoxal kinase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,715,0200GT 8.9% 217.3 / 9.0 67intergenic (+38/+6)gstA/pdxYglutathione S‑transferase GstA/pyridoxal kinase 2
Reads supporting (aligned to +/- strand):  ref base G (29/32);  new base T (2/4);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

ATACGCGGTGAAACTGAATCTGGAAGGGTTAGAGCACATTGCAGCATTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGAGCAGCCACCACTTGCAGCTCATATTCCTGCATTGCTTTGGTGGTCACCATGATTTCGTAGACTGCAGCGGTCACATGTTCCAGCGCCTCCTG  >  NC_000913/1714871‑1715169
                                                                                                                                                     |                                                                                                                                                     
aTACGCGGTGAAACTGAATCTGGAAGGGTTAGAGCACATTGCAGCATTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATgag                                                                                                                                                       <  2:429946/150‑1 (MQ=255)
aTACGCGGTGAAACTGAATCTGGAAGGGTTAGAGCACATTGCAGCATTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATgag                                                                                                                                                       <  2:424515/150‑1 (MQ=255)
 tACGCGGTGAAACTGAATCTGGAAGGGTTAGAGCACATTGCAGCATTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATgaga                                                                                                                                                      >  1:486617/1‑150 (MQ=255)
 tACGCGGTGAAACTGAATCTGGAAGGGTTAGAGCACATTGCAGCATTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATgaga                                                                                                                                                      <  1:1402431/150‑1 (MQ=255)
     cgGTGAAACTGAATCTGGAAGGGTTAGAGCACATTGCAGCATTTATGCTACGTATGGCTGAACGTCCGGAAGTACAAGTCGCGCTGTCAGCGGAAGGCTTAAATTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGTAATGAGATAtt                                                                                                                                                  <  2:461827/150‑1 (MQ=255)
        tGAAACTGAATCTGGAAGGGTTAGAGCACATTGCAGCATTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCa                                                                                                                                               <  1:104552/150‑1 (MQ=255)
         gAAACTGAATCTGGAAGGGTTAGAGCACATTGCAGCATTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCag                                                                                                                                              <  1:1190646/150‑1 (MQ=255)
         gAAACTGAATCTGGAAGGGTTAGAGCACATTGCAGCATTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCag                                                                                                                                              >  1:1205348/1‑150 (MQ=255)
           aaCTGAATCTGGAAGGGTTAGAGCACATTGCAGCATTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCagag                                                                                                                                            >  2:117702/1‑150 (MQ=255)
             cTGAATCTGGAAGGGTTAGAGCACATTGCAGCATTTATGCAACGTATGGCTGCACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCt                                                                                                                                          >  2:1521634/1‑150 (MQ=255)
             cTGAATCTGGAAGGGTTAGAGCACATTGCAGCATTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCt                                                                                                                                          >  2:1537906/1‑150 (MQ=255)
             cTGAATCTGGAAGGGTTAGAGCACATTGCAGCATTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGAGCTGTAAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTTTGGGACGAAATGAGATATTTCAGAGCt                                                                                                                                          >  1:540254/1‑150 (MQ=255)
              tGAATCTGGAAGGGTTAGAGCACATTGCAGCATTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCtt                                                                                                                                         >  1:472390/1‑150 (MQ=255)
                 aTCTGGAAGGGTTAGAGCACATTGCAGCATTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTttgt                                                                                                                                      <  1:541306/150‑1 (MQ=255)
                   cTGGAAGGGTTAGAGCACATTGCAGCATTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTttgttg                                                                                                                                    <  2:467330/150‑1 (MQ=255)
                   cTGGAAGGGTTAGAGCACATTGCAGCATTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTttgttg                                                                                                                                    <  2:1206463/150‑1 (MQ=255)
                    tagaatcgTTAGAGCACATTGAAGAATTTATGCAACGTATCGCTGAACGTCCGGAAGTACAAAACACGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTgc                                                                                                                                   >  2:1058345/8‑150 (MQ=255)
                       aaGGGTTAGAGCACATTGCAGCATTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCt                                                                                                                                <  2:486617/150‑1 (MQ=255)
                          ggTTAGAGCACATTGCAGCATTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGaa                                                                                                                             >  1:620532/1‑150 (MQ=255)
                           gTTAGAGCACATTGCAGCATTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAg                                                                                                                            <  2:620532/150‑1 (MQ=255)
                              agcgCACATTGCAGCATTTATGAAACGTATGCCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTAGTTGTGGGCCGAAATGAGATATTTCACAGCTTTGTTGCGCTGAATTaa                                                                                                                         >  2:1363772/4‑150 (MQ=255)
                                 gCACATTGCAGCATTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGt                                                                                                                      >  2:158827/1‑150 (MQ=255)
                                 gCACATTGCAGCATTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGt                                                                                                                      >  1:638787/1‑150 (MQ=255)
                                        gcagcaTTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTttgt                                                                                                                                      <  2:986789/127‑1 (MQ=255)
                                        gcagcaTTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTttgt                                                                                                                                      >  1:986789/1‑127 (MQ=255)
                                        gcagcaTTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGtt                                                                                                               >  1:244510/1‑150 (MQ=255)
                                        gcagcaTTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGtt                                                                                                               <  2:1447785/150‑1 (MQ=255)
                                        gcagcaTTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGtt                                                                                                               <  2:74334/150‑1 (MQ=255)
                                        gcagcaTTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGtt                                                                                                               >  1:556131/1‑150 (MQ=255)
                                         cagcaTTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGttt                                                                                                              <  2:537484/150‑1 (MQ=255)
                                          agcaTTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTg                                                                                                             <  1:876241/150‑1 (MQ=255)
                                           gcaTTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGATATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTgg                                                                                                            <  1:655933/150‑1 (MQ=255)
                                           gcaTTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTgg                                                                                                            >  1:359490/1‑150 (MQ=255)
                                            caTTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTATTAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGc                                                                                                           <  2:614852/150‑1 (MQ=255)
                                            caTTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGATATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGc                                                                                                           <  1:485862/150‑1 (MQ=255)
                                            caTTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGc                                                                                                           <  1:1511331/150‑1 (MQ=255)
                                              tttATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGATATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCaa                                                                                                         <  1:1501929/150‑1 (MQ=255)
                                              tttATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCaa                                                                                                         >  1:1613549/1‑150 (MQ=255)
                                                 aTGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGATATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATAc                                                                                                      >  2:200824/1‑150 (MQ=255)
                                                 aTGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATAc                                                                                                      <  2:366643/150‑1 (MQ=255)
                                                            ggCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGagca                                                                                           >  1:1532993/1‑150 (MQ=255)
                                                                  aCGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGAGCAGccacc                                                                                     <  2:621876/150‑1 (MQ=255)
                                                                  aCGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGAGCAGccacc                                                                                     <  1:1048833/150‑1 (MQ=255)
                                                                        ggAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGAGCAGCCACCACTTGc                                                                               >  1:1047608/1‑150 (MQ=255)
                                                                        ggAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGAGCAGCCACCACTTGc                                                                               >  1:1067283/1‑150 (MQ=255)
                                                                              aCAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGa                                                                                              <  2:461889/129‑1 (MQ=255)
                                                                              aCAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGa                                                                                              >  1:461889/1‑129 (MQ=255)
                                                                                   aCGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGTGCCGAAATGATATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGAGCAGCCACCACTTGCAGCTCATATTc                                                                    <  1:1209089/150‑1 (MQ=255)
                                                                                        cTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTAGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTTCGCTTAAGTAATGTTCTGGTTTGGCAATACGATCCTGAGCAGCCACCACTTGCAGCTCATATTCCTGCa                                                               <  1:1021463/150‑1 (MQ=255)
                                                                                                       ttAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGAGCAGCCCCCACTTGCCGCTCATATTCCTGCATTGCTTTGGTGGTcc                                                >  1:81151/1‑149 (MQ=255)
                                                                                                         aaaGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGAGCAGCCACCACTTGCAGCTCATATTCCTGCATTGCTTTGGTGGTCAcc                                              >  1:605057/1‑150 (MQ=255)
                                                                                                             tAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGATATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGAGCAGCCACCACTTGCAGCTCATATTCCTGCATTGCTTTGGTGGTCACCATGa                                          >  1:989757/1‑150 (MQ=255)
                                                                                                                  gCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGAGCAGCCACCACTTGCAGCTCATATTCCTGCATTGCTTTGGTGGTCACCATGATTTCg                                     <  1:158827/150‑1 (MQ=255)
                                                                                                                  gCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCAGAGCAGCCCCCACTTGCAGCTCATATTCCTGCATTGCTTTGGTGGTCCCCATGATTTCg                                     >  2:67186/1‑150 (MQ=255)
                                                                                                                  gCCATCGGCCCACAGGGCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGAGCAGCCACCACTTGCAGCTCATATTCCTGCATTGCTTTGGTGGTCACCATGATTTCg                                     <  2:1336772/150‑1 (MQ=255)
                                                                                                                       cGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGAGCAGCCACCACTTGCAGCTCATATTCCTGCATTGCTTTGGTGGTCACCATGATTTCGTAGAc                                <  1:409855/150‑1 (MQ=255)
                                                                                                                       cGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGAGCAGCCACCACTTGCAGCTCATATTCCTGCATTGCTTTGGTGGTCACCATGATTTCGTAGAc                                <  1:1229654/150‑1 (MQ=255)
                                                                                                                         gCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGAGCAGCCACCACTTGCAGCTCATATTCCTGCATTGCTTTGGTGGTCACCATGATTTCGTAGACTg                              <  2:472390/150‑1 (MQ=255)
                                                                                                                          cccACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGAGCAGCCACCACTTGCAGCTCATATTCCTGCATTGCTTTGGTGGTCACCATGATTTCGTAGACTGc                             <  1:940874/150‑1 (MQ=255)
                                                                                                                           ccACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGAGCAGCCACCACTTGCAGCTCATATTCCTGCATTGCTTTGGTGGTCACCATGATTTCGTAGACTGCa                            >  1:791408/1‑150 (MQ=255)
                                                                                                                               aGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGAGCAGCCACCACTTGCAGCTCATATTCCTGCATTGCTTTGGTGGTCACCATGATTTCGTAGACTGCAGCgg                        <  2:1067283/150‑1 (MQ=255)
                                                                                                                                 gTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGAGCAGCCACCACTTGCAGCTCATATTCCTGCATTGCTTTGGTGGTCACCATGATTTCGTAGACTGCAGCGGTc                      <  1:1318676/150‑1 (MQ=255)
                                                                                                                                   cGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGAGCAGCCACCACTTGCAGCTCATATTCCTGCATTGCTTTGGTGGTCACCATGATTTCGTAGACTGCAGCGGTcac                    >  1:1235584/1‑150 (MQ=255)
                                                                                                                                    cctgtgGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGAGCAGCCACCACTTGCAGCTCATATTCCTGCATTGCTTTGGTGGTCACCATGATTTCGTAGACTGCAGCGGTcaca                   <  1:1330208/148‑1 (MQ=255)
                                                                                                                                       gtgGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGAGCAGCCACCACTTGCAGCTCATATTCCTGCATTGCTTtgg                                                      >  2:581942/1‑112 (MQ=255)
                                                                                                                                       gtgGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGAGCAGCCACCACTTGCAGCTCATATTCCTGCATTGCTTtgg                                                      <  1:581942/112‑1 (MQ=255)
                                                                                                                                        tgGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGAGCAGCCACCACTTGCAGCTCATATTCCTGCATTGCTTTGGTGGTCACCATGATTTCGTAGACTGCAGCGGTCACATGtt               >  1:198507/1‑150 (MQ=255)
                                                                                                                                         gggCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGAGCAGCCACCACTTGCAGCTCATATTCCTGCATTGCTTTGGTGGTCACCATGATTTCGTAGACTGCAGCGGTCACATGTTc              <  1:1017407/150‑1 (MQ=255)
                                                                                                                                         gggCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGAGCAGCCACCACTTGCAGCTCATATTCCTGCATTGCTTTGGTGGTCACCATGATTTCGTAGACTGCAGCGGTCACATGTTc              >  1:296844/1‑150 (MQ=255)
                                                                                                                                            ccGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGAGCAGCCACCACTTGCAGCTCATATTCCTGCATTGCTTTGGTGGTCACCATGATTTCGTAGACTGCAGCGGTCACATGTTCCAg           >  1:421616/1‑150 (MQ=255)
                                                                                                                                                aaTGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGAGCAGCCACCACTTGCAGCTCATATTCCTGCATTGCTTTGGTGGTCACCATGATTTCGTAGACTGCAGCGGTCACATGTTCCAGCGcc       <  2:1374915/150‑1 (MQ=255)
                                                                                                                                                     gaTATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGAGCAGCCACCACTTGCAGCTCATATTCCTGCATTGCTTTGGTGGTCACCATGATTTCGTAGACTGCAGCGGTCACATGTTCCAGCGCCTCCTg  >  2:446718/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATACGCGGTGAAACTGAATCTGGAAGGGTTAGAGCACATTGCAGCATTTATGCAACGTATGGCTGAACGTCCGGAAGTACAAGACGCGCTGTCAGCGGAAGGCTTAAAGTAATAGCCATCGGCCCACAGGTCGTTGTGGGCCGAAATGAGATATTTCAGAGCTTTGTTGCGCTGAAGTAATGTTCTGGTTTGGCAATACGATCCTGAGCAGCCACCACTTGCAGCTCATATTCCTGCATTGCTTTGGTGGTCACCATGATTTCGTAGACTGCAGCGGTCACATGTTCCAGCGCCTCCTG  >  NC_000913/1714871‑1715169

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