Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913298,9600AC48.9% 103.8 / 36.4 89Y659D (TAT→GAT) paoCPaoABC aldehyde oxidoreductase, Moco‑containing subunit
Reads supporting (aligned to +/- strand):  ref base A (32/13);  new base C (0/43);  total (32/57)
Fisher's exact test for biased strand distribution p-value = 7.94e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACAC  >  NC_000913/298811‑299107
                                                                                                                                                        |                                                                                                                                                   
gTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATAc                                                                                                                                                    >  1:65349/1‑151 (MQ=255)
  tataCACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATAccc                                                                                                                                                  >  2:997417/1‑151 (MQ=255)
   ataCACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATAcccc                                                                                                                                                 <  2:168654/151‑1 (MQ=255)
      cacCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGcc                                                                                                                                              <  1:185508/151‑1 (MQ=255)
       acCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCa                                                                                                                                             <  1:1065122/151‑1 (MQ=255)
          gcgTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCatat                                                                                                                                          <  1:1033843/151‑1 (MQ=255)
               ggCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTg                                                                                                                                      >  2:264812/1‑150 (MQ=255)
                gCGATAGCCGCGCTCCCGCCGCACAGGCCAAGCCCACCGACACTTTTGGCCTTCATCGGGGAGAATAGGGGGTCGTTACCATCCAGAAAAACCCCCCCCTTTTTTGGGATACCCGCATGAACCGGCCCCTCATACCCCGCCATATCGTGAt                                                                                                                                    <  2:1383790/151‑1 (MQ=255)
                gCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGAt                                                                                                                                    >  2:575584/1‑151 (MQ=255)
                 cGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGAtt                                                                                                                                   >  2:530934/1‑151 (MQ=255)
                 cGATAGCCGCGCTCACCCCGCACAGGCCCAGCTCACCGACCCCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGAtt                                                                                                                                   <  1:298977/151‑1 (MQ=255)
                  gATGGCCGCGCTCACGCCGCACAGGCCCAGCTCCCCGACCCCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCACCCAGAAAAATCCCCTCCTTTTTGGGGATATCCGCATGAACCGGCACCTCATCCCCCGCCATATCGTGATTa                                                                                                                                  <  2:1296400/151‑1 (MQ=255)
                  gATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGAtt                                                                                                                                   >  2:191632/1‑150 (MQ=255)
                      gCCGCGCTCACGCCGCACAGGCCCAGCTCCCCAACCCCTTTGGCCTTCATGGGGGAGGATAGGGGGTCGTTATCATCCAGAAAAACCCCCTCCTGTTTTGGAATATCCGCATGAACCGGCCCCTCATCCCCCGCCATATCGTGATTAACGa                                                                                                                              <  2:1330927/151‑1 (MQ=255)
                       ccGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGa                                                                                                                              >  2:1406177/1‑150 (MQ=255)
                         gcgcTCGCGCCGCACGGGCCCACCTCACCGCCCCCTTTGCCCTTCATCGGGGAGGATAGGGGGTCGGTATCACCCAGAAAAACCCCCCCCTGTTTGGGAATACCCGCATAAACCGCCACCTCATCCCCCGCCATATCGTGATTAACGAAGt                                                                                                                           <  2:995188/151‑1 (MQ=255)
                         gcgcTCACGCCGCACAGGCCCAGCTCACCGCCCCCTTTGGCTTTCATGGGGGAGGATATGGGGTCGGTATCATCCAGAAAAATCCCCTCCTGTTTTGGGATATCCGCATGACCCGGCACCTCATCCCCCGCCATATCGTGATTAACGAAGt                                                                                                                           <  1:1195051/151‑1 (MQ=255)
                               ccgccgcACAGGCCCAGCTCACCGACCCCTTTGGCCTTCATGGGGGAGGATATGGGGTCGGTATCACCCAGGAAAATCCCCTCCTTTTTGGGGATATCCGCATGAACCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCa                                                                                                                     <  1:1227648/150‑1 (MQ=255)
                                  ccgaACAGGCCCAGCCCACCGACCCTTTTGGCCTTCTTCGGGAAGAATAGGGGGTCGATATCATCCAGGAAAATCCCCTCCTTTTTTGGAATCTCCGAAAGAACCGGCCCCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAAc                                                                                                                  <  2:1164510/151‑1 (MQ=255)
                                      agaGGCCCAGCTCCCCGACCCCTTTGGCCTTCATCGGGGAGGATATGGGGTGGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATACCCGCATGAACCGGCACCTTATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTc                                                                                                              <  1:1215293/149‑1 (MQ=255)
                                      acaGGCCCAGCTCCCCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGAAAAATCCCCTCCTGTTTTGGGATATCCGCATGAACCGGCCCCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTc                                                                                                              <  1:1263574/151‑1 (MQ=255)
                                         ggCCCAGCTCACCGACACCTTTCGCCTTCATCGGGTAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCTCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATc                                                                                                           >  2:586779/1‑151 (MQ=255)
                                          gCCCAGCTCACCGACCCCTTTGGCCTTCATGGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCCCCTGTTTGGGGATATCCGCATGACCCGGCCCCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATcc                                                                                                          <  1:1088370/151‑1 (MQ=255)
                                           cccAGCTCACCGACCCTTTTGCCCTTCAGCGGGGAGGATATGGGGTCGGTATCACCCAGAAAAACCCCCTCCTTTTTGGGGATACCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCa                                                                                                         <  2:1120656/151‑1 (MQ=255)
                                            ccAGCTCACCGACCCCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCACCCAGAAAAACCCCCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCAc                                                                                                        <  1:1039440/151‑1 (MQ=255)
                                              aGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCg                                                                                                      >  2:1172748/1‑151 (MQ=255)
                                                 cccccaacccttttGCCCTTCAGGGGGGAGAATATGGGTTCTATATCACCCAGAAAATTCCCCCCCTTTTTTGGAATACCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCa                                                                                                   <  2:1050668/140‑1 (MQ=255)
                                                  cACCGACACCTTGGGCCTTCAGGGGGGAGGATAGGGGGGCGGTACCACCCAGAAAAATCCCCTCCGTTTTGGGGATACCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAg                                                                                                  <  2:237593/151‑1 (MQ=255)
                                                    ccGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGct                                                                                                >  2:1133087/1‑151 (MQ=255)
                                                      gACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcc                                                                                              >  2:1434267/1‑151 (MQ=255)
                                                         cccTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGAAAAATCCCCTCCTGTTTTGGGATATCCGCATGAACCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCtcctcc                                                                                           <  1:150610/150‑1 (MQ=255)
                                                          ccTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCtcctcc                                                                                           >  2:454117/1‑150 (MQ=255)
                                                          ccTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGCCACATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCtcctcc                                                                                           >  1:269839/1‑150 (MQ=255)
                                                            tttGCCTTTCATCGGGGAGGATATGGGGTCGGTATCACCCAGAAAAACCCCCCCCTTTTTGGGAATACCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCa                                                                                          <  2:1157604/149‑1 (MQ=255)
                                                             ttGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATc                                                                                        >  2:1413710/1‑150 (MQ=255)
                                                              gggCCTTCTGCGGGGAGGATATGGGTGCGGTATCACCCAGGAAAATCCCCTCCTGTTTGGGGATACCCGCATGAACCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCa                                                                                       <  2:142908/149‑1 (MQ=255)
                                                                     caggggggAGGATAGGGGGTCGGTATCACCCAGAAAAACCCCCCCCTGTTTTGGGATACCCGAATGACCCGGCACCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCcgc                                                                                <  2:233402/146‑1 (MQ=255)
                                                                     cATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCcgcg                                                                               >  1:1012435/1‑151 (MQ=255)
                                                                          gggagtcattGGGGACGGTATAACCCGGGAAACCCCCCCCCGGTTTTGGATAGCCCTAATGAACCGGCCCCTCAGGCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccat                                                                          <  2:1175868/142‑1 (MQ=255)
                                                                          gggAGGATATGGGTTGGAAATCCCCCAGAAAACTCCCCTCCTGTTTGGGCATCTCCACATAAACCGCCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccat                                                                          <  2:639232/151‑1 (MQ=255)
                                                                          gggAGGATATGGGGTCGGTATCATCCAGAAAAACCCCCTCCTGTTTTGGGATATCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccat                                                                          <  1:142660/151‑1 (MQ=255)
                                                                           ggcggatggggtgtgggTACCACCCGGGAAACTCCCCCCCCGTTTGGGTATAGCCGAACGACCGGGCCCCTCACACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatg                                                                         <  1:106900/136‑1 (MQ=255)
                                                                              ggaTAGGGGGTCGGTATCACCCAGAAAAACCCCCTCCTGTTTTGGGATACCCGCAGAACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgccc                                                                      <  1:717853/151‑1 (MQ=255)
                                                                               gaTATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgccc                                                                      >  1:1414554/1‑150 (MQ=255)
                                                                                atatGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgcccat                                                                    >  2:1209039/1‑151 (MQ=255)
                                                                                  agggtgTCGGTATCATCGAGAAAAATCCCCCCCTTTTTTGGTATACCCGTATAACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATa                                                                   <  2:862747/145‑1 (MQ=255)
                                                                                   tGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGt                                                                 >  1:544292/1‑151 (MQ=255)
                                                                                    ggggTCGGTATCACCCAGAAAAATCCCCCCCTTTTTTGGGATACCCGCATGACCCGGCCCCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTc                                                                <  2:811901/151‑1 (MQ=255)
                                                                                    ggggTCGGTATCACCCAGAAAAACCCCCCCCTGTTTTGGGATACCCGAATGAACCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTc                                                                <  2:1462430/151‑1 (MQ=255)
                                                                                     gggTCGTTATCACCCAGAAAAATCACCTCCTTTTTTGGGATATCCGCAGGACCGGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCa                                                               <  1:925204/151‑1 (MQ=255)
                                                                                     gggCCGTTATCACCAAGGAAAACCCCCCCCTTTTTGGGAATACCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCa                                                               <  2:1045682/151‑1 (MQ=255)
                                                                                       gTCGGTATCATCCAGAAAAACCCCCTCCTTTTTTGGGATACCCGCATGAACCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCAtt                                                             <  1:884473/151‑1 (MQ=255)
                                                                                         cGGTATCATCCAGAAAAATCCCCTCCTGTTTTGGGATACCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTg                                                            <  1:202534/150‑1 (MQ=255)
                                                                                         cGGTATCACCCAGGAAAACCCCCCCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTgc                                                           <  1:575586/151‑1 (MQ=255)
                                                                                             atcatcCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGcc                                                        >  1:1307807/1‑150 (MQ=255)
                                                                                               catcCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAAt                                                     >  2:414145/1‑151 (MQ=255)
                                                                                               cacccAGAAAAACCCCCTCCTTTTTGGGGATACCCGCATAACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAAt                                                     <  1:695100/148‑1 (MQ=255)
                                                                                                  ccAGAAAAATCCCCTCCTGTTTTGGGATATCCGCATGACCCGGCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGa                                                   <  1:997424/150‑1 (MQ=255)
                                                                                                  ccAGAAAAATCACCTCCTGTTTGGGGATACCCGCATAACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGAc                                                  <  2:500812/151‑1 (MQ=255)
                                                                                                    agaaaaaTCCCCTCGTTTTTGGGGATATCCGATTAAACGGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCt                                                <  2:885935/148‑1 (MQ=255)
                                                                                                     gaaaaaTCCCCTCCTGTTTTGGGATATCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTg                                               <  1:125212/149‑1 (MQ=255)
                                                                                                        aaaTCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCt                                            >  1:388822/1‑151 (MQ=255)
                                                                                                         aaTCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTg                                           >  1:520079/1‑151 (MQ=255)
                                                                                                            cccctcctTTTTTGGGATATCCGCAGGAACCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgc                                        <  1:1451366/149‑1 (MQ=255)
                                                                                                             aCCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgcg                                       >  2:756843/1‑151 (MQ=255)
                                                                                                               cccctTTTTTGGAATACCCCAATAACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCa                                     <  2:1207344/149‑1 (MQ=255)
                                                                                                                 cctGTTTTGGGATATCCGCATGAACCGGCCCCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGt                                   <  1:129098/151‑1 (MQ=255)
                                                                                                                   tttttcGGGATTTCCGCTTAACCCGGCCCTCCTCCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGtt                                  <  1:157594/148‑1 (MQ=255)
                                                                                                                     tttgGGGATACCCGCATGACCCGCCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGtttt                                <  2:368302/150‑1 (MQ=255)
                                                                                                                      tttGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCg                              >  2:871678/1‑151 (MQ=255)
                                                                                                                      tttGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCg                              >  2:1361533/1‑151 (MQ=255)
                                                                                                                       ttGGGAATACCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTc                               <  2:308670/150‑1 (MQ=255)
                                                                                                                        tGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCgg                             >  2:577042/1‑150 (MQ=255)
                                                                                                                           gATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTc                         >  2:1354714/1‑151 (MQ=255)
                                                                                                                                atcgccaggAACCGGCGCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAgg                      <  1:1304878/143‑1 (MQ=255)
                                                                                                                                 cccccATGACCCGGCACCTAATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAgg                      <  1:264808/146‑1 (MQ=255)
                                                                                                                                  ccGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGAt                    <  2:532321/151‑1 (MQ=255)
                                                                                                                                   cGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATg                   >  2:1418708/1‑151 (MQ=255)
                                                                                                                                    gcaacaaCCGCCACCACATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGc                  <  1:123669/146‑1 (MQ=255)
                                                                                                                                    gCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGc                  <  1:1406172/151‑1 (MQ=255)
                                                                                                                                        gAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTcc              >  1:1020115/1‑151 (MQ=255)
                                                                                                                                         atccggcccttattACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCt             <  1:43688/138‑1 (MQ=255)
                                                                                                                                           ccGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgc           <  1:667469/151‑1 (MQ=255)
                                                                                                                                            cGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgc           >  2:1368391/1‑150 (MQ=255)
                                                                                                                                             gccACCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgcag         <  2:858537/149‑1 (MQ=255)
                                                                                                                                              cccccTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgcagc        <  2:1276049/148‑1 (MQ=255)
                                                                                                                                              ccACCTAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgcagc        <  2:568498/150‑1 (MQ=255)
                                                                                                                                                aCCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcac      >  1:140213/1‑151 (MQ=255)
                                                                                                                                                     aTACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcacacac  >  1:1392370/1‑150 (MQ=255)
                                                                                                                                                        |                                                                                                                                                   
GTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACAC  >  NC_000913/298811‑299107

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