Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,311,859 A→G S62P (TCT→CCT)  alsB ← D‑allose ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,311,8590AG86.8% 1213.3 / 151.1 394S62P (TCT→CCT) alsBD‑allose ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (32/20);  new base G (203/139);  total (235/159)
Fisher's exact test for biased strand distribution p-value = 8.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

Only 100 of 405 total aligned reads displayed.

TTCACTGAGGATAATGGAGCGAAGGCGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAGGGTTTTCAATACGAC  >  NC_000913/4311759‑4311958
                                                                                                    |                                                                                                   
 tCACTGAGGATAATGGAGCGAAGGCGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGagag                                                                                                    <  1:7945556/101‑1 (MQ=255)
  cACTGAGGATAATGGAGCGAAGGCGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAgggg                                                                                                   >  2:7946022/1‑101 (MQ=255)
    cTGAGGATAATGGAGCGAAGGCGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCa                                                                                                 <  1:7945542/101‑1 (MQ=255)
      gAGGATAATGGAGCGAAGGCGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCaaa                                                                                               >  1:7946029/1‑101 (MQ=255)
       aGGATAATGGAGCGAAGGCGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCaaaa                                                                                              >  2:7946033/1‑101 (MQ=255)
        ggATAATGGAGCGAAGGCGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCaaaaa                                                                                             <  2:7945911/101‑1 (MQ=255)
          aTAATGGAGCGAAGGCGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAata                                                                                           <  2:7945380/101‑1 (MQ=255)
           tAATGGAGCGAAGGCGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAatat                                                                                          <  1:7945338/101‑1 (MQ=255)
           tAATGGAGCGAAGGCGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAatat                                                                                          >  2:7946038/1‑101 (MQ=255)
             aTGGAGCGAAGGCGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCa                                                                                        <  2:7945583/101‑1 (MQ=255)
              tGGAGCGAAGGCGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCaa                                                                                       <  1:7945469/101‑1 (MQ=255)
              tGGAGCGAAGGCGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCaa                                                                                       >  2:7946044/1‑101 (MQ=255)
              tGGAGCGAAGGCGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCaa                                                                                       <  2:7945709/101‑1 (MQ=255)
                gAGCGAAGGCGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACg                                                                                     >  1:7946046/1‑101 (MQ=255)
                 aGCGAAGGCGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGc                                                                                    <  1:7945793/101‑1 (MQ=255)
                   cGAAGGCGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTg                                                                                  <  2:7945413/101‑1 (MQ=255)
                    gAAGGCGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGa                                                                                 <  2:7945477/101‑1 (MQ=255)
                    gAAGGCGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGa                                                                                 >  1:7946057/1‑101 (MQ=255)
                    gAAGGCGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGa                                                                                 >  1:7946052/1‑101 (MQ=255)
                    gAAGGCGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGa                                                                                 >  2:7946054/1‑101 (MQ=255)
                        gCGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGcc                                                                             >  1:7946064/1‑101 (MQ=255)
                          gATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCa                                                                           >  1:7946069/1‑101 (MQ=255)
                           aTACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAg                                                                          <  1:7945021/101‑1 (MQ=255)
                           aTACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAg                                                                          >  2:7946072/1‑101 (MQ=255)
                            tACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAgt                                                                         >  1:7946076/1‑101 (MQ=255)
                             aCCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAgtg                                                                        <  2:7945390/101‑1 (MQ=255)
                              ccTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAgtgt                                                                       >  2:7946080/1‑101 (MQ=255)
                               cTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGtt                                                                      <  1:7945596/101‑1 (MQ=255)
                               cTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGtt                                                                      >  1:7946083/1‑101 (MQ=255)
                                tttGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGttt                                                                     >  1:7946087/1‑101 (MQ=255)
                                tttGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGttt                                                                     <  2:7945248/101‑1 (MQ=255)
                                   gTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTg                                                                  <  1:7945458/101‑1 (MQ=255)
                                   gTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTg                                                                  <  1:7945724/101‑1 (MQ=255)
                                    tAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGc                                                                 >  2:7946092/1‑101 (MQ=255)
                                     aaTTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCt                                                                >  2:7946097/1‑101 (MQ=255)
                                      aTTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCtt                                                               >  1:7946099/1‑101 (MQ=255)
                                       tttttATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTc                                                              <  2:7945451/101‑1 (MQ=255)
                                         tttATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCAt                                                            >  1:7946105/1‑101 (MQ=255)
                                         tttATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCAt                                                            >  1:7946109/1‑101 (MQ=255)
                                         tttATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCAt                                                            >  2:7946111/1‑101 (MQ=255)
                                           tattaCTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCt                                                          <  1:7945664/101‑1 (MQ=255)
                                            attaCTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCtt                                                         >  2:7946116/1‑101 (MQ=255)
                                             ttaCTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTc                                                        >  1:7946121/1‑101 (MQ=255)
                                              taCTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCa                                                       >  1:7946128/1‑101 (MQ=255)
                                              taCTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCa                                                       <  2:7945530/101‑1 (MQ=255)
                                              taCTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCa                                                       >  2:7946125/1‑101 (MQ=255)
                                               aCTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCaa                                                      >  2:7946131/1‑101 (MQ=255)
                                                   agaTCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGcc                                                  >  1:7946137/1‑101 (MQ=255)
                                                    gaTCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCt                                                 >  2:7946140/1‑101 (MQ=255)
                                                     aTCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCtt                                                >  2:7946145/1‑101 (MQ=255)
                                                        ttCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCttttt                                             >  2:7946153/1‑101 (MQ=255)
                                                        ttCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCttttt                                             >  1:7946150/1‑101 (MQ=255)
                                                         tCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCtttttt                                            >  1:7946156/1‑101 (MQ=255)
                                                           aaaTAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCa                                          >  1:7946161/1‑101 (MQ=255)
                                                           aaaTAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCa                                          >  2:7946163/1‑101 (MQ=255)
                                                            aaTAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCat                                         <  1:7945804/101‑1 (MQ=255)
                                                            aaTAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCat                                         <  1:7945702/101‑1 (MQ=255)
                                                             aTAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCata                                        >  1:7946169/1‑101 (MQ=255)
                                                              tAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCatat                                       <  1:7945931/101‑1 (MQ=255)
                                                              tAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCatat                                       >  1:7946174/1‑101 (MQ=255)
                                                               aaCTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATc                                      >  2:7946177/1‑101 (MQ=255)
                                                                aCTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCt                                     >  2:7946181/1‑101 (MQ=255)
                                                                aCTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCt                                     <  2:7945741/101‑1 (MQ=255)
                                                                  tGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTAc                                   <  1:7945495/101‑1 (MQ=255)
                                                                  tGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTAc                                   <  2:7945600/101‑1 (MQ=255)
                                                                  tGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTAc                                   >  2:7946185/1‑101 (MQ=255)
                                                                   gCAATTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTAcc                                  <  2:7945720/101‑1 (MQ=255)
                                                                     aaTTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCa                                >  2:7946193/1‑101 (MQ=255)
                                                                     aaTTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCa                                >  1:7946191/1‑101 (MQ=255)
                                                                      aTTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCaa                               >  2:7946198/1‑101 (MQ=255)
                                                                      aTTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCaa                               <  1:7945504/101‑1 (MQ=255)
                                                                      aTTGAGATTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCaa                               <  1:7945632/101‑1 (MQ=255)
                                                                         gagaTTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAAt                            >  2:7946204/1‑101 (MQ=255)
                                                                          agaTTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCaaaa                             <  1:7945446/99‑1 (MQ=255)
                                                                          agaTTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATg                           <  2:7945611/101‑1 (MQ=255)
                                                                           gaTTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATgg                          <  2:7945514/101‑1 (MQ=255)
                                                                           gaTTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATgg                          <  1:7945366/101‑1 (MQ=255)
                                                                            aTTGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATggg                         >  1:7946208/1‑101 (MQ=255)
                                                                            aTTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATggg                         >  2:7946210/1‑101 (MQ=255)
                                                                             ttGAAAATCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGt                        <  2:7945549/101‑1 (MQ=255)
                                                                               gAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTg                      >  1:7946216/1‑101 (MQ=255)
                                                                                aaaaTCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTgg                     <  2:7945648/101‑1 (MQ=255)
                                                                                aaaaTCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTgg                     >  2:7946218/1‑101 (MQ=255)
                                                                                aaaaTCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTgg                     >  1:7946221/1‑101 (MQ=255)
                                                                                 aaaTCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGga                    <  2:7945572/101‑1 (MQ=255)
                                                                                 aaaTCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGga                    >  1:7946225/1‑101 (MQ=255)
                                                                                  aaTCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGgag                   <  1:7945643/101‑1 (MQ=255)
                                                                                  aaTCGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGgag                   <  2:7945769/101‑1 (MQ=255)
                                                                                     cGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAgg                >  2:7946233/1‑101 (MQ=255)
                                                                                     cGCCTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAgg                >  1:7946232/1‑101 (MQ=255)
                                                                                      gCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAggg               >  1:7946238/1‑101 (MQ=255)
                                                                                       ccTTCTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAGGGt              >  2:7946243/1‑101 (MQ=255)
                                                                                        cTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAggg               >  1:7946248/1‑99 (MQ=255)
                                                                                        cTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAGGGtt             <  2:7945682/101‑1 (MQ=255)
                                                                                         ttCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGt                        >  2:7946249/1‑89 (MQ=255)
                                                                                           cTGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAGGGTTTTc          <  2:7945800/101‑1 (MQ=255)
                                                                                            tGAAGGAGGGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAGGGTTTTCa         >  2:7946255/1‑101 (MQ=255)
                                                                                            tGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGgaga                  >  1:7946258/1‑92 (MQ=255)
                                                                                                 gaggGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAGGGTTTTCAATACg    <  2:7945631/101‑1 (MQ=255)
                                                                                                   gagGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAGGGTTTTCAATACGAc  >  1:7946265/1‑101 (MQ=255)
                                                                                                    |                                                                                                   
TTCACTGAGGATAATGGAGCGAAGGCGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAGGGTTTTCAATACGAC  >  NC_000913/4311759‑4311958

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