Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913298,9600AC49.3% 92.6 / 33.1 69Y659D (TAT→GAT) paoCPaoABC aldehyde oxidoreductase, Moco‑containing subunit
Reads supporting (aligned to +/- strand):  ref base A (32/3);  new base C (0/34);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.50e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCAC  >  NC_000913/298810‑299103
                                                                                                                                                           |                                                                                                                                               
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:72581/1‑151 (MQ=255)
                                                                                                                       tGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGttt                             >  2:1152811/1‑151 (MQ=255)
                                                                                                                         ttgtGGGATACCCGAATGACCCGCCCCCTAACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGtt                              <  2:26017/145‑1 (MQ=255)
                                                                                                                          tttGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCg                          >  2:270525/1‑151 (MQ=255)
                                                                                                                           ttGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCgg                         >  1:385817/1‑151 (MQ=255)
                                                                                                                           ttGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCgg                         >  2:1082861/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                  >  1:810606/1‑150 (MQ=255)
                                                                                                                                     tCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAgg                  >  2:1354167/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)
                                                                                                                                                    cccTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCTCAGTTTTCGGATTCAGGATGCGTCCTg        <  1:1265862/144‑1 (MQ=255)
                                                                                                                                                    aCCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcac  >  1:92925/1‑151 (MQ=255)
                                                                                                                                                           |                                                                                                                                               
CGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCAC  >  NC_000913/298810‑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: