Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,793,667 T→G 9.1% E148A (GAA→GCA) ‡ btuE ← glutathione peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,793,6670TG 9.1% 244.2 / 13.7 88E148A (GAA→GCA) ‡btuEglutathione peroxidase
Reads supporting (aligned to +/- strand):  ref base T (39/41);  new base G (4/4);  total (43/45)
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.97e-01

GTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACT  >  NC_000913/1793518‑1793816
                                                                                                                                                     |                                                                                                                                                     
gTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAAttttt                                                                                                                                                       <  1:1302406/150‑1 (MQ=255)
 tAGATTCCGCTACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCGAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGTGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTc                                                                                                                                                      <  1:1715862/150‑1 (MQ=255)
 tAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTc                                                                                                                                                      <  1:357430/150‑1 (MQ=255)
   taTTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCaa                                                                                                                                                    <  2:1233895/149‑1 (MQ=255)
     ttCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGATAAACGCTGGATGACTTTTACGTCCATGCCAACCAGGAATTTTTCaaaa                                                                                                                                                  >  1:304981/1‑150 (MQ=255)
       ccGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCCGGAATTTTTCAAAAtt                                                                                                                                                >  2:1227792/1‑150 (MQ=255)
        cGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTc                                                                                                                                               >  2:786728/1‑150 (MQ=255)
         gCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTcc                                                                                                                                              >  2:1384928/1‑150 (MQ=255)
            aCATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATa                                                                                                                                           <  1:1224909/150‑1 (MQ=255)
            aCATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTGAAAAATTcctgg                                                                                                                                           <  1:1346776/150‑4 (MQ=255)
             cATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACTATGGGATCTTCCGGCGTCATATCCGGGGAAACACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTACATaa                                                                                                                                          >  2:1435388/1‑150 (MQ=255)
             cATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTGAAAAATTcctggt                                                                                                                                          <  1:551445/150‑5 (MQ=255)
             cATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTGAAAAATTcctggt                                                                                                                                          <  1:1438741/150‑5 (MQ=255)
             cATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTATCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTGAAAAATTcctggt                                                                                                                                          <  1:1040787/150‑5 (MQ=255)
               tCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATaaaa                                                                                                                                        >  1:1169233/1‑150 (MQ=255)
                   gTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATatc                                                                                                                                    >  2:268582/1‑150 (MQ=255)
                     aaCTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATatcat                                                                                                                                  <  2:448861/150‑1 (MQ=255)
                          cATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCggg                                                                                                                             >  1:1244344/1‑150 (MQ=255)
                             cacaATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTAc                                                                                                                          <  2:1090014/150‑1 (MQ=255)
                              acaATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTGCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACa                                                                                                                         >  2:1219603/1‑150 (MQ=255)
                                        ttattTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACg                                                                                                               >  2:455101/1‑150 (MQ=255)
                                          attTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCGCTGCCAACCAGGAATTTTTCAAAATTCCATTAAATATCATCCGGGTACAGCGGTGCACGGc                                                                                                             >  2:1032344/1‑150 (MQ=255)
                                            tttGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCt                                                                                                           <  1:372283/150‑1 (MQ=255)
                                              tGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTATCCGTCACTGCCAACCAGGAATATTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCttt                                                                                                         >  1:510247/1‑150 (MQ=255)
                                                ccAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATGTCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGc                                                                                                       >  2:525245/1‑150 (MQ=255)
                                                ccAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGc                                                                                                       <  2:123014/150‑1 (MQ=255)
                                                 cAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCt                                                                                                      <  2:506806/150‑1 (MQ=255)
                                                 cAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCt                                                                                                      <  1:1539440/150‑1 (MQ=255)
                                                   aCGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGa                                                                                                    <  1:85838/150‑1 (MQ=255)
                                                    cGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAATCGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTATGCTGAc                                                                                                   >  1:1666154/1‑150 (MQ=255)
                                                    cGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGAc                                                                                                   >  1:133802/1‑150 (MQ=255)
                                                      ccAGTTTAATGCTTTCCATCACAATGGGAGCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTACCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCa                                                                                                 <  2:1346170/150‑1 (MQ=255)
                                                       cAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACAAGGAATTTTTCAAAATTCCATAAAATATCATCAGGGTACAGCGGTGCACGGCCTTTGCTGACCAt                                                                                                >  2:1360952/1‑150 (MQ=255)
                                                         gTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATAc                                                                                              >  2:528522/1‑150 (MQ=255)
                                                           ttAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACgg                                                                                            <  1:628566/150‑1 (MQ=255)
                                                              aTGCTTTACATCACAATGGGATCTTCCGGCGTAATATGAGGGGAAAAACGCTGGAGGTCTTTTTCGTCCCTGCCTACCAGGAATTTTTCATAATTCCTTAAAATATCTTCCGGGTAGAGCGGTGCACGGCCTTTGCTGTCCATACTGGCa                                                                                         <  1:905010/150‑1 (MQ=255)
                                                               tGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGTAATTTTTCAAAAATCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCAt                                                                                        <  2:885924/150‑1 (MQ=255)
                                                                 cTTTCCATCACAATGGGATCTTCCGGCGTCATGTCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAg                                                                                      <  1:525245/150‑1 (MQ=255)
                                                                 cTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAg                                                                                      <  1:1555949/150‑1 (MQ=255)
                                                                      cATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATcc                                                                                 <  1:889669/150‑1 (MQ=255)
                                                                      cATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATcc                                                                                 <  2:1480018/150‑1 (MQ=255)
                                                                       aTCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCg                                                                                <  2:806795/150‑1 (MQ=255)
                                                                         cacaATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTGCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGct                                                                              <  1:1219603/150‑1 (MQ=255)
                                                                          acaATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGctc                                                                             >  1:796115/1‑150 (MQ=255)
                                                                           caATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGctct                                                                            <  1:212509/150‑1 (MQ=255)
                                                                                  aTCTTCTGGGGTTATATCCCGGGAAAAACGCTGGATGAATTTTACGTCCCTGCAAACATGGAATTATTCAAAATTCCATAAAATATTATCAGTGTACAGGGGTGAACGGCCTTTGCTGACCATACTGGCATAGAATCCGCTCTCTTACgg                                                                     <  2:1089248/150‑1 (MQ=255)
                                                                                   tCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCACTCTCTTCCGgc                                                                    <  2:1037587/150‑1 (MQ=255)
                                                                                    cTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGgcg                                                                   >  2:930522/1‑150 (MQ=255)
                                                                                    cTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGgcg                                                                   <  2:926747/150‑1 (MQ=255)
                                                                                     ttCCGGCGTCATATCCGGGGAAAAACGCGGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGgcgc                                                                  <  1:1354653/150‑1 (MQ=255)
                                                                                       ccGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGcgcgc                                                                >  2:1344431/1‑149 (MQ=255)
                                                                                        cGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGAc                                                               >  1:1191415/1‑150 (MQ=255)
                                                                                                tatCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCg                                                       >  2:933804/1‑150 (MQ=255)
                                                                                                     ggggAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCt                                                  <  1:1061931/150‑1 (MQ=255)
                                                                                                     ggggAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCt                                                  >  2:1050704/1‑150 (MQ=255)
                                                                                                        gAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCg                                               >  1:1663067/1‑150 (MQ=255)
                                                                                                             aCGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGcaat                                          >  1:1129650/1‑150 (MQ=255)
                                                                                                                cTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGcaatcaa                                       <  1:930522/150‑1 (MQ=255)
                                                                                                                 gggATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGcaatcaat                                      <  1:436917/149‑1 (MQ=255)
                                                                                                                  ggATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAAtt                                     >  1:1649338/1‑150 (MQ=255)
                                                                                                                     tGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCACACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAAttttt                                  <  2:331184/150‑1 (MQ=255)
                                                                                                                      gACTTTTCCGTCCCTGCCACCCAGGAATTTATCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCTACTGCGGTCGGCGCTGCGGCAATCAATTTTTg                                 <  1:299533/150‑1 (MQ=255)
                                                                                                                      gACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGGGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTg                                 <  1:716692/150‑1 (MQ=255)
                                                                                                                      gACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTg                                 <  2:796115/150‑1 (MQ=255)
                                                                                                                       aCTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGa                                >  1:877504/1‑150 (MQ=255)
                                                                                                                        cTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCGTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGCCTGCGGTCGGCGCTGCGGCAATCAATTTTTGAt                               >  1:63610/1‑150 (MQ=255)
                                                                                                                          tttCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATAc                             <  1:1347238/150‑1 (MQ=255)
                                                                                                                          tttCCGTCCCTCCCAACCAGGAATTTTTCAAAATTCCAAAAAATATCATCCGGGTACAGCGCTGCAAGGCCTATGCTGACAAAACGTGAATAGAATCCGCTCTCTTACGGCGCGTCAGCGGTCGGCGCTGCGGCAATCAATATTTGATAc                             >  1:269407/1‑150 (MQ=255)
                                                                                                                              cGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCAGTCGGCGCTGCGGCAATCAATTTTTGATACAGCg                         >  1:1533183/1‑150 (MQ=255)
                                                                                                                               gTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCgg                        >  2:1233780/1‑150 (MQ=255)
                                                                                                                                tCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGAGGCAATCAATTTTTGATACAGCGGa                       >  1:76666/1‑150 (MQ=255)
                                                                                                                                  ccTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATg                     <  2:1663067/150‑1 (MQ=255)
                                                                                                                                   cTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATgg                    <  2:192451/150‑1 (MQ=255)
                                                                                                                                    tGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGGACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGc                   <  2:1243651/150‑1 (MQ=255)
                                                                                                                                      ccAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGt                 <  2:1169233/150‑1 (MQ=255)
                                                                                                                                       cAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTc                >  1:1297465/1‑150 (MQ=255)
                                                                                                                                        aaCCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTcc               >  2:1010471/1‑150 (MQ=255)
                                                                                                                                        aaCCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTcc               <  1:937981/150‑1 (MQ=255)
                                                                                                                                         ttatggAATTTTGAAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACTGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCTATCAATTTTTGATACAGCGGCTGGCGTCCt              >  1:973751/5‑150 (MQ=255)
                                                                                                                                         ttatggAATTTTGAAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGAGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTccg              >  2:684493/5‑149 (MQ=255)
                                                                                                                                         aCCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCt              >  1:1458026/1‑150 (MQ=255)
                                                                                                                                           cAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGAACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTc            >  1:1041806/1‑150 (MQ=255)
                                                                                                                                           atggAATTTTGAAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTc            >  1:409851/3‑150 (MQ=255)
                                                                                                                                            tggAATTTTGAAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCg           >  1:1612909/2‑150 (MQ=255)
                                                                                                                                             ggAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGc                                                <  1:976119/112‑1 (MQ=255)
                                                                                                                                             ggAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGc                                                >  2:976119/1‑112 (MQ=255)
                                                                                                                                              gAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGcc         >  1:853477/1‑150 (MQ=255)
                                                                                                                                               aaTTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCa        <  2:1324790/150‑1 (MQ=255)
                                                                                                                                                aTTTTTCAAAATTCCAGAATATATCATCCGGGTACAGCGGTGGACGGCCTTTGCTGACCATACGGGCAGAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCGGCGGCAATCAATTTTTGATACAGCGGATGGCGTCGTTCGCCAt       >  1:246720/1‑150 (MQ=255)
                                                                                                                                                   tttCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGGCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTaa    <  2:1363952/150‑1 (MQ=255)
                                                                                                                                                    ttCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAAc   >  2:800909/1‑150 (MQ=255)
                                                                                                                                                     tCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACt  >  1:806461/1‑150 (MQ=255)
                                                                                                                                                     tCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACt  >  1:130046/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACT  >  NC_000913/1793518‑1793816

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