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

TACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTG  >  NC_000913/1376359‑1376645
                                                                                                                                                   |                                                                                                                                             
tACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCt                                                                                                                                           >  2:606268/1‑151 (MQ=255)
  ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCT‑c                                                                                                                                         <  2:703781/150‑1 (MQ=255)
     ccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCT‑CGCCg                                                                                                                                     <  1:558388/151‑1 (MQ=255)
        aggacctccgggcgTACCTTGGAGGTCGGCCGCTAAAAGGTGAGGGAGGCCCTGGGCGTACCGGAGCAGCCGGCCGCCACTGCCGGGCGCGTTGTGTATTTCGAACCGTTAAA‑CCGCTATCAGGCTAGTAGGATCAACACCCT‑CGCCGAt                                                                                                                                   <  1:128575/139‑1 (MQ=255)
        aTGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCT‑CGCCGATg                                                                                                                                  >  1:807458/1‑151 (MQ=255)
          gACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCACGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCT‑CGCCGATGcc                                                                                                                                <  2:65135/151‑1 (MQ=255)
                 ccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCT‑CGCCGATGCCCGCCGtt                                                                                                                         >  2:806817/1‑151 (MQ=255)
                   gcgTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCT‑CGCCGATGCCCGCCGTTAc                                                                                                                       >  2:998989/1‑151 (MQ=255)
                          atGGACAGCGACAGCAAAATGGTGAGTGATTCCCTGCGCGTACGGGAACAGGTCGCCGCGCGTACCGGAACCGTGGAGTATCTCGAACCGTTAAA‑CCGCTATAAGGATCATATGATCAACACCCT‑CGCCGATGCCCGCCGTTAAAtcgtcg                                                                                                                >  1:877202/2‑151 (MQ=255)
                                  cgcccGAAAAATGGTGAGTGATCCCCGGCGCGTACTGGAACGGGTCGCCGCGCTTCCCGGAACCGTGGTGTTTCTCGACCCGTAAAA‑CCGCTATCAGGATCATATGACCAACCCCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGa                                                                                                         <  1:978722/147‑1 (MQ=255)
                                  cGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGAt                                                                                                        >  2:1049901/1‑151 (MQ=255)
                                       acAAAAGGGTGAGTGATCCCCTGCGCGTACTGGAACGGGTCGCCGCGCGTCCCGGACCCGTGGTGTATCTCGACCCTTTAAC‑CCTCTATCAGAATAATATGACCAACCCCCT‑CCCCGATGCCCGCCGTTACATCGTCGAAAACGa                                                                                                         <  1:1085938/144‑1 (MQ=255)
                                         aacaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                 >  2:796544/4‑151 (MQ=255)
                                         aaaaTGGTGGGTTTTCCCCTGCGCGTACTGGAACAGTTCGCGGCGCGTCCGGGACCCGTGTTTTA‑CCCCAACCCTTAAACCCGATAAAAGAATAATATGATCAACCCCCC‑CGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTaaa                                                                                                  <  2:250919/150‑1 (MQ=255)
                                           aaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACAt                                                                                               >  2:717547/1‑151 (MQ=255)
                                              ggTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTCCGGGAACGGTGGTTTATCTCAACCCGTTAAC‑CCGCTACCAGAATAATATGACCAACCCCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTa                                                                                            <  2:138960/151‑1 (MQ=255)
                                              ggTGAGTGATTCCCTGCGCGTACGGGAACAGGTCGCCGCGCGTACCGGACCCGTGGTGTATCCCGACCCGTAAAA‑CCGCTATCAGGATCATATGATCAACCCCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATg                                                                                              <  2:526234/149‑1 (MQ=255)
                                                tgagtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                           >  1:895095/1‑150 (MQ=255)
                                                      aTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                    >  1:671089/1‑151 (MQ=255)
                                                       ttCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATc                                                                                   >  1:381070/1‑151 (MQ=255)
                                                       tccccTGCGGGTAGTGAAACGGGTCGCCGCGCGTACCGAACCCGTGGTTTATCTCGACCCGTTAAC‑CCGCTATCAGGATCATAGGATCAACCCCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGa                                                                                        <  2:861011/144‑1 (MQ=255)
                                                        tCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCg                                                                                  >  2:235716/1‑151 (MQ=255)
                                                         cccTGCGCGTACTGGAACAGGTCGCCGCGCGTCCCGGACCGGTGGTGTATCTCGAACCTTAAAC‑CCGCTACCAGGATCATATGATCAACCCCCC‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATc                                                                                   <  2:901312/149‑1 (MQ=255)
                                                             gcgcGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAt                                                                             >  2:977213/1‑151 (MQ=255)
                                                                 gTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                         >  2:360589/1‑151 (MQ=255)
                                                                   aCTGAAACGGGTCGCCGCGCGTCCGGGAACGGTGGTGTATCTCAAACCTTTAAC‑CCGCTATAAGGATCATATGACCAACCCCCC‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                         <  2:742631/149‑1 (MQ=255)
                                                                    cTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                      >  2:31202/1‑151 (MQ=255)
                                                                          cAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATg                                                                 >  2:202118/1‑150 (MQ=255)
                                                                            ggcgcccgcGCGTCCGGGAACGGGGGTGTACCCCAACCCGTTAAC‑CCCCTATCGGGATCATATGACCACCCCCCC‑CGCTGATGGCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa                                                                <  1:282682/144‑1 (MQ=255)
                                                                            ggTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAat                                                              >  2:1030949/1‑151 (MQ=255)
                                                                               cgccggGCGTACCGGACCCGTGGTGTATCTCAACCCGTAAAC‑CCGCTATCAGGATAATATGACCACCCCCCC‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAata                                                             <  2:972339/149‑1 (MQ=255)
                                                                               cgccgcGCGAACCGGAACCGTGTGTTGTCTCGAACCGTAAAA‑CCGCTTTCAGGATCATTTGA‑CAACCCCCTCCGCCCATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAatat                                                            <  2:295966/150‑1 (MQ=255)
                                                                                     gcgTCCCGGAACCGTGGTGTATCTC‑AACCCTTAAACCCGCTATCAGGATCATATGACCACCCCCCC‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                      <  2:1084651/150‑1 (MQ=255)
                                                                                        tACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                      >  2:674683/1‑147 (MQ=255)
                                                                                        tACCGAACCCGTGTTGTATCTCAACCCGTTAAC‑CCCCTAACAGGTTCATATGATCAACCCCCC‑CCCCGGTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAagaaga                                                   <  1:410327/150‑1 (MQ=255)
                                                                                         aCCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                  >  2:773104/1‑150 (MQ=255)
                                                                                          cgggACCCGTGGTTTATCTCACCCCGTTAAC‑CCGCTTTCAGGCCAAAATAACCACCACCAG‑CGCCGACCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                  <  1:320849/147‑1 (MQ=255)
                                                                                          cggaaaCGGTGGTGTATCTCGACCCGTTAAC‑CCGCTATAAGGATCATATGAAAAACCCCCC‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGaagaag                                                    <  2:34996/143‑1 (MQ=255)
                                                                                          ccGGACCCGTGGTGTATCTCGACCCGTAAAA‑CCGCTATAAGGATCATATGACCAACCCCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                  <  2:344282/149‑1 (MQ=255)
                                                                                            ggAACGGGGGTGTCTCCCCAACCGTAAAC‑CCCCTATCAGAAAAATATAATCAACCCCCC‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                  <  1:344428/147‑1 (MQ=255)
                                                                                              aaCGGTGTTGTATCTCGAACCGTTAAC‑CCGCTATCAGGATCATATGATCACCCCCCC‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                  <  1:557908/145‑1 (MQ=255)
                                                                                              aaCCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTg                                            >  2:905963/1‑151 (MQ=255)
                                                                                                ccGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgc                                          >  1:707414/1‑151 (MQ=255)
                                                                                                  gtggtgTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                        >  1:133103/1‑151 (MQ=255)
                                                                                                  gtggtgTATCGCGACCCGTTAAC‑CCGCTATGAGGAACATATAACCACCCCCCG‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcg                                         <  1:67261/150‑1 (MQ=255)
                                                                                                    ggtgTATCGC‑AACCCTTAAACCCGCTATCAGGATAATATGAACAACCCCCC‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                        <  2:626538/149‑1 (MQ=255)
                                                                                                    ggtgTATCCCAAACCTTAAAC‑CCGTTATCAGAATCATATGATAAACCCCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAg                                      <  2:647043/151‑1 (MQ=255)
                                                                                                     gttgtttcccgaCCGTTAAACCCCCTTTCAGGATCATATAATCAACCCCCC‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCt                                             <  1:954338/133‑1 (MQ=255)
                                                                                                        tATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgcgc                                  >  1:92381/1‑151 (MQ=255)
                                                                                                           cccGACCCGTAAAA‑CCGCTATCAGGATCATATGAAAAACCCCCC‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                    <  2:548282/144‑1 (MQ=255)
                                                                                                            tcGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGc                               >  1:754988/1‑150 (MQ=255)
                                                                                                            cccaaCCGTTAAACCCCCTAACAGTAACATATAACAAACCCCCC‑CGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                    <  1:605998/143‑1 (MQ=255)
                                                                                                              gAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                            >  1:136315/1‑151 (MQ=255)
                                                                                                               aaCCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                            >  2:631280/1‑150 (MQ=255)
                                                                                                                aCCCTTAAACCCGCTATCAGGATAATATGACCAACCCCCC‑CGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCt                                 <  2:872414/145‑1 (MQ=255)
                                                                                                                  cGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACa                         >  2:141615/1‑150 (MQ=255)
                                                                                                                       aacccGCTATCAGGATCATATAAACAACCCCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAcc                      <  2:1056496/146‑1 (MQ=255)
                                                                                                                       aacccGCTATCAGAATCATATGATCAACCCCCC‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcg                     <  2:131036/147‑1 (MQ=255)
                                                                                                                            gCTATCAGGATCATATGATCAACACCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt               >  1:602042/1‑151 (MQ=255)
                                                                                                                                 agggATCATATGACCAACCCCCC‑CCCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg          <  1:310134/149‑1 (MQ=255)
                                                                                                                                  aGGATCATATGATCAACACCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTCAACATGTCCAGATTATCGGCGATTTCTATCACATGAATATCGAAGGAGACAACCTGGCGCATGCCCTGCACGACAACCGCGACCTGCTc           >  1:366055/1‑149 (MQ=255)
                                                                                                                                  aGGATCATATGATCAACACCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCgg         >  2:188148/1‑151 (MQ=255)
                                                                                                                                  aGGATCATATGATCAACACCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCgg         >  2:582581/1‑151 (MQ=255)
                                                                                                                                        atatGATCAACACCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgt   >  2:860412/1‑151 (MQ=255)
                                                                                                                                         tatGATCAACACCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgtg  >  2:112646/1‑151 (MQ=255)
                                                                                                                                             acaaaCCCCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgtg  <  2:1050488/144‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
TACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCT‑CGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTG  >  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: ‑