Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913298,9600AC50.7% 76.2 / 113.7 70Y659D (TAT→GAT) paoCPaoABC aldehyde oxidoreductase, Moco‑containing subunit
Reads supporting (aligned to +/- strand):  ref base A (28/6);  new base C (0/35);  total (29/41)
Fisher's exact test for biased strand distribution p-value = 8.02e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.00e-11

ATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCAC  >  NC_000913/298814‑299103
                                                                                                                                                       |                                                                                                                                               
ataCACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATAcccc                                                                                                                                             <  2:644099/151‑1 (MQ=255)
  acacCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCg                                                                                                                                            <  1:1345842/150‑1 (MQ=255)
  acacCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGc                                                                                                                                           <  1:377873/151‑1 (MQ=255)
    acCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGcc                                                                                                                                          >  1:1294220/1‑150 (MQ=255)
        cgTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATc                                                                                                                                     >  1:43259/1‑151 (MQ=255)
         gTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCg                                                                                                                                    >  1:402257/1‑151 (MQ=255)
          ttGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCa                                                                                                                                         >  2:1223203/1‑145 (MQ=255)
          ttGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTCTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCa                                                                                                                                         <  1:1222361/145‑1 (MQ=255)
                atggCCGCGCTCACGCCGCACAGGCCCAGCCCACCGACCCCTTTGG‑CCTTCATCGGGGAGGATAGGGGTTGGGTATCACCCAGGAAAATCCCCTCCTGTTT‑GGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATCCCCCGCCATATCGTGAtt                                                                                                                               <  2:79895/146‑1 (MQ=255)
                   gCCGCGCTCCCGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCCCCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CCCCTCATCCCCCGCCATATCGTGATTAAc                                                                                                                            <  2:988219/149‑1 (MQ=255)
                   gCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGCGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACg                                                                                                                           >  1:1069645/1‑150 (MQ=255)
                    ccGCCCTCCCGCCGCACAGCCCCAGCTCCCCGCCCCCTTGGG‑CCTTCATCGGGGAGGATAGGGGGTCGGTATCACCCAGAAAAACCCCCTCCTTTTT‑TGGA‑ATATCCGC‑AGGAACCGC‑CCCCTCATCCCCCGCCATATCGTGATTAACg                                                                                                                           <  1:882702/149‑1 (MQ=255)
                     cccgcaggcccgccacgGGCCCACCTGCCCGACACCTTTGC‑CCTTCATCGGGAAGGTAAGGGGTTCGGTACAAACCGGAAAACTCCCCTCCTTTGG‑GGGG‑TTGCCCGT‑ATGCCCCGC‑CACCTAATACCCCGCCATATCGTGATTAACGAAg                                                                                                                        <  1:1016620/138‑1 (MQ=255)
                         cTCCCGCCGCACGGGCCCACCTCCCGGACCCTTTTGG‑CCTTCACCGGGGAGGATAGGGGGTCGGTTTCATCCAGAAAAATCCCCTCCTGTTT‑TGGG‑ATATCCGA‑AAGAACCGG‑CCCCTAATCCCCCGCCATATCGTGATTAACGAAGTAg                                                                                                                     <  2:585918/150‑1 (MQ=255)
                           cACGCCGCACAGGCCCAGCTCACCGACCCCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCACCCAGAAAAACCCCCCCCTGTTT‑TGGG‑ATATCCGC‑AGGAACCGG‑CACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGccc                                                                                                                  <  2:984854/151‑1 (MQ=255)
                           cACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGccc                                                                                                                  >  1:636096/1‑151 (MQ=255)
                              cccgcATAAGCCCCGGCCACCGACACTTTCGG‑CCTTCGTGGGGGAGGCTGAGGTTGTGGAAACCCCCAGGAAAATCCCCTCCTTTTG‑TGAG‑AT‑TCCGCTAGGTACGGC‑CACAGAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCa                                                                                                                 <  1:19106/148‑1 (MQ=255)
                                cgcACAGGCCCAGCCCCCCGACCCTTTGGG‑CCTTCATCGGGGAGGATAGGGGGTGGTTATCATCCAGAAAAATCCCCTCCTTTTT‑TGGG‑ATATCCGC‑ATGACCCGC‑CCCCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAAc                                                                                                              <  1:1306045/150‑1 (MQ=255)
                                   acaGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTc                                                                                                          >  1:928543/1‑151 (MQ=255)
                                      ggCCCAGCTCCCCGACCCCTTTGG‑CCTTCATCGGGGAGGATAGGGGGTCGGTATCACCCAGAAAAACCCCCTCCTGTTT‑TGGG‑ATACCCGC‑ATGAACCGG‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACgg                                                                                                            <  2:1086778/146‑1 (MQ=255)
                                      ggCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATAGGGGGTCGGTATCATCCAGAAAAATCCCCTCCTTTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTc                                                                                                          <  2:420866/148‑1 (MQ=255)
                                        cccAGCCCCCCGCCCCCTTTGG‑CCTTCATCGGGGAGGATAGGGGGTCGGTATCATCCAGGAAAATCCCCTCCTTTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATcc                                                                                                      <  2:916002/150‑1 (MQ=255)
                                        cccAGCCCCCCGCCACCTTTGG‑CCTTAAGGGGGGGGGATAGGGGGTCGGTATCACCCAGAAAAACCCCCTCCTGTTT‑TGGA‑AAACCCGC‑ATGAACCGG‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCa                                                                                                     <  1:27753/151‑1 (MQ=255)
                                            gCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGc                                                                                                 >  1:751242/1‑151 (MQ=255)
                                               cACCGCCCCTTTTGC‑CCTTCATGGGGGAGGATAGGGGGTCGTTACCACCAAGAAAACTCCCCTCTTTTTT‑GGGT‑ATATCCGC‑ATGAACCGG‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCa                                                                                               <  2:265051/150‑1 (MQ=255)
                                                 ccgcccccTTTGC‑CCTCCATGGGGGAGGATAGGGGGCATGTACCACCCAGGAAAATCCCCTCTTGTTT‑TGGG‑TTATCCGC‑ATGACCGGC‑CCCCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGct                                                                                            <  1:651930/145‑1 (MQ=255)
                                                  cGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctc                                                                                           >  1:121353/1‑151 (MQ=255)
                                                      cccTTTGGCCCTTCAGGGGGGGGGATGGGGGGGCGTTATCACCCA‑GAAAAACCCCCCCTTTTT‑TGGGTATACCCGA‑AGGACCGGG‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcc                                                                                          <  1:73765/148‑1 (MQ=255)
                                                            gg‑CCTTCATGGGGGAGGAATGGGGGTAGATAGCACCCAGAAAAACCCCCTCCTGTTT‑GGGG‑ATACCCGA‑ATGACCCGG‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAg                                                                                  <  1:1269810/150‑1 (MQ=255)
                                                             gcccTTCATGGGGGAGGATAGGGGGTCGTTATCACCCAGAAAAATCCCCCCCTGTTT‑GGGG‑ATACCCGC‑ATGACCCGC‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGt                                                                                 <  2:838740/149‑1 (MQ=255)
                                                             gcccTTCAGCGGGGAGGATAGGGGTTCGTTATCCCCCAGAAAAATCCCCTCCTTTTT‑TGGG‑ATATCCGC‑ATGACCCCC‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCa                                                                                   <  2:477981/147‑1 (MQ=255)
                                                                 ttCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCcg                                                                             >  1:1273360/1‑151 (MQ=255)
                                                                       ggggAGGATAGGGGGTGGTTATCACCAAGGAAAATCCCCTCCTGTTT‑TGGG‑ATACCCGC‑ATAACCCGC‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCcgc                                                                            <  2:409509/146‑1 (MQ=255)
                                                                       ggggAGAATTGGGGGTGGGTACCA‑CCAAGAAAACCCCCCCCTTTTTGTGGT‑ATATCCGG‑TTAACCCGG‑CCCCTAACCCCCCGCCATATCGTGATTAGCGAAGTAGCCCAAGCGGTCATCCTCCGCCAGCTCCTCCATCAGTGc                                                                               <  1:213657/143‑1 (MQ=255)
                                                                       ggggAGAATTGGGGGCCGGTTTCACCCGGAAAAACCCCCCCCTTTTG‑AGGG‑ATA‑CCGC‑AAGAACCGGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCcgc                                                                            <  2:121372/146‑1 (MQ=255)
                                                                        gggAGGATAGGGGGTGGGTATCACCCAGAAAAATCCCCTCTTGTTT‑TGGG‑ATACCCGC‑ATGACCCGG‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccat                                                                      <  2:468769/151‑1 (MQ=255)
                                                                        gggAGGATAGGGGGTCGGTATCACCCAGAAAAATCCCCCCCTGTTT‑TGGG‑ATATCCGC‑ATGACCCGC‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccat                                                                      <  2:1233507/151‑1 (MQ=255)
                                                                                  ggggTGGGTATCACCCAGAAAAATCCCCTCCTTTTT‑TGGG‑ATACCCGC‑ATAACCCGG‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGt                                                             <  2:1320131/150‑1 (MQ=255)
                                                                                  ggggTCGTTATCACCGAGGAAAACCCCCCCTTTTTT‑TGGG‑ATATCCAC‑AGGACCGGC‑CCCCTCAACCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgcccat                                                                <  2:94029/147‑1 (MQ=255)
                                                                                  ggggTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTc                                                            >  1:1368278/1‑151 (MQ=255)
                                                                                         gTACCACCCAGGAAACTCCCCTCCTGTTT‑GGGT‑ATACCCGC‑ATGACCCGG‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTgcgc                                                     <  1:696477/151‑1 (MQ=255)
                                                                                          tATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGcc                                                    >  2:1229410/1‑151 (MQ=255)
                                                                                                ccAGAAAAATCCCCTCTTTTTT‑TGGG‑ATACCCGC‑ATGACCCGG‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTgcgc                                                     <  2:1065382/144‑1 (MQ=255)
                                                                                                ccAGAAAAATCCCCTCCTTTTT‑GGGG‑ATCTCCGC‑ATGACCCGC‑CACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGAc                                              <  2:87355/151‑1 (MQ=255)
                                                                                                  aGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCt                                            >  1:1178795/1‑151 (MQ=255)
                                                                                                     aaaaTCCCCTCCTGTTT‑GGGG‑ATCCCCCC‑ATGACCCGC‑CCCCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATg                                                <  1:102854/144‑1 (MQ=255)
                                                                                                       aaccccctcctGTTT‑TGGG‑ATATCCGA‑ATGACCCGC‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTg                                       <  2:528718/146‑1 (MQ=255)
                                                                                                        accccctcttGTTT‑TGGG‑ATACCCGC‑AGAACCCGG‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTg                                       <  2:585478/146‑1 (MQ=255)
                                                                                                         tCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcg                                     >  1:1065351/1‑151 (MQ=255)
                                                                                                         tCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcg                                     >  2:658472/1‑151 (MQ=255)
                                                                                                          cACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcg                                     >  1:581665/1‑150 (MQ=255)
                                                                                                           aCCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgcg                                   >  1:395882/1‑151 (MQ=255)
                                                                                                             ctcctGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgcgc                                  >  2:726682/1‑150 (MQ=255)
                                                                                                             ctcctGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgcgc                                  >  1:1240883/1‑150 (MQ=255)
                                                                                                             ctcctGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAATCATTGCGCCAATGACCTGGCTgcgcgc                                  >  2:690327/1‑150 (MQ=255)
                                                                                                                 tGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGttt                             >  1:1381215/1‑151 (MQ=255)
                                                                                                                   ttt‑tGTT‑ATCTCCGC‑TTGTACCGG‑CACCTTTCTCCCCGCCATATCGTGCGTGACGAGGTAACCCAGCCGGTCATCTACCGGCAGCTCCTGCAGCAGGGCGGCGGCCTTGCCCATAGTCATTGCGCCAAGGTCCTGGCGGCGCGCGGttt                             >  1:56818/1‑149 (MQ=255)
                                                                                                                    tt‑tGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCg                          >  1:159232/1‑151 (MQ=255)
                                                                                                                       tGGG‑ATATCCGC‑ATAAACCGG‑CACCTCATACCCCGCCATATCGTGCTTAACGAAGTAGACCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCAGGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGGTTTCgg                         >  1:492882/1‑150 (MQ=255)
                                                                                                                             tatCCGC‑ATGAACCGG‑CAACTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCACTGACCTGGCCGCGCGCAGTTTTCGGATTc                     >  1:376573/1‑149 (MQ=255)
                                                                                                                                ccGC‑ATGAACCGG‑CCCCTCATCCCCCGCCATATCGTGATCAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCa                    <  2:392565/147‑1 (MQ=255)
                                                                                                                                ccGC‑ATGAACCGC‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAgg                  <  2:395803/149‑1 (MQ=255)
                                                                                                                                  gC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGc              >  2:1044251/1‑151 (MQ=255)
                                                                                                                                     catcaaccccccccTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATg               <  1:351050/137‑1 (MQ=255)
                                                                                                                                         aCCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTg        >  1:450130/1‑151 (MQ=255)
                                                                                                                                          ccGG‑CCCGTCAACCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCt         <  1:1132032/149‑1 (MQ=255)
                                                                                                                                            cgcccccTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTcc          <  1:397130/142‑1 (MQ=255)
                                                                                                                                            cgcccccTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCt         <  1:719108/143‑1 (MQ=255)
                                                                                                                                                aCCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcac  >  2:757136/1‑151 (MQ=255)
                                                                                                                                                 ccTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcac  >  1:871811/1‑150 (MQ=255)
                                                                                                                                                   tcccaCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcac  <  1:914293/144‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                               
ATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCAC  >  NC_000913/298814‑299103

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