Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,376,5050AC49.2% 91.0 / 24.0 67T158P (ACC→CCC) ycjRputative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (31/3);  new base C (0/33);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 5.46e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCAT  >  NC_000913/1376355‑1376648
                                                                                                                                                          |                                                                                                                                               
cGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAaca                                                                                                                                                 <  2:1070741/151‑1 (MQ=255)
   ttACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAaca                                                                                                                                                 <  2:811400/148‑1 (MQ=255)
   ttACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACAccc                                                                                                                                              >  2:834795/1‑151 (MQ=255)
    tACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCt                                                                                                                                             >  1:73419/1‑151 (MQ=255)
       cgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGCGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCg                                                                                                                                           >  2:690495/1‑150 (MQ=255)
        gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGcc                                                                                                                                         <  2:499771/151‑1 (MQ=255)
          cgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCg                                                                                                                                        >  1:112386/1‑150 (MQ=255)
              gACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGcc                                                                                                                                   >  1:321781/1‑151 (MQ=255)
                   cgccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCccgcc                                                                                                                               <  2:177814/150‑1 (MQ=255)
                       gcgtggcgggggcGGCGACCGAAAAGTGGTGATTGCTCCCCTGCGCGTACGGAAACGGGCCGCCGCGCGTACGGGCCCGGTGGTGTACCCCAACCTTTAAACCCGTTATCAGAATCATATGACCACCCCCCTCGCCGATGCCCGCCGtt                                                                                                                            <  1:84560/137‑1 (MQ=255)
                       gcgTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTAc                                                                                                                          >  1:108209/1‑151 (MQ=255)
                       gcgTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCAGATGCCCGCCGTTAc                                                                                                                          >  2:490623/1‑151 (MQ=255)
                          tAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtc                                                                                                                       >  1:371094/1‑151 (MQ=255)
                             cgtGGACGGCGACCAAAAATTGTGGATTGATCCCCTGCGCTTACTGAAACGGGCCCCCGCCCGTCCGGGAGCCGTGTTTTTTCTCGACCTTTTAACCCGCTAGCAGAATAATATGATCACCCCCCCCGCCGATGCCCGCCGTTACAtcg                                                                                                                      <  1:953437/147‑1 (MQ=255)
                             ccgggACGGCGCCCGAAAAATGGTGAGTGATCCCCGGCGCGTACTGGAACAGGTCGCCGCGCGTCCCGAACCCGTGGTGTATCTCGACCCGTAAACCCGCTATCAGAATAATATGATAAACCCCCTCGCCGATGCCCGCCGTTACAtcgtc                                                                                                                    <  2:91797/148‑1 (MQ=255)
                                 gACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGaa                                                                                                                 >  2:527551/1‑150 (MQ=255)
                                    gccgcccaaaaaaGGGTGGGTATTTCCCGGCGCTTCCGGGAACGGGCCGCCGCGCGTACCGAAACGGTGGTTTATCCCAACCCTTAAACCCGTTACAAGAATAAAATGACAACCCCCCCCGCGGAGGCCCGCCGTTACATCGTCGAAAAc                                                                                                              <  1:638694/141‑1 (MQ=255)
                                     gCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACg                                                                                                             >  2:951702/1‑150 (MQ=255)
                                         ccaaaaaaGGGTGAGTTTTCCCCGGCGGGTACGGGAAGAGCCCCCCGCGCGTCCCAGAACTGGGGTTTATCTGAAAGAGTAAAACCGCTATCAGAATCATATGACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGa                                                                                                            <  2:407593/143‑1 (MQ=255)
                                         ccGAAAAAGGGTGGGTGATCCCCTGCGCGTACGGGAACGGGCCGCCGCGCGTACCGGACCGGTGGTGTATCCCAACCCGTTAACCCGCTATCAGGATAATAGGACAACCCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGa                                                                                                            <  1:1082458/147‑1 (MQ=255)
                                            cAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTaa                                                                                                      >  2:1016869/1‑150 (MQ=255)
                                               aaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACAt                                                                                                  >  1:290833/1‑151 (MQ=255)
                                                aTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATg                                                                                                 >  2:815217/1‑151 (MQ=255)
                                                 tGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt                                                                                                >  2:206770/1‑151 (MQ=255)
                                                   gggagtgaTCCCCTGCGCGTACTGAAACAGGCCGCCGCGCGTACCGGACCCGTGGTGTACCTCGACCCTTAAAACCGTTATAAGGATAATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                              <  1:1274732/149‑1 (MQ=255)
                                                        ttatccccGGCGCGTACTGGAACAGGTCCCCGCGCGTCCCGGACCCGTGGTTTATCCCAACCCGTAAACCCGCTATCAGGATAATATAATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAt                                                                                          <  1:204699/145‑1 (MQ=255)
                                                        tgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAtt                                                                                         >  2:710092/1‑151 (MQ=255)
                                                         gaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAtt                                                                                         >  1:215949/1‑150 (MQ=255)
                                                           ttCCCTGCGCGTACTGGACCAGGCCGCCGCGCGTCCCGGACCCGGGGTGTATCTCGACCCGTTAAACCCCTATCAGGATCATATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAtt                                                                                         <  2:460271/148‑1 (MQ=255)
                                                           ttCCCTGCGCGTACTGGAACGGGTCGCCGCGCGTCCCGGACCGGTGGTGTATCTCGACCCTTTAACCCGCTACCAGGATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                       <  2:606885/150‑1 (MQ=255)
                                                           tccccGGCGCGTACGGGAACAGGCCGCCGCGCGTCCGGGACCGGTGGTTTATCTCGACCCGTAAACCCGCTATCAGGATAATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAtt                                                                                         <  2:1168672/146‑1 (MQ=255)
                                                             cccGGCGCGTACGGGAACAGGTCGCCGCGCGTCCGGAAACGGTGGTGTATCTCGACCCTTTAACCCGCTATCAGGATAATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATc                                                                                      <  2:977305/149‑1 (MQ=255)
                                                              cttgggcttacggaaaCAGGCCGCCGCGCGCACGGAACCGGTGGTGTACCTCAACCCTTAAACCCCCAACAAGGATAAAATAACAACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAg                                                                                            <  1:840805/128‑1 (MQ=255)
                                                              ccggcgcGTACTGAAACGGGTCGCCGCGCGTCCCGGACCGGTGGTGTATCCCGACCCTTTAACCCGCTATCAGGATAATATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGc                                                                                   <  1:1096036/148‑1 (MQ=255)
                                                              ccTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGc                                                                                   >  1:165832/1‑151 (MQ=255)
                                                                         tGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                         >  1:1177705/1‑150 (MQ=255)
                                                                            aaCAGGTCGCCGCGCGTCCCGAACCGGGGGTGTATCTCAACCGTTTAACCCGCTATAAGGATAATATGACAACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                         <  2:1220236/147‑1 (MQ=255)
                                                                            aaCAGGCCCCCGCGCGTCCCGAACCGGTGTTTTTTCTCGACCCTTAAACCCGCTATCAGGATAATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcaca                                                                      <  2:544396/150‑1 (MQ=255)
                                                                                ggcccccgcGCGAACCGGACCCGTGGTTTATCTCGACCCGTTAACCCCTTATCAGGATAATATGCCCAACCCCCTCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcaca                                                                      <  1:1342255/141‑1 (MQ=255)
                                                                                ggTCGCCGCGCGTCCCGGACCGGTGGTGTATCTCGACCCTTTAAACCGCTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa                                                                   <  2:1304724/149‑1 (MQ=255)
                                                                                   cgccgcGCGTACCGGACCCGTGGTGTATCTCGAACCTTAAACCCGCTATCAGGATAATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAatat                                                               <  2:496906/150‑1 (MQ=255)
                                                                                   ccccgcGCGACCCGAACCTGTGGTTTTCCTCAACCCGTTAACCCGTTTTGAGGATCATTTGATCAACCCCCTCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAata                                                                <  1:1155472/147‑1 (MQ=255)
                                                                                     ccgcGCGTCCCGGAACGGTGGTGTCTCTCAACCCTTTAACCCGCTACAAGGATAATATGATAACCCCCCTCGCCATTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATc                                                              <  1:421640/149‑1 (MQ=255)
                                                                                      cgcgcgTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                            >  1:868373/1‑150 (MQ=255)
                                                                                        cgcgTCCGGGACCGGTGGTGTATCTCGACCTTTAAACCCGCTATCAGGATCATATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaag                                                          <  2:793153/150‑1 (MQ=255)
                                                                                            tACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                         >  1:298915/1‑147 (MQ=255)
                                                                                            tACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                         >  1:703442/1‑147 (MQ=255)
                                                                                              ccGAAACCGTGTTGTATCTCGACCCTTTAACCCGCTATCAGGATCATATAATAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATa                                                    <  1:970437/150‑1 (MQ=255)
                                                                                                 ggaccggtggtgTATCTCGACCCGTTAACCCGCTATCAGTATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAcc                                                 <  2:1064783/145‑1 (MQ=255)
                                                                                                 ggaactgggtggtATCCCAACCCGTTAACCCGCTATCAGGATTATATGTCCAACACCCTCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAc                                                  <  2:460136/139‑1 (MQ=255)
                                                                                                 gAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCt                                                >  2:541758/1‑151 (MQ=255)
                                                                                                       gtggtgTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                           >  1:835344/1‑151 (MQ=255)
                                                                                                           tgtATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                       >  1:210493/1‑151 (MQ=255)
                                                                                                           tgtATCTCGAACCGTAAACCCGCTATCAGGATCATAGGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAgg                                        <  2:104831/150‑1 (MQ=255)
                                                                                                            gtATCTCGACCCGTAAACCCCCTATCAGGATAATATGATCAACCCCCCCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                       <  2:148134/150‑1 (MQ=255)
                                                                                                                    aaCTGTTAACCCATTATCAGTATAATATAACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                               <  2:195384/150‑1 (MQ=255)
                                                                                                                     acccttaaaaCCCCTATCAGGATCATATGACCAACCCCCTCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCa                                 <  2:215952/141‑1 (MQ=255)
                                                                                                                     aCCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc                             >  2:264053/1‑151 (MQ=255)
                                                                                                                             aacccGCTATCAGGATCATATAACCAACCCCCGCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAcc                         <  1:613028/146‑1 (MQ=255)
                                                                                                                                     aTCAGAATAATATAACCACCCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTg                 <  1:831547/149‑1 (MQ=255)
                                                                                                                                       cAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg             >  2:814619/1‑151 (MQ=255)
                                                                                                                                       cAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg             >  2:983871/1‑151 (MQ=255)
                                                                                                                                          ggATCATATGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTg                 <  2:1197447/145‑1 (MQ=255)
                                                                                                                                          ggATAATATAACCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg             <  2:211264/149‑1 (MQ=255)
                                                                                                                                             tCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAt        >  1:1003199/1‑151 (MQ=255)
                                                                                                                                                tatGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc              <  1:206767/142‑1 (MQ=255)
                                                                                                                                                   gATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCat  >  2:250281/1‑151 (MQ=255)
                                                                                                                                                    aTCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCat  >  1:830528/1‑150 (MQ=255)
                                                                                                                                                     tCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtg       >  2:355463/1‑144 (MQ=255)
                                                                                                                                                          |                                                                                                                                               
CGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCAT  >  NC_000913/1376355‑1376648

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