Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,275,6110AC39.5% 137.5 / 13.8 77E110D (GAA→GACgarDD‑galactarate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (36/10);  new base C (0/30);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 9.16e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTTACGCCGTGCGTGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGA  >  NC_000913/3275462‑3275751
                                                                                                                                                        |                                                                                                                                            
ccGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTTACGCCGTGCGTGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAAc                                                                                                                                             <  2:794158/151‑1 (MQ=255)
 cGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTTACGCCGTGCGTGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAAcc                                                                                                                                            <  1:376555/151‑1 (MQ=255)
  ggCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTTACGCCGTGCGGGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAAccc                                                                                                                                           <  2:245372/151‑1 (MQ=255)
   gCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTTACGCCGTGCGTGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCt                                                                                                                                          <  1:505702/151‑1 (MQ=255)
   gCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTTACGCCGTGCGTGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCt                                                                                                                                          <  1:443063/151‑1 (MQ=255)
          gTGAAATTATTCGTTATGGCGAAGTGATTGGTTACGCCGTGCGTGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTAccgcc                                                                                                                                   >  1:239286/1‑151 (MQ=255)
            gAAATTATTCGTTATGGCGAAGTGATTGGTTACGCCGTGCGTGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACcgccgc                                                                                                                                 >  2:120047/1‑151 (MQ=255)
              aaTTATTCGTTATGGCGAAGTGATTGGTTACGCCGTGCGTGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTg                                                                                                                               >  1:382026/1‑151 (MQ=255)
               attattCGTTATGGCGAAGTGATTGGTTACGCCGTGCGTGCAATCCCACGCGAAAGCTGGATCGCCGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGAAACCAAAGTCCCGGAACCCTTACCGCCGCTgg                                                                                                                              >  2:812451/1‑151 (MQ=255)
                ttattCGTTATGGCGAAGTGATTGGTTACGCCGTGCGTGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTgg                                                                                                                              >  1:876718/1‑150 (MQ=255)
                       ttATGGCGAGGTGTTGGTTCGCGCCGGGGGTGCAACCCCCCGGGAAAGCTGATCCAACAAATAAATGTTTGTACTCCGGGAAGCGCCGCCTTTACACACGCTGCCACTGGCAACAAAAGTCCCGGACCCCTTACCGCCGCTGGaa                                                                                                                            <  2:622347/145‑1 (MQ=255)
                       ttATGGCGAAGTGATTGGTTACGCCGTGCGTGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATac                                                                                                                      >  2:732791/1‑151 (MQ=255)
                         aTGGCGAAGTGATTGGTTACGCCGTGCGTGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATacac                                                                                                                    >  1:116890/1‑151 (MQ=255)
                               aaGTGATTGGTTACGCCGTGCGTGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTg                                                                                                               >  1:782053/1‑150 (MQ=255)
                                 gTGATGGGTTACGCCGTGCGTGCAACCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTCCCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAg                                                                                                             <  1:538605/150‑1 (MQ=255)
                                       gtttACGCCGTGCGTGCAACCCCACGCGGAAGCTGGATCGACAAATAAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATc                                                                                                      <  1:600366/149‑1 (MQ=255)
                                         ttACGCCGTGCGTGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGc                                                                                                    >  1:43261/1‑151 (MQ=255)
                                          ttccGCCGGGGGGGCAATCCCACGGGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCCCCGTTACACACGCTGCCACTGGCAACAAAAGTCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCg                                                                                                     <  2:325081/147‑1 (MQ=255)
                                          tACGCCGTGCGTGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGc                                                                                                    >  1:1010175/1‑150 (MQ=255)
                                             cgcctggcgggcAACCCCCCGCGGAAGCTGGACCAACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACCCCCGCTGCCACTGCCACCCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCaa                                                                                                  <  2:65584/144‑1 (MQ=255)
                                             ccccgtccgtGCATCCCCACGGGGAAGCTGGATCGCCAATGCATTGGTTATATACCGGAAAGCCCCGCTGTTGCACCCGCTGCCACTGGCAACCAAATTCCCGGACCCTTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGc                                                                                                    <  1:1034077/146‑1 (MQ=255)
                                             cGCCGTGCGTGCAACCCCACGGGGAAGCTGGATCGCGAAACCAATGGTTGTACTCCCGGAAGCGCCGCCGTTACACCCGCTGCCACTGGCAACCAAAGTCCCGGACCCCTTCCCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAAt                                                                                                 <  1:88469/151‑1 (MQ=255)
                                              gCCGTGCGTGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATg                                                                                                >  1:728249/1‑151 (MQ=255)
                                               ccGTGCGTGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCACCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGc                                                                                               >  1:311459/1‑151 (MQ=255)
                                                 gtgcgtgaAACCCCACGGGGAAGCTGGATCGACGAATCAATGTTTGTACTCCGGGAGGCGCCGCCTTTACACACGCTGCCACGGGCAACCAAAGTCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCg                                                                                             <  1:671392/151‑1 (MQ=255)
                                                         aaTCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGc                                                                                     >  2:203712/1‑151 (MQ=255)
                                                             ccgcggggAAGCTGGACCGACGAATAAAGGTTGGTACCCCGGGAGGCGCCGCCGTTCCCCACGCTGCCAGTGGCACCCAAAGCCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTgg                                                                                 <  2:963690/145‑1 (MQ=255)
                                                             ccACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTgg                                                                                 >  2:930304/1‑151 (MQ=255)
                                                               aCGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACCCCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGc                                                                               >  2:318772/1‑151 (MQ=255)
                                                                 gcgGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACTCCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCAc                                                                             >  2:922212/1‑151 (MQ=255)
                                                                    gAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCa                                                                           >  2:523602/1‑150 (MQ=255)
                                                                     aaGGTGGATCGACGAATCATTGGTTGTACACCGGGAAGCGCCGCCGTTACACCCGCTGCCACTGGCACCCAAAGTCCGGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCaaa                                                                         <  1:24518/151‑1 (MQ=255)
                                                                     aaGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCaa                                                                          >  1:11860/1‑150 (MQ=255)
                                                                       gCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCaaaa                                                                        >  1:99741/1‑150 (MQ=255)
                                                                         tGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAAcc                                                                     >  2:568574/1‑151 (MQ=255)
                                                                          ggATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCt                                                                    >  1:763928/1‑151 (MQ=255)
                                                                          ggAGCGACAAATCAAGGGTTGTCCTCCGGAAGCCGCCCCCGTTCCACACGCTGCCATTGGCAACCAAAGTCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCaaaa                                                                        <  2:131060/147‑1 (MQ=255)
                                                                            aTCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGc                                                                  >  1:440619/1‑151 (MQ=255)
                                                                              cgacgaATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTc                                                                >  1:924618/1‑151 (MQ=255)
                                                                                   caTCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTAt                                                           >  2:707060/2‑151 (MQ=255)
                                                                                   aaTCAATGGTTGTACACCGGAAAGCGCCGCCGTTACCCACGCTGCCACGGGCAACCAAAGTCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTAt                                                           <  2:515698/151‑1 (MQ=255)
                                                                                       aaTGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATcacc                                                       >  1:300588/1‑151 (MQ=255)
                                                                                                actacCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTc                                              >  1:404416/1‑151 (MQ=255)
                                                                                                 ctacCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTcc                                             >  1:294542/1‑151 (MQ=255)
                                                                                                    ccGGAGCCGCCGCCTTTACACCCGCTGCCACTGGCACCCAAAGTCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTc                                              <  2:971105/147‑1 (MQ=255)
                                                                                                         aaGCCCCGCCTTTACCCACGCTGCCACGGGCAACCAAAGTCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACtgtg                                       <  2:53303/150‑1 (MQ=255)
                                                                                                             cgccgcCGTTACACCCGCTGCCACTGGCAACCAAAGTCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGc                                   <  1:288745/151‑1 (MQ=255)
                                                                                                              gccgccGTTACCCACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCCGGAAGGATAAACCTTTGAGGGCTATCGAAATGCCGATGTCTGCGTGGGCACCAAAAACCTGCTAGGTATCACCACCAGCGTCCACTGTGTGGc                                   >  2:853655/1‑150 (MQ=255)
                                                                                                                 gccgTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGc                               >  2:18612/1‑151 (MQ=255)
                                                                                                                   cgtaaccctcGCTGCCGCTGGCAACCAAAGCCCCGCACCCCTTCCCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAgg                                <  2:300395/139‑1 (MQ=255)
                                                                                                                     ttACCCACGCGCCCACTGGCAACCAAAGCCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCgt                             <  1:221523/149‑1 (MQ=255)
                                                                                                                     ttACACACTCTCCCACTGGCAACAAACGACCCGTATCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGc                                   <  2:814252/143‑1 (MQ=255)
                                                                                                                      tACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCgtgg                           >  2:113166/1‑150 (MQ=255)
                                                                                                                        cacacGCTGCCGCTGTCACCCAAAGTCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCgtggtg                         <  2:969714/150‑1 (MQ=255)
                                                                                                                          cccccTGCCCCTGGAAACCAAATCCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGTGCTATCGCAATGCCGATGGCAGCGTGTGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTTTGTGGCATGCttggtacac                      <  2:198687/147‑5 (MQ=255)
                                                                                                                           ccGCTGCCACTGGCACCCAAAGTCCCGGACCCTTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACt                     <  2:381868/150‑1 (MQ=255)
                                                                                                                             gCTGCCACTGGCAACCACAGTCCCGGAACCCTTACCGCCGCCGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGCGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGCGGCAGGCGCGGCGGACTa                    >  2:104882/1‑150 (MQ=255)
                                                                                                                              cggcaaCTGCCAACCAAATCCCGGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTTTGGCAGGCGTGGTGGACTAt                   <  2:99712/145‑1 (MQ=255)
                                                                                                                              cTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATg                  >  1:377352/1‑151 (MQ=255)
                                                                                                                              cTGCCACTGGCAACCAAAGTCCCGAACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTAt                   <  1:732770/150‑1 (MQ=255)
                                                                                                                                gCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATg                  <  1:895728/149‑1 (MQ=255)
                                                                                                                                 ccACTGGCAACCAAAGCCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATg                  <  1:826311/148‑1 (MQ=255)
                                                                                                                                 caagtGGCAACCAAGGTCCCGGACCCTTTTCCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCCCTGTGTGGCAGGCGTGGTGGACTAt                   <  1:105952/143‑1 (MQ=255)
                                                                                                                                  cgcTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATgtagt              <  1:94395/149‑1 (MQ=255)
                                                                                                                                      ggaaaCCAAACCCCCGACCCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATgtagt              <  2:1042838/144‑1 (MQ=255)
                                                                                                                                       gCAACGAAAGCCCCAGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCTCCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTATTaaaaa         <  2:919227/151‑1 (MQ=255)
                                                                                                                                       gCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTaaaa          >  1:259166/1‑150 (MQ=255)
                                                                                                                                         aaCCAAAGTCGCGAAACCGTTACCGCGGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATTACCACCATCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTaaaa          <  2:190238/148‑1 (MQ=255)
                                                                                                                                          aCCAAAGTCCCGAACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCa      <  1:651241/151‑1 (MQ=255)
                                                                                                                                           ccAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGACGTAAAAATCAt     >  2:796727/1‑151 (MQ=255)
                                                                                                                                           caaaaGTCACGGACCCCTTACCGCCGCTGGAAGGATACTCCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCGGCTCAGTATCACCAACAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAAt        <  2:239216/146‑1 (MQ=255)
                                                                                                                                             aaatgcccGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTg   <  2:116853/146‑1 (MQ=255)
                                                                                                                                             aaaGTCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCAtt    <  1:360940/150‑1 (MQ=255)
                                                                                                                                             aaaGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTg   >  2:444401/1‑151 (MQ=255)
                                                                                                                                             aaaGCCCCGAACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATc       <  2:300512/147‑1 (MQ=255)
                                                                                                                                              aaGTCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGa  <  1:807526/151‑1 (MQ=255)
                                                                                                                                               aGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGa  >  1:373768/1‑150 (MQ=255)
                                                                                                                                                        |                                                                                                                                            
CCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTTACGCCGTGCGTGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGA  >  NC_000913/3275462‑3275751

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