Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,376,5050AC42.4% 113.5 / 13.8 66T158P (ACC→CCC) ycjRputative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (33/5);  new base C (0/28);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.39e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTG  >  NC_000913/1376359‑1376645
                                                                                                                                                   |                                                                                                                                             
tACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCt                                                                                                                                           >  2:606268/1‑151 (MQ=255)
  ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTc                                                                                                                                         <  2:703781/150‑1 (MQ=255)
     ccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCg                                                                                                                                     <  1:558388/151‑1 (MQ=255)
        aggacctccgggcgTACCTTGGAGGTCGGCCGCTAAAAGGTGAGGGAGGCCCTGGGCGTACCGGAGCAGCCGGCCGCCACTGCCGGGCGCGTTGTGTATTTCGAACCGTTAAACCGCTATCAGGCTAGTAGGATCAACACCCTCGCCGAt                                                                                                                                   <  1:128575/139‑1 (MQ=255)
        aTGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATg                                                                                                                                  >  1:807458/1‑151 (MQ=255)
          gACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCACGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGcc                                                                                                                                <  2:65135/151‑1 (MQ=255)
                 ccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGtt                                                                                                                         >  2:806817/1‑151 (MQ=255)
                   gcgTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTAc                                                                                                                       >  2:998989/1‑151 (MQ=255)
                          atGGACAGCGACAGCAAAATGGTGAGTGATTCCCTGCGCGTACGGGAACAGGTCGCCGCGCGTACCGGAACCGTGGAGTATCTCGAACCGTTAAACCGCTATAAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTAAAtcgtcg                                                                                                                >  1:877202/2‑151 (MQ=255)
                                  cgcccGAAAAATGGTGAGTGATCCCCGGCGCGTACTGGAACGGGTCGCCGCGCTTCCCGGAACCGTGGTGTTTCTCGACCCGTAAAACCGCTATCAGGATCATATGACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGa                                                                                                         <  1:978722/147‑1 (MQ=255)
                                  cGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGAt                                                                                                        >  2:1049901/1‑151 (MQ=255)
                                       acAAAAGGGTGAGTGATCCCCTGCGCGTACTGGAACGGGTCGCCGCGCGTCCCGGACCCGTGGTGTATCTCGACCCTTTAACCCTCTATCAGAATAATATGACCAACCCCCTCCCCGATGCCCGCCGTTACATCGTCGAAAACGa                                                                                                         <  1:1085938/144‑1 (MQ=255)
                                         aacaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                 >  2:796544/4‑151 (MQ=255)
                                         aaaaTGGTGGGTTTTCCCCTGCGCGTACTGGAACAGTTCGCGGCGCGTCCGGGACCCGTGTTTTACCCCAACCCTTAAACCCGATAAAAGAATAATATGATCAACCCCCCCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTaaa                                                                                                  <  2:250919/150‑1 (MQ=255)
                                           aaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACAt                                                                                               >  2:717547/1‑151 (MQ=255)
                                              ggTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTCCGGGAACGGTGGTTTATCTCAACCCGTTAACCCGCTACCAGAATAATATGACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTa                                                                                            <  2:138960/151‑1 (MQ=255)
                                              ggTGAGTGATTCCCTGCGCGTACGGGAACAGGTCGCCGCGCGTACCGGACCCGTGGTGTATCCCGACCCGTAAAACCGCTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATg                                                                                              <  2:526234/149‑1 (MQ=255)
                                                tgagtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                           >  1:895095/1‑150 (MQ=255)
                                                      aTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                    >  1:671089/1‑151 (MQ=255)
                                                       ttCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATc                                                                                   >  1:381070/1‑151 (MQ=255)
                                                       tccccTGCGGGTAGTGAAACGGGTCGCCGCGCGTACCGAACCCGTGGTTTATCTCGACCCGTTAACCCGCTATCAGGATCATAGGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGa                                                                                        <  2:861011/144‑1 (MQ=255)
                                                        tCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCg                                                                                  >  2:235716/1‑151 (MQ=255)
                                                         cccTGCGCGTACTGGAACAGGTCGCCGCGCGTCCCGGACCGGTGGTGTATCTCGAACCTTAAACCCGCTACCAGGATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATc                                                                                   <  2:901312/149‑1 (MQ=255)
                                                             gcgcGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAt                                                                             >  2:977213/1‑151 (MQ=255)
                                                                 gTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                         >  2:360589/1‑151 (MQ=255)
                                                                   aCTGAAACGGGTCGCCGCGCGTCCGGGAACGGTGGTGTATCTCAAACCTTTAACCCGCTATAAGGATCATATGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                         <  2:742631/149‑1 (MQ=255)
                                                                    cTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                      >  2:31202/1‑151 (MQ=255)
                                                                          cAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATg                                                                 >  2:202118/1‑150 (MQ=255)
                                                                            ggcgcccgcGCGTCCGGGAACGGGGGTGTACCCCAACCCGTTAACCCCCTATCGGGATCATATGACCACCCCCCCCGCTGATGGCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa                                                                <  1:282682/144‑1 (MQ=255)
                                                                            ggTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAat                                                              >  2:1030949/1‑151 (MQ=255)
                                                                               cgccggGCGTACCGGACCCGTGGTGTATCTCAACCCGTAAACCCGCTATCAGGATAATATGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAata                                                             <  2:972339/149‑1 (MQ=255)
                                                                               cgccgcGCGAACCGGAACCGTGTGTTGTCTCGAACCGTAAAACCGCTTTCAGGATCATTTGACAACCCCCTCCGCCCATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAatat                                                            <  2:295966/150‑1 (MQ=255)
                                                                                     gcgTCCCGGAACCGTGGTGTATCTCAACCCTTAAACCCGCTATCAGGATCATATGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                      <  2:1084651/150‑1 (MQ=255)
                                                                                        tACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                      >  2:674683/1‑147 (MQ=255)
                                                                                        tACCGAACCCGTGTTGTATCTCAACCCGTTAACCCCCTAACAGGTTCATATGATCAACCCCCCCCCCGGTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAagaaga                                                   <  1:410327/150‑1 (MQ=255)
                                                                                         aCCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                  >  2:773104/1‑150 (MQ=255)
                                                                                          cgggACCCGTGGTTTATCTCACCCCGTTAACCCGCTTTCAGGCCAAAATAACCACCACCAGCGCCGACCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                  <  1:320849/147‑1 (MQ=255)
                                                                                          cggaaaCGGTGGTGTATCTCGACCCGTTAACCCGCTATAAGGATCATATGAAAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGaagaag                                                    <  2:34996/143‑1 (MQ=255)
                                                                                          ccGGACCCGTGGTGTATCTCGACCCGTAAAACCGCTATAAGGATCATATGACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                  <  2:344282/149‑1 (MQ=255)
                                                                                            ggAACGGGGGTGTCTCCCCAACCGTAAACCCCCTATCAGAAAAATATAATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                  <  1:344428/147‑1 (MQ=255)
                                                                                              aaCGGTGTTGTATCTCGAACCGTTAACCCGCTATCAGGATCATATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                  <  1:557908/145‑1 (MQ=255)
                                                                                              aaCCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTg                                            >  2:905963/1‑151 (MQ=255)
                                                                                                ccGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgc                                          >  1:707414/1‑151 (MQ=255)
                                                                                                  gtggtgTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                        >  1:133103/1‑151 (MQ=255)
                                                                                                  gtggtgTATCGCGACCCGTTAACCCGCTATGAGGAACATATAACCACCCCCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcg                                         <  1:67261/150‑1 (MQ=255)
                                                                                                    ggtgTATCGCAACCCTTAAACCCGCTATCAGGATAATATGAACAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                        <  2:626538/149‑1 (MQ=255)
                                                                                                    ggtgTATCCCAAACCTTAAACCCGTTATCAGAATCATATGATAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAg                                      <  2:647043/151‑1 (MQ=255)
                                                                                                     gttgtttcccgaCCGTTAAACCCCCTTTCAGGATCATATAATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCt                                             <  1:954338/133‑1 (MQ=255)
                                                                                                        tATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgcgc                                  >  1:92381/1‑151 (MQ=255)
                                                                                                           cccGACCCGTAAAACCGCTATCAGGATCATATGAAAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                    <  2:548282/144‑1 (MQ=255)
                                                                                                            tcGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGc                               >  1:754988/1‑150 (MQ=255)
                                                                                                            cccaaCCGTTAAACCCCCTAACAGTAACATATAACAAACCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                    <  1:605998/143‑1 (MQ=255)
                                                                                                              gAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                            >  1:136315/1‑151 (MQ=255)
                                                                                                               aaCCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                            >  2:631280/1‑150 (MQ=255)
                                                                                                                aCCCTTAAACCCGCTATCAGGATAATATGACCAACCCCCCCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCt                                 <  2:872414/145‑1 (MQ=255)
                                                                                                                  cGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACa                         >  2:141615/1‑150 (MQ=255)
                                                                                                                       aacccGCTATCAGGATCATATAAACAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAcc                      <  2:1056496/146‑1 (MQ=255)
                                                                                                                       aacccGCTATCAGAATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcg                     <  2:131036/147‑1 (MQ=255)
                                                                                                                            gCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt               >  1:602042/1‑151 (MQ=255)
                                                                                                                                 agggATCATATGACCAACCCCCCCCCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg          <  1:310134/149‑1 (MQ=255)
                                                                                                                                  aGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTCAACATGTCCAGATTATCGGCGATTTCTATCACATGAATATCGAAGGAGACAACCTGGCGCATGCCCTGCACGACAACCGCGACCTGCTc           >  1:366055/1‑149 (MQ=255)
                                                                                                                                  aGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCgg         >  2:188148/1‑151 (MQ=255)
                                                                                                                                  aGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCgg         >  2:582581/1‑151 (MQ=255)
                                                                                                                                        atatGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgt   >  2:860412/1‑151 (MQ=255)
                                                                                                                                         tatGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgtg  >  2:112646/1‑151 (MQ=255)
                                                                                                                                             acaaaCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgtg  <  2:1050488/144‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
TACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTG  >  NC_000913/1376359‑1376645

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