Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913298,9530AC43.5% 111.2 / 23.7 70V661G (GTG→GGG) paoCPaoABC aldehyde oxidoreductase, Moco‑containing subunit
Reads supporting (aligned to +/- strand):  ref base A (32/7);  new base C (0/30);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 3.31e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCAC  >  NC_000913/298803‑299103
                                                                                                                                                           |                                                                                                                                                      
ccGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCa                                                                                                                                                        <  1:993554/151‑1 (MQ=255)
     ggCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTc                                                                                                                                                    <  2:512726/150‑1 (MQ=255)
       cGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCAt                                                                                                                                                  >  1:1035391/1‑150 (MQ=255)
       cGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATa                                                                                                                                                 <  2:70515/151‑1 (MQ=255)
        gTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATAc                                                                                                                                                >  1:1294600/1‑151 (MQ=255)
            taCACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCg                                                                                                                                            >  1:61483/1‑151 (MQ=255)
             acacCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGc                                                                                                                                           >  1:1366328/1‑151 (MQ=255)
                ccGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCCCCGACCCCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCACCCAGAAAAATCCCCTCCTTTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCat                                                                                                                                        <  2:765384/151‑1 (MQ=255)
                      tGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTg                                                                                                                                  >  2:409893/1‑151 (MQ=255)
                       ggCGATCGCCGCGCTCACGCCGCACGGCCCCAGCTCACCGACCCCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCACCCAGAAAAACCCCCCCCTTTTTGGGGATATCCGCATGACCCGGCACCTCATCCCCCGCCATATCGt                                                                                                                                   <  2:1199028/149‑1 (MQ=255)
                               ccGCCCTCACGCCGAACGGCCCCAGCCCCCCGACCCCTTTGGCCTTCTCCGGGGGGGATTTGGGGCCGGTATCCTCCAGGAAAACCCCCTCGTTTTTTGTGATACCCCCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGaa                                                                                                                         <  1:1239363/151‑1 (MQ=255)
                                cgcgCTCCCGCCGCACAGGCCCAGCCCCCCGACCCTTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGAAAAACCCCCTCCTTTTTGGGGATATCCGCATGAACCGGCCCCTCATCCCCCGCCATATCGTGATTAACg                                                                                                                           <  2:1282879/148‑1 (MQ=255)
                                   gcTTACCCCGCCCCGGCCCAGCCACCCTCCCCTTTTGGCCTTAACGGGGGGGGATATGGGTTCGTTATCAACCAGAAAAATCCCCTCTTTTTTTGGGATACCCGCATGAACCGGCACCTAATCCCCCGCCATATCGTGATTAACGAAGTa                                                                                                                      <  1:210583/150‑1 (MQ=255)
                                      ccctccgcccggccccAGCTCCCCGACCCCTTGGTCCTTCAGCGGGGAGTATTGGGGTAGGGTATAATCCAGGAAATTCCCCTCCTGTTTGGGGATTTCCCCGTAACCCGGCTCCTCACCCCCCGCCATATCGTGATTAAc                                                                                                                            <  1:800100/128‑1 (MQ=255)
                                      cACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGccc                                                                                                                  >  1:1335656/1‑151 (MQ=255)
                                        cgccgcACGGGCCCAGCTCACCGACCCTTTTGGCCTTCAGGGGGGAGGATAGGGGTTCGGTATCATCCAGAAAAACCCCCCCCTGTTTGGGGATACCCGAATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCaa                                                                                                                <  2:336451/151‑1 (MQ=255)
                                        cgccgcACAGGCCCAGCTCCCCGACCCCTTTGGCCTTCATCGGGGAGGATAGGGGGTCGGTATCATCCAGGAAAATCACCCCCTGTTTTGGGATATCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCa                                                                                                                 <  1:458696/150‑1 (MQ=255)
                                        cgccgcACAGGCCCAGCCCCCCGACCCCTTTGGCCTTCACGGGGGAGGATATGGGTTCGGTATCACCCAGAAAAATCACCCCCTTTTTGGGAATACCCGCATAAACCGGCCCCTAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCa                                                                                                                 <  2:362929/150‑1 (MQ=255)
                                        ccccgcACAGGCCCAGCTCCCCGACACCTTTGCCTTTCATCGGGGGGGATAGGGGGTCGGTATCCCCCAGAAAAACCCCCCCCTTTTTGGGGATACCCGCATCACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGccct                                                                                                                 <  2:226225/148‑2 (MQ=255)
                                          ccgcACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAAc                                                                                                              >  1:302495/1‑151 (MQ=255)
                                          ccgcACAGGCCCAGCTCACCAACCCTTTTGGCTTTCATCGGGGAGGATATGGGGTCGTTATCACCCAGAAAAATCCCCTCTTTTTTTGGGATACCCGCATGACCCGCCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAAc                                                                                                              <  2:1212741/151‑1 (MQ=255)
                                           cgcACAGGCCCAGCCCCCCGACCCCTTTGGCCTTCATCGGGGAGGATAGGGGGTCGGTATCATCCAGAAAAACCCCCCCCTGTTTTGGGATATCCGCATGACCCGGCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACg                                                                                                             <  1:904326/151‑1 (MQ=255)
                                             cacaGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGt                                                                                                           >  1:1408878/1‑151 (MQ=255)
                                               cggccccAGCCCCCCGACACCTTTGGCTTTCATGGGGGAGGATAGGGGGTCGGTATCATCCAGAAAAACCCCCTCCTGTTTTGGGATACCCGCATGAACCGCCCCCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCaaaag                                                                                                             <  1:402819/143‑3 (MQ=255)
                                                aGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCAt                                                                                                        >  2:705650/1‑151 (MQ=255)
                                                 ggCCCAGCTCACCGACCCCTTTGGCCTTCAGCGGGGAGGATAGGGGGTCGGTATCACCCAGAAAAACCCCCTCCTGTTTTGGGATACCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTc                                                                                                          <  2:649977/148‑1 (MQ=255)
                                                 ggCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATc                                                                                                       >  2:996069/1‑151 (MQ=255)
                                                  gCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATcc                                                                                                      >  2:1208990/1‑151 (MQ=255)
                                                    ccAGCCCACCGCCGCCTTTGGCCTTCAGGGGGGGGAATGGGGGTTTGGTTTCTCCAAGAAAAATCCCTCCCTTTTTTGGGATACCCCTATGAACCGGCCTCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATcc                                                                                                      <  1:15302/149‑1 (MQ=255)
                                                        ccccccGACCCCTTTGGCCTTCATCGGGGAGGATATGGGGTCGTTATCATCAAGAAAAACCCCCTCCTTTTTTGGAATACCCGCATGAACCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCg                                                                                                  <  2:242284/145‑1 (MQ=255)
                                                          cccccgccccTTTGGGCCTTCATGGGGGAGGATGTGGGGTGGTTATCACCCAGAAAAATCCCCCCTTTTTTTGGAATACCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCa                                                                                               <  1:1358314/143‑1 (MQ=255)
                                                                  ccTTTGCCTTTCATGGGGGAGGATTGGGGGTCGTTATCACCCAGAAAAACCCCCTCCTTTTTGGGGATACCCGCATGAACCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcct                                                                                         <  1:416780/148‑1 (MQ=255)
                                                                      ttGCCCTTAATGGGGGGGAATGGGGGGTCGTTATCACCAAGAAAAACCCCCTCTTTTTTGGGGATTCCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATc                                                                                    <  2:1423891/150‑1 (MQ=255)
                                                                        tgcccTTCATGGGGGAGAATATGGGGTGGGTATCACCCAGAAAAATCCCCTCCTTTTTTGGGATACCCGCAGAACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCa                                                                                   <  1:77317/147‑1 (MQ=255)
                                                                           ccTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGc                                                                               >  2:1218913/1‑151 (MQ=255)
                                                                             ttCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCcg                                                                             >  1:1303427/1‑151 (MQ=255)
                                                                               caggggggAGGATAGGGGGTCGGTACCACCAAGAAAAACCCCCCCTTTTTTGGGGATACCCGCATGACCCGGCCCCTAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGcc                                                                              <  2:1391577/144‑1 (MQ=255)
                                                                                   ggggAGGATATGGGGTCGTTATCACCCAGAAAAACCCCCTCCTGTTTTGGGATATCCGCATGACCCGGCCCCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgccca                                                                       <  1:1236092/151‑1 (MQ=255)
                                                                                    gggAGGATAGGGGGTCGTTATCACCCAGAAAAACCCCCTCTTTTTTTGGAATACCCGCATAACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccat                                                                      <  1:99869/151‑1 (MQ=255)
                                                                                       aggaTATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgcc                                                                   >  2:845458/1‑151 (MQ=255)
                                                                                           tatGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgcccat                                                                >  2:875071/1‑150 (MQ=255)
                                                                                              ggggTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTc                                                            >  1:362932/1‑151 (MQ=255)
                                                                                              ggggCCGGTACCATCAAGAAAACTCCCTCCTTTTTTGGGAATCCCCCCAACACCCGCCCCCTTTCCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTc                                                            <  1:46964/151‑1 (MQ=255)
                                                                                                 gTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCAtt                                                         >  1:1217646/1‑151 (MQ=255)
                                                                                                     ggTATCACCCAGGAAAATCCCCTCCTTTTTTGGAATACCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTgc                                                       <  1:347614/150‑1 (MQ=255)
                                                                                                        atcatcCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGcc                                                    >  1:1349403/1‑150 (MQ=255)
                                                                                                          catcCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAAt                                                 >  1:666627/1‑151 (MQ=255)
                                                                                                             ccAGAAAAATCCCCTCCTTTTTGGGAATACCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATg                                                <  1:1192428/149‑1 (MQ=255)
                                                                                                               aGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCt                                            >  1:716873/1‑151 (MQ=255)
                                                                                                                gaaaaaTCCCCTCTTTTTTTGGGATCCCCGCAGAACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATg                                                <  1:416711/144‑1 (MQ=255)
                                                                                                                  aaacccccctcttTTTTTGGGATACCCGCATAACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGa                                               <  2:1274067/138‑1 (MQ=255)
                                                                                                                  aaaaTCCCCTCTTGTTTTGGAATACCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGc                                         <  1:409950/151‑1 (MQ=255)
                                                                                                                   aaacccccccccgTTTTGGGAATACCCGCAGAACCCGGCCCCTAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATg                                                <  2:1469600/132‑1 (MQ=255)
                                                                                                                        aCCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgc                                    >  2:545054/1‑150 (MQ=255)
                                                                                                                         cctcttTTTTGGGGATACCCGCATGACCCGCCCCCTAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgtgc                                    <  1:798527/149‑1 (MQ=255)
                                                                                                                            cctGTTTGGGGATACCCCCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCa                                 <  1:830755/149‑1 (MQ=255)
                                                                                                                              tGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGttt                             >  2:1152811/1‑151 (MQ=255)
                                                                                                                              tGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGttt                             >  2:72581/1‑151 (MQ=255)
                                                                                                                                ttgtGGGATACCCGAATGACCCGCCCCCTAACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGtt                              <  2:26017/145‑1 (MQ=255)
                                                                                                                                 tttGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCg                          >  2:270525/1‑151 (MQ=255)
                                                                                                                                  ttGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCgg                         >  2:1082861/1‑151 (MQ=255)
                                                                                                                                  ttGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCgg                         >  1:385817/1‑151 (MQ=255)
                                                                                                                                   tGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCgg                         >  2:646558/1‑150 (MQ=255)
                                                                                                                                      ggaataccccaaTGACCCGCCCCCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCg                          <  2:1366225/137‑1 (MQ=255)
                                                                                                                                         atatCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTc                     >  1:330757/1‑150 (MQ=255)
                                                                                                                                          gccccgcgggacccGGCCCCAAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTc                     <  2:1035291/137‑1 (MQ=255)
                                                                                                                                            tCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAgg                  >  2:1354167/1‑150 (MQ=255)
                                                                                                                                            tCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAgg                  >  1:810606/1‑150 (MQ=255)
                                                                                                                                             ccGCATGACCCGCCCCCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGAt                <  2:1335549/151‑1 (MQ=255)
                                                                                                                                               gCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGc              >  1:866963/1‑151 (MQ=255)
                                                                                                                                                       cGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgc       >  2:1434655/1‑150 (MQ=255)
                                                                                                                                                           aCCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcac  >  1:92925/1‑151 (MQ=255)
                                                                                                                                                           |                                                                                                                                                      
CCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCAC  >  NC_000913/298803‑299103

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: