Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,793,668 C→A 8.9% E148* (GAA→TAA) ‡ btuE ← glutathione peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,793,6680CA 8.9% 248.3 / 13.7 90E148* (GAA→TAA) ‡btuEglutathione peroxidase
Reads supporting (aligned to +/- strand):  ref base C (42/40);  new base A (4/4);  total (46/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 = 1.83e-01

TAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACT  >  NC_000913/1793519‑1793816
                                                                                                                                                     |                                                                                                                                                    
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:1438741/150‑5 (MQ=255)
            cATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTGAAAAATTcctggt                                                                                                                                          <  1:551445/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)
                                                                                                                                                     |                                                                                                                                                    
TAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACT  >  NC_000913/1793519‑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: ‑