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

GCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGC  >  NC_000913/1376356‑1376646
                                                                                                                                                         |                                                                                                                                             
gCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAacac                                                                                                                                              >  2:410061/1‑151 (MQ=255)
  ttACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAacac                                                                                                                                              <  1:941966/149‑1 (MQ=255)
    aCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTc                                                                                                                                          <  2:1040130/151‑1 (MQ=255)
       gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGcc                                                                                                                                       <  1:55301/151‑1 (MQ=255)
           aTGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATg                                                                                                                                   <  1:1162336/151‑1 (MQ=255)
                  cgccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCccgccg                                                                                                                            <  2:660740/151‑1 (MQ=255)
                  cgccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCccgccg                                                                                                                            >  2:1214173/1‑151 (MQ=255)
                   gccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGt                                                                                                                           >  1:1029555/1‑151 (MQ=255)
                     cgcgTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGtt                                                                                                                          >  2:352782/1‑150 (MQ=255)
                               ggACGGCGCCCAAAAAATGGGGGGTAATCCCCTGCGCGTACTGGAACAGTCCGCCGCGCGTCCCGGACCGGTGGTGTATCTCAACCCTTAAACCCGTTATCAGTATAAAATGACCACCCCCCCCCCCGATGCCCGCCGTTACAtcgtcg                                                                                                                 <  1:841855/149‑1 (MQ=255)
                               ggACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGa                                                                                                                >  2:1337518/1‑150 (MQ=255)
                                            aaaaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                  >  2:301978/1‑151 (MQ=255)
                                            aaaaGGGGGAGTGATCCCCTGCGCGTACGGGAACAGGTCGCCGCGCGTACCGGACCCGTGGTGTATCTCGAACCGTAAAACCGCTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTaaa                                                                                                   <  2:214254/150‑1 (MQ=255)
                                             aaaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                  >  1:1256511/1‑150 (MQ=255)
                                              aaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACAt                                                                                                >  2:922008/1‑151 (MQ=255)
                                                 ggggagttaTCCCCGGCGGGTTCTGAAACGGGCCCCCGCGCGTCCCGGACCCGTGTTGTATCCCGACCCGTTAACCCCCTATCAGGATCATATGACCACCCCCCCCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt                                                                                              <  1:987282/147‑1 (MQ=255)
                                                    gagttaTTCCCGGCGCTTTCGGAAACAGGTCGCCGCGCGTACGGGACCGGGGGTGTATCTCGACCCGTTAACCCGCTATGATTATTATATAATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACa                                                                                           <  1:1259536/150‑1 (MQ=255)
                                                              cTGCGGTTACTGGAACAGGCCCCCGCGCGACCCGAACCGGGGGTTTACCTCGAACCGTTAACCCCCTACCAGAATCAAATGACCACCCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGc                                                                                 <  1:175820/150‑1 (MQ=255)
                                                              cTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                >  2:1073180/1‑151 (MQ=255)
                                                                gcgcGTACGGGAACAGGGCGCCGCCCGTACGGGACCCGTGGTTTATCTCAACCCTTTAACCCGTTATCAGGATAATATTTACAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                <  2:734457/149‑1 (MQ=255)
                                                                 cgcgTACGGGAACAGGCCCCCGCGCGTACCGGAACCGTGGTTTATCTCGAACCGTTAAACCGTTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGa                                                                               <  2:307217/149‑1 (MQ=255)
                                                                   cgTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAttt                                                                            >  1:18223/1‑150 (MQ=255)
                                                                     tatggcaaagggcccccgcGCGTGCGGCAGCGGTGGTTTATCCCAACCCTTTAACCCGCTATAAGGTTCATAGGATCACCACCCCCGCAGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                          <  1:1148730/135‑1 (MQ=255)
                                                                     tagggaaaCAGGTCCCCGCGCGTCCCGCACCGGGGTTTTATAACGAACCTTAAACCCGCTATAAGAATCTAATGATCACCACCCTCGCTGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGa                                                                               <  2:831379/139‑1 (MQ=255)
                                                                      aCTGGAACAGGTCGCCGCGCGTACCGGACCCGTGGTGTATCCCGAACCGTTAACCCGCTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTAt                                                                        <  2:881059/151‑1 (MQ=255)
                                                                           aaCAGGCCCCGGCGCGTCCCGGAACGGGGTGGTTCCCCGACCCGTAAACCCGTTATCAGGATAATATGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcaca                                                                    <  1:262475/150‑1 (MQ=255)
                                                                               ggTCGCCGCGCGTACCGGACCCGTGGTGTATCCCGACCCGTAAACCCGTTAAAAGGATAATATGACAACCCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAat                                                               <  1:770268/151‑1 (MQ=255)
                                                                                         cgTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                       >  2:15524/1‑149 (MQ=255)
                                                                                            cccgaacccGTGGTGTATCTCAAACCTTTAACCCGTTATCAGGATCATATAATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATa                                                  <  2:583372/144‑1 (MQ=255)
                                                                                              cggaagccttggttTACCTCGACCGTTAAACCCCTTTTCAGATAAATATGATCACCCCCCCCGCCTATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAata                                                              <  2:148122/129‑1 (MQ=255)
                                                                                               ggAACCGGGGTTTATCCCGACCCTTTAACCCGCTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATa                                                  <  2:963792/148‑1 (MQ=255)
                                                                                                 aaCGGTGGTGTATCTCGAACCGTTAACCCGCTATCAGGATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTg                                             <  2:598007/151‑1 (MQ=255)
                                                                                                  aCCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTgg                                            >  2:733757/1‑151 (MQ=255)
                                                                                                   ccGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgc                                           >  2:566902/1‑151 (MQ=255)
                                                                                                      gtggtgTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                         >  1:88960/1‑151 (MQ=255)
                                                                                                      ggggtgTCTCCCAAGCCGTTAACCCGCTATAAGAAAAATATAACCACCCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATaa                                                 <  1:135362/141‑1 (MQ=255)
                                                                                                       gtggtatacctccaCCGTTAAACCCATTATCAGTATAATAGGATCAACCCCCCCTCCTTGTCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCt                                              <  1:591820/133‑1 (MQ=255)
                                                                                                          tgtATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                     >  2:778968/1‑151 (MQ=255)
                                                                                                            tATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgcgc                                   >  2:989522/1‑151 (MQ=255)
                                                                                                             aTCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCt                                  >  1:1164317/1‑151 (MQ=255)
                                                                                                                tcGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTg                                 >  2:162605/1‑149 (MQ=255)
                                                                                                                 ccgaacGTTAAACCCGATATCGGGATAAATGAACAACCCCCCCCGCCGTGCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCa                               <  1:873306/146‑1 (MQ=255)
                                                                                                                  caaCCCTTAAACCCCTTATAGGAATAATAGGAACACCCCCCCCCCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCt                                  <  1:1055213/146‑1 (MQ=255)
                                                                                                                   aaCCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                             >  1:1158442/1‑150 (MQ=255)
                                                                                                                    aCCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc                           >  1:643022/1‑151 (MQ=255)
                                                                                                                    aCCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc                           >  1:146476/1‑151 (MQ=255)
                                                                                                                      cGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACa                          >  1:718930/1‑150 (MQ=255)
                                                                                                                      cGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACa                          >  1:150663/1‑150 (MQ=255)
                                                                                                                           aacccccTACCAGAATCATATGAACACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAc                        <  1:73878/142‑1 (MQ=255)
                                                                                                                            aacccgttacaaGGATAATAGAACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcg                      <  1:298915/139‑1 (MQ=255)
                                                                                                                            aCCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                    >  2:141759/1‑150 (MQ=255)
                                                                                                                             agccgttacaaGGAACATATGACAAACCCCCTCGCGGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc                     <  2:345025/140‑1 (MQ=255)
                                                                                                                              ccGCTATCAGTATAATATGACCAACCCCCTCGCCTATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc                     <  1:4245/148‑1 (MQ=255)
                                                                                                                               cGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAcc                 >  1:651618/1‑151 (MQ=255)
                                                                                                                                gCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt                >  2:499497/1‑151 (MQ=255)
                                                                                                                                       aGGATCAAATGACAACCCCCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTg               <  2:902315/146‑1 (MQ=255)
                                                                                                                                        ggATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGt         >  2:474474/1‑151 (MQ=255)
                                                                                                                                        ggATCATATGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc            <  2:706704/148‑1 (MQ=255)
                                                                                                                                        ggATCATATAATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCt             <  1:1311091/147‑1 (MQ=255)
                                                                                                                                          aTCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc        >  1:970803/1‑150 (MQ=255)
                                                                                                                                          aTCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc        >  2:937456/1‑150 (MQ=255)
                                                                                                                                            cATAGGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtg     <  2:542471/151‑1 (MQ=255)
                                                                                                                                             atatGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc            <  2:731656/143‑1 (MQ=255)
                                                                                                                                             atatGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgt    >  2:682712/1‑151 (MQ=255)
                                                                                                                                             atatGACCACCCCCCCCCCCGACGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc            <  1:210828/143‑1 (MQ=255)
                                                                                                                                              tatGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgtg   >  1:307116/1‑151 (MQ=255)
                                                                                                                                              tatGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgtg   >  1:594021/1‑151 (MQ=255)
                                                                                                                                              tatGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtg     >  2:941619/1‑149 (MQ=255)
                                                                                                                                               atGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc  >  1:212818/1‑151 (MQ=255)
                                                                                                                                                tGATCAACCCCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgtg   <  2:836667/149‑1 (MQ=255)
                                                                                                                                                         |                                                                                                                                             
GCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGC  >  NC_000913/1376356‑1376646

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