Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,603,543 T→C F157L (TTT→CTT)  b1516 → putative LACI‑type transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,603,5430TC99.2% 817.6 / ‑5.2 248F157L (TTT→CTT) b1516putative LACI‑type transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (127/119);  total (128/120)
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.98e-01

Only 100 of 272 total aligned reads displayed.

ACTACATTAATCAGGGAACGCCCGCCCAGTTAGGAGGTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGTTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGCCAAAGAGCATCCAGGCTGGGAA  >  NC_000913/1603451‑1603632
                                                                                            |                                                                                         
actacATTAATCAGGGAACGCCCGCCCAGTTAGGAGGTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGt                                                                                           >  2:4607330/1‑93 (MQ=255)
  tacATTAATCAGGGAACGCCCGCCCAGTTAGGAGGTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGctt                                                                                         <  2:5178516/93‑1 (MQ=255)
   acATTAATCAGGGAACGCCCGCCCAGTTAGGAGGTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGtttt                                                                                        >  2:6987550/1‑93 (MQ=255)
   acATTAATCAGGGAACGCCCGCCCAGTTAGGAGGTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGcttt                                                                                        >  2:8215672/1‑93 (MQ=255)
     aTTAATCAGGGAACGCCCGCCCAGTTAGGAGGTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTc                                                                                      >  1:6340042/1‑93 (MQ=255)
       tAATCAGGGAACGCCCGCCCAGTTAGGAGGTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTa                                                                                    <  1:5658352/93‑1 (MQ=255)
         aTCAGGGAACGCCCGCCCAGTTAGGAGGTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACt                                                                                  <  1:6637125/93‑1 (MQ=255)
         aTCAGGGAACGCCCGCCCAGTTAGGAGGTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGAGCTTTTCTACt                                                                                  >  2:2506451/1‑93 (MQ=255)
          tCAGGGAACGCCCGCCCAGTTAGGAGGTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTc                                                                                 <  1:1159366/93‑1 (MQ=255)
            aGGGAACGCCCGCCCAGTTAGGAGGTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCaa                                                                               >  2:6833921/1‑93 (MQ=255)
             gggAACGCCCGCCCAGTTAGGAGGTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAg                                                                              >  1:8640683/1‑93 (MQ=255)
             gggAACGCCCGCCCAGTTAGGAGGTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAg                                                                              >  1:7634948/1‑93 (MQ=255)
              ggAACGCCCGCCCAGTTAGGAGGTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGc                                                                             >  1:2354127/1‑93 (MQ=255)
              ggAACGCCCGCCCAGTTAGGAGGTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGc                                                                             <  2:1827706/93‑1 (MQ=255)
                aaCGCCCGCCCAGTTAGGAGGTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGccc                                                                           <  2:8805101/93‑1 (MQ=255)
                aaCGCCCGCCCAGTTAGGAGGTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGccc                                                                           >  2:7458234/1‑93 (MQ=255)
                 aCGCCCGCCCAGTTAGGAGGTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGcccc                                                                          >  2:3526075/1‑93 (MQ=255)
                 aCGCCCGCCCAGTTAGGAGGTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGcccc                                                                          <  2:4940625/93‑1 (MQ=255)
                 aCGCCCGCCCAGTTAGGAGGTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGcccc                                                                          <  1:1564374/93‑1 (MQ=255)
                 aCGCCCGCCCAGTTAGGAGGAATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGcctc                                                                          >  2:7749166/1‑91 (MQ=255)
                  cgcccgccCAGTTAGGAGGTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCa                                                                         <  2:3751734/93‑1 (MQ=255)
                    cccgcctAGTTAGGAGGTATGTTGGTGGATATGGCGCCGCTTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCAcc                                                                       <  2:59813/93‑1 (MQ=255)
                      cgcccAGTTAGGAGGTATGTTGGTGGATGTGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTCCTCAAGCCCCACCGt                                                                     <  1:1611436/93‑1 (MQ=255)
                         ccAGTTAGGAGGTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTAc                                                                  <  1:5151856/93‑1 (MQ=255)
                         ccAGTTAGGAGGTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTAc                                                                  <  2:1120242/93‑1 (MQ=255)
                              tAGGAGGTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGAcc                                                             >  1:4951842/1‑93 (MQ=255)
                              tAGGAGGTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGAcc                                                             <  1:3904250/93‑1 (MQ=255)
                                ggaggTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCaa                                                           >  2:321999/1‑93 (MQ=255)
                                ggaggTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCaa                                                           <  1:2649830/93‑1 (MQ=255)
                                ggaggTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCaa                                                           >  2:8393645/1‑93 (MQ=255)
                                 gaggTATGTTGGTGGATATGGCGGCGCGTCAGGTGGATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGACCCACCGTTACGGACCaaa                                                          >  2:3708146/1‑93 (MQ=255)
                                 gaggTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCaaa                                                          >  2:5621430/1‑93 (MQ=255)
                                 gaggTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCaaa                                                          >  1:5632981/1‑93 (MQ=255)
                                   ggTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAAc                                                        >  1:8539620/1‑93 (MQ=255)
                                   ggTATGTTGGTGGATAGGGAGGGGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTCCTCAAGCCCCACCGTTACGGACCAAAAc                                                        <  1:1526655/93‑1 (MQ=255)
                                    gTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAAcc                                                       >  1:3763523/1‑93 (MQ=255)
                                     tATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCa                                                      <  2:4698547/93‑1 (MQ=255)
                                       tGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGt                                                    <  2:4162312/93‑1 (MQ=255)
                                         ttGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTgg                                                  >  2:7213399/1‑93 (MQ=255)
                                           ggtggATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGt                                                <  1:7980398/93‑1 (MQ=255)
                                           ggtggATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGt                                                <  1:2424985/93‑1 (MQ=255)
                                           ggtggATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGt                                                <  1:6144945/93‑1 (MQ=255)
                                              ggATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGaa                                             >  1:7510016/1‑93 (MQ=255)
                                               gATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGaaa                                            >  1:6735562/1‑93 (MQ=255)
                                               gATAGGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGaaa                                            <  2:777288/93‑1 (MQ=255)
                                                atatGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAaag                                           <  1:6874660/93‑1 (MQ=255)
                                                  atGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCTCTTTTCTCCTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAaagaa                                         <  1:4722653/93‑1 (MQ=255)
                                                  atGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAaagaa                                         <  2:1201669/93‑1 (MQ=255)
                                                  atGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAaagaa                                         >  2:202372/1‑93 (MQ=255)
                                                      cggcgCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGaagcga                                     <  1:4413406/93‑1 (MQ=255)
                                                        gcgcGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGaagcgaaa                                   >  2:6271604/1‑93 (MQ=255)
                                                        gcgcGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGaagcgaaa                                   >  1:5466436/1‑93 (MQ=255)
                                                         cgggTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGaagcgaaag                                  <  1:1097871/90‑1 (MQ=255)
                                                         cgcgTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGaagcgaaag                                  <  2:5312621/93‑1 (MQ=255)
                                                          gcgTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAATaagcgaaagc                                 <  1:7763325/93‑1 (MQ=255)
                                                          gcgTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGaagcgaaagc                                 >  1:2784205/1‑93 (MQ=255)
                                                           cgTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGaagcgaaagcg                                >  1:3303074/1‑93 (MQ=255)
                                                            gTCAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGaagcgaaagcga                               >  2:2816083/1‑93 (MQ=255)
                                                              cAGGTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCACAACCAGTGGGTGAAAGAAgcgaaagagaaa                             >  1:3603579/1‑93 (MQ=255)
                                                                ggTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAAt                           >  2:3370799/1‑93 (MQ=255)
                                                                ggTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAAt                           <  2:6052551/93‑1 (MQ=255)
                                                                 gTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATc                          >  1:2159564/1‑93 (MQ=255)
                                                                 gTGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATc                          >  1:5097722/1‑93 (MQ=255)
                                                                  tGAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCg                         <  2:4075347/93‑1 (MQ=255)
                                                                   gAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGc                        >  1:3013197/1‑93 (MQ=255)
                                                                   gAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGc                        <  1:5966626/93‑1 (MQ=255)
                                                                   gAATAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGc                        >  1:6506140/1‑93 (MQ=255)
                                                                    aaTAAAGATAAAGCCATAGTCGTGCTTTTTTATTCATGCCCCTTTGTTACGTACCAAAACCAGTGGGTGAAAGAAGTGAATGCTAAAATCGct                       >  2:8061381/1‑92 (MQ=255)
                                                                    aaTAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGcc                       >  1:8443395/1‑93 (MQ=255)
                                                                    aaTAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGcc                       >  2:1349217/1‑93 (MQ=255)
                                                                    aaTAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGAGGGTGAAGGAAGCGAAAGCGAAAATCGcc                       >  2:7630053/1‑93 (MQ=255)
                                                                     aTAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAGACCAGTGGGTGAAAGAAGCGAAAGCGAGAATCGCCa                      >  2:2945544/1‑93 (MQ=255)
                                                                      tAAAGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGCCaa                     >  1:143524/1‑93 (MQ=255)
                                                                      aaaaGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGCCaa                     <  1:5835364/92‑1 (MQ=255)
                                                                       aaaGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGCCaaa                    <  2:1165034/93‑1 (MQ=255)
                                                                        aaGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGCCAAag                   >  1:6948760/1‑93 (MQ=255)
                                                                        aaGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGCCAAag                   >  1:8358921/1‑93 (MQ=255)
                                                                        aaGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGCCAAag                   >  1:4517247/1‑93 (MQ=255)
                                                                        aaGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCCCCAAag                   >  2:6190164/1‑93 (MQ=255)
                                                                         aGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGCCAAaga                  <  2:6080219/93‑1 (MQ=255)
                                                                         aGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGCCAAaga                  >  1:7286325/1‑93 (MQ=255)
                                                                         aGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGCCAAaga                  <  2:4317703/93‑1 (MQ=255)
                                                                         aGACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGCCAAaga                  >  2:4375883/1‑93 (MQ=255)
                                                                          gACAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGCCAAagag                 >  2:6651986/1‑93 (MQ=255)
                                                                           aCAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGCCAAAGAGc                <  1:651804/93‑1 (MQ=255)
                                                                           aCAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGCCAAAGAGc                <  1:196777/93‑1 (MQ=255)
                                                                            cAAAGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGCCAAAGAGCa               >  1:2507621/1‑93 (MQ=255)
                                                                             aaaGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGCCAAAGAGCAt              >  1:3687350/1‑93 (MQ=255)
                                                                               aGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGCCAAAGAGCATcc            >  2:8654089/1‑93 (MQ=255)
                                                                               aGCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGCCAAAGAGCATcc            >  1:8223401/1‑93 (MQ=255)
                                                                                gCCAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGCCAAAGAGCATCCa           <  2:128285/93‑1 (MQ=255)
                                                                                  cAAAGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGCCAAAGAGCATCCAgg         >  1:6084029/1‑93 (MQ=255)
                                                                                   aaaGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGCCAAAGAGCATCCAGGc        >  2:6382528/1‑93 (MQ=255)
                                                                                    aaGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGCCAAAGAGCATCCAGGCt       >  1:8797746/1‑93 (MQ=255)
                                                                                    aaGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGCCAAAGAGCATCCAGGCt       >  2:5113241/1‑93 (MQ=255)
                                                                                     aGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGCCAAAGAGCATCCAGGCTg      >  2:4859332/1‑93 (MQ=255)
                                                                                     aGTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGCCAAAGAGCATCCAGGCTg      <  1:2197806/93‑1 (MQ=255)
                                                                                      gTCGCGCTTTTCTCCTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGCCAAAGAGCATCCAGGCTgg     <  1:819635/93‑1 (MQ=255)
                                                                                      gTCGCGCTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGCCAAAGAGCATCCAGGCTgg     >  2:1402906/1‑93 (MQ=255)
                                                                                         gcgcTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGCCAAAGAGCATCCAGGCTGGGaa  <  1:3436634/93‑1 (MQ=255)
                                                                                            |                                                                                         
ACTACATTAATCAGGGAACGCCCGCCCAGTTAGGAGGTATGTTGGTGGATATGGCGGCGCGTCAGGTGAATAAAGACAAAGCCAAAGTCGCGTTTTTCTACTCAAGCCCCACCGTTACGGACCAAAACCAGTGGGTGAAAGAAGCGAAAGCGAAAATCGCCAAAGAGCATCCAGGCTGGGAA  >  NC_000913/1603451‑1603632

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