Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,506,790 C→T 5.2% A369T (GCT→ACT)  fryC ← putative enzyme IIC component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,506,7900CT 5.2% 269.2 / 3.2 77A369T (GCT→ACT) fryCputative enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base C (38/35);  new base T (2/2);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CTTCTTACAGGCTATCGATTAACAATTTGCCTTTACGGAACATCATCAGACGCAGGAACACCACCATCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAGCGCAAACGGAATCGCACCCTCACTGATCCCCATAAACGCCAGGAACA  >  NC_000913/2506643‑2506939
                                                                                                                                                   |                                                                                                                                                     
cttcttACAGGCTATCGATTAACAATTTGCCTTTACGGAACATCATCAGACGCAGGAACACCACCATCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCgg                                                                                                                                                     <  2:1231472/150‑1 (MQ=255)
 ttcttACAGGCTATCGATTAACAATTTGCCTTTACGGAACATCATCAGACGCAGGAACACCACCATCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGa                                                                                                                                                    >  1:1104378/1‑150 (MQ=255)
 ttcttACAGGCTATCGATTAACAATTTGCCTTTACGGAACATCATCAGACGCAGGAACACCACCATCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGa                                                                                                                                                    >  1:1593635/1‑150 (MQ=255)
    ttACAGGCTATCGATTAACAATTTGCCTTTACGGAACATCATCAGACGCAGGAACACCACCATCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTc                                                                                                                                                 <  1:212409/150‑1 (MQ=255)
    ttACAGGCTATCGATTAACAATTTGCCTTTACGGAACATCATCAGACGCAGGAACACCACCATCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTc                                                                                                                                                 <  1:1083129/150‑1 (MQ=255)
        aGGCTATCGATTAACAATTTGCCTTTACGGAACATCATCAGACGCAGGAACACCACCATCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCggg                                                                                                                                             <  1:1676885/150‑1 (MQ=255)
                       aaTTTGCCTTTACGGAACATCATCAGACGCAGGAACACCACCATCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTgcg                                                                                                                              <  1:1549047/150‑1 (MQ=255)
                       aaTTTGCCTTTACGGAACATCATCAGACGCAGGAACACCACCATCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTgcg                                                                                                                              <  1:397123/150‑1 (MQ=255)
                           tGCCTTTACGGAACATCATCAGACGCAGGAACACCACCATCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCcca                                                                                                                          <  1:1036302/150‑1 (MQ=255)
                           tGCCTTTACGGAACATCATCAGACGCAGGAACACCACCATCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCcca                                                                                                                          <  2:171132/150‑1 (MQ=255)
                           tGCCTTTACGGAACATCATCAGACGCAGGAACACCACCATCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCcca                                                                                                                          <  1:585280/150‑1 (MQ=255)
                           tGCCTTTACGGAACATCATCAGACGCAGGAACACCACCATCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCcca                                                                                                                          <  2:690781/150‑1 (MQ=255)
                              cTTTACGGAACATCATCAGACGCAGGAACACCACCATCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCccagcc                                                                                                                       <  1:136541/150‑1 (MQ=255)
                               tttACGGAACATCATCAGACGCAGGAACACCACCATCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCccagcca                                                                                                                      <  1:520313/150‑1 (MQ=255)
                                 tACGGAACATCATCAGACGCAGGAACACCACCATCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGa                                                                                                                    >  1:38398/1‑150 (MQ=255)
                                         atcatcAGACGCAGGAACACCACCATCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAgcggcg                                                                                                            >  2:1713061/1‑150 (MQ=255)
                                             tcaGACGCAGGAACACCACCATCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTg                                                                                                        >  1:1065469/1‑150 (MQ=255)
                                                gACGCAGGATCACCACCATCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGc                                                                                                     <  1:274896/150‑1 (MQ=255)
                                                gACGCAGGAACACCACCATCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGc                                                                                                     <  2:488101/150‑1 (MQ=255)
                                                     aGGAACACCACCATCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACa                                                                                                >  1:1125380/1‑150 (MQ=255)
                                                     aGGAACACCACCATCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACa                                                                                                >  2:1517914/1‑150 (MQ=255)
                                                        aaCACCACCATCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATa                                                                                             <  1:471687/150‑1 (MQ=255)
                                                        aaCACCACCATCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATa                                                                                             <  1:471592/150‑1 (MQ=255)
                                                         acacCACCATCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAg                                                                                            <  2:628706/150‑1 (MQ=255)
                                                             caccaTCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGcc                                                                                        <  2:1260905/150‑1 (MQ=255)
                                                             caccaTCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGcc                                                                                        >  2:265203/1‑150 (MQ=255)
                                                              accaTCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCg                                                                                       >  2:357053/1‑150 (MQ=255)
                                                              accaTCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCg                                                                                       >  2:357054/1‑150 (MQ=255)
                                                              accaTCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCg                                                                                       >  2:924941/1‑150 (MQ=255)
                                                                 aTCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGAcc                                                                                    >  1:1228366/1‑150 (MQ=255)
                                                                 aTCAGAGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGAcc                                                                                    >  1:1228382/1‑150 (MQ=255)
                                                                    aGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCata                                                                                 >  1:755871/1‑150 (MQ=255)
                                                                    aGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCata                                                                                 >  1:433088/1‑150 (MQ=255)
                                                                     gcgcAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCatat                                                                                >  1:839930/1‑150 (MQ=255)
                                                                       gcAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATac                                                                              <  1:746237/150‑1 (MQ=255)
                                                                           ttATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATacgacg                                                                          >  1:1602891/1‑150 (MQ=255)
                                                                                  tGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGc                                                                   >  1:1545928/1‑150 (MQ=255)
                                                                                  tGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGc                                                                   >  1:1676617/1‑150 (MQ=255)
                                                                                     ccccAGCGCGATCCCCGCCATATAGACGCCCTGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGt                                                                <  1:1550255/150‑1 (MQ=255)
                                                                                      cccAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTg                                                               >  2:404739/1‑150 (MQ=255)
                                                                                       ccAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGTGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGa                                                              >  1:1384891/1‑150 (MQ=255)
                                                                                       ccAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGTGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGa                                                              >  1:523532/1‑150 (MQ=255)
                                                                                       ccAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGa                                                              >  1:245010/1‑150 (MQ=255)
                                                                                        cAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGAt                                                             >  1:108143/1‑150 (MQ=255)
                                                                                              gATCCCCGCCATATAGACGCCCAGGTAAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTAGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCt                                                       >  1:1037833/1‑150 (MQ=255)
                                                                                               aTCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCtt                                                      >  2:783805/1‑150 (MQ=255)
                                                                                                    cGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAg                                                 <  2:1390472/150‑1 (MQ=255)
                                                                                                    cGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAg                                                 <  2:1390532/150‑1 (MQ=255)
                                                                                                     gCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAgc                                                <  2:637437/150‑1 (MQ=255)
                                                                                                     gCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAgc                                                <  2:494496/150‑1 (MQ=255)
                                                                                                     gCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAgc                                                <  2:1494667/150‑1 (MQ=255)
                                                                                                          ataGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAGCGCaaa                                           >  1:700905/1‑150 (MQ=255)
                                                                                                            cgACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCAGATACGACGGAATGGCGGTGATGGGGCTTTCCAGCGCAAACg                                         >  1:396106/2‑150 (MQ=255)
                                                                                                            aGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAGCGCAAACg                                         >  1:232218/1‑150 (MQ=255)
                                                                                                               cGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAGCGCAAACGGaa                                      <  1:803149/150‑1 (MQ=255)
                                                                                                                gCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAGCGCAAACGGAAt                                     >  2:1395196/1‑150 (MQ=255)
                                                                                                                gCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAGCGCAAACGGAAt                                     >  2:1395150/1‑150 (MQ=255)
                                                                                                                 cccAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAGCGCAAACGGAATc                                    <  2:1562095/150‑1 (MQ=255)
                                                                                                                      gTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAGCGCAAACGGAATCGCAcc                               <  1:467856/150‑1 (MQ=255)
                                                                                                                           tGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAGCGCAAACGGAATCGCACCCTCAc                          <  2:860413/150‑1 (MQ=255)
                                                                                                                            gACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAGCGCAAACGGAATCGCACCCTCACt                         >  1:758505/1‑150 (MQ=255)
                                                                                                                             aCCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAGCGCAAACGGACTCGCACCCTCACTg                        >  2:238950/1‑150 (MQ=255)
                                                                                                                             aCCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAGCGCAAACGGAATCGCACCCTCACTg                        >  1:1610824/1‑150 (MQ=255)
                                                                                                                             aCCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAGCGCAAACGGAATCGCACCCTCACTg                        >  1:269448/1‑150 (MQ=255)
                                                                                                                             aCCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAGCGCAAACGGAATCGCACCCTCACTg                        >  2:613198/1‑150 (MQ=255)
                                                                                                                               cAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAGCGCAAACGGAATCGCACCCTCACTGAt                      <  2:839930/150‑1 (MQ=255)
                                                                                                                               cAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAGCGCAAACGGAATCGCACCCTCACTGAt                      <  1:1635556/150‑1 (MQ=255)
                                                                                                                               cAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAGCGCAAACGGAATCGCACCCTCACTGAt                      <  2:1464771/150‑1 (MQ=255)
                                                                                                                                   ggCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAGCGCAAACGGAATCGCACCCTCACTGATcccc                  <  2:250361/150‑1 (MQ=255)
                                                                                                                                   ggCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAGCGCAAACGGAATCGCACCCTCACTGATcccc                  <  2:709088/150‑1 (MQ=255)
                                                                                                                                     ccATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAGCGCAAACGGAATCGCACCCTCACTGATCCCCAt                <  2:623084/150‑1 (MQ=255)
                                                                                                                                       aTGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAGCGCAAACGGAATCGCACCCTCACTGATCCCCATaa              >  2:731026/1‑150 (MQ=255)
                                                                                                                                       aTGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAGCGCAAACGGAATCGCACCCTCACTGATCCCCATaa              >  1:101540/1‑150 (MQ=255)
                                                                                                                                         gCCCAGATAGTGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAGCGCAAACGGAATCGCACCCTCACTGATCCCCATAAAc            <  2:523532/150‑1 (MQ=255)
                                                                                                                                         gCCCAGATAGTGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAGCGCAAACGGAATCGCACCCTCACTGATCCCCATAAAc            <  2:1384891/150‑1 (MQ=255)
                                                                                                                                           ccAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAGCGCAAACGGAATCGCACCCTCACTGATCCCCATAAACGc          >  2:937793/1‑150 (MQ=255)
                                                                                                                                                 aGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAGCGCAAACGGAATCGCACCCTCACTGATCCCCATAAACGCCAGGaa    >  1:411667/1‑150 (MQ=255)
                                                                                                                                                   cGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAGCGCAAACGGAATCGCACCCTCACTGATCCCCATAAACGCCAGGAACa  <  2:758505/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CTTCTTACAGGCTATCGATTAACAATTTGCCTTTACGGAACATCATCAGACGCAGGAACACCACCATCAGCGCAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGACAATAGCGCCGACCATATACGACGGAATGGCGGTGATGGGGCTTTCCAGCGCAAACGGAATCGCACCCTCACTGATCCCCATAAACGCCAGGAACA  >  NC_000913/2506643‑2506939

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