Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,903,6920CA43.1% 101.7 / 45.7 72intergenic (+318/+1055)ygcE/queEputative kinase/7‑carboxy‑7‑deazaguanine synthase; queosine biosynthesis
Reads supporting (aligned to +/- strand):  ref base C (3/38);  new base A (31/0);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 2.69e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

GTTTTTACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAA  >  NC_000913/2903542‑2903839
                                                                                                                                                        |                                                                                                                                                        
gTTTTTACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAAc                                                                                                                                                          <  1:55586/151‑1 (MQ=255)
      aCAAGGAATATTAACTACACTCGATGTCTTTAGAATAAGAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCACT‑c                                                                                                                                                    <  2:1049463/150‑1 (MQ=255)
      aCAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TTTAAAAACa                                                                                                                                                        >  2:315215/1‑147 (MQ=255)
        aaGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTT‑TTTAAAAACAAA‑CAtt                                                                                                                                                 >  1:1070307/1‑151 (MQ=255)
          ggAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TTTGAAAACACT‑c                                                                                                                                                    >  2:221411/1‑146 (MQ=255)
          ggAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAg                                                                                                                                               <  1:599842/151‑1 (MQ=255)
            aaTATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTT‑TGTTAAAAAACT‑CAtt                                                                                                                                                 >  1:848967/1‑147 (MQ=255)
               atTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTT‑TTTAAAAACAAC‑CATTAAcccc                                                                                                                                           >  2:139306/1‑148 (MQ=255)
                  aaCTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TTTAAAAAAACT‑CAtttaactccaa                                                                                                                                        >  1:1187721/1‑141 (MQ=255)
                   aCTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAAACAAA‑CATTTACCT‑C‑a‑aaa                                                                                                                                    >  1:760163/1‑150 (MQ=255)
                       cacTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAg                                                                                                                               <  1:1170316/151‑1 (MQ=255)
                          tCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TTTTAAAAAAAA‑CACTAACCT‑C‑c‑aaacaaaaaa                                                                                                                             >  1:1212981/1‑147 (MQ=255)
                               gTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACtt                                                                                                                       <  2:989603/151‑1 (MQ=255)
                                     aGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTCTAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGa                                                                                                                  <  2:1060399/150‑1 (MQ=255)
                                        aTAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑Ttta                                                                                                             <  1:139300/151‑1 (MQ=255)
                                           aaaGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTTAAAACACC‑CCTTTATCT‑C‑A‑AAACAAAAAAAAAAttttcgattttttt                                                                                                           >  2:150481/1‑139 (MQ=255)
                                           aaaCACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTT‑TTTTAAAAAACC‑CCTTATTCT‑C‑A‑AAAAAAAAAAAAAACTACAGA‑tttttt                                                                                                           >  1:579555/1‑150 (MQ=255)
                                            caGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTT‑TTTTAAAACAAT‑CAATTTCCTCC‑A‑AAAAAAAAAAAAAATTCCACA‑ttattt                                                                                                           >  1:1050447/2‑146 (MQ=255)
                                            aaGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TTTAAAAACACC‑CATTTTCCC‑CAA‑AAAAAAAAAAAAAATTCCAAA‑ttttt                                                                                                            >  2:779913/1‑147 (MQ=255)
                                              gACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTT‑TGTAAAAACACC‑CATTAAACT‑C‑A‑AAAAAAAAAAAAAATTCCCGA‑Attttttttt                                                                                                       >  2:337586/1‑147 (MQ=255)
                                               aCAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTa                                                                                                      <  2:372090/151‑1 (MQ=255)
                                                    aCGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTAACCaa                                                                                                 <  2:635516/151‑1 (MQ=255)
                                                    aCGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACCTAACCa                                                                                                  <  1:31191/150‑1 (MQ=255)
                                                     cGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑GTA‑TTTAAAAACCCC‑CAATAGCCT‑C‑A‑AAAAAAAAAAAAAAATCCCAA‑TTTTTTACTTAACCaa                                                                                                 >  2:907516/1‑150 (MQ=255)
                                                       aaTCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTAACCAATGa                                                                                              <  2:511621/151‑1 (MQ=255)
                                                       aaTCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTAACCAATGa                                                                                              <  2:597128/151‑1 (MQ=255)
                                                        aTCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGATTT‑TTTAAAAAAACC‑CA‑TAACCT‑C‑C‑AAAAAAAAACAAAAATTCACA‑TTTTTTAATTTACCaa                                                                                                 >  2:1054214/1‑147 (MQ=255)
                                                           gCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTT‑TTTTAAAACAAA‑CAATTTACT‑C‑A‑AAAAAAAAAAAAAAATTCCGA‑Tttatttttttaacaaaaaa                                                                                             >  1:1104550/1‑134 (MQ=255)
                                                           gCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTAGATTAAAAAAACT‑CATTATCCA‑C‑C‑AAAAAAAAAAAATTTCCCAGA‑Ttaatatttgtac                                                                                                    >  1:652997/1‑134 (MQ=255)
                                                                     aaTTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTaaa                                                                                <  2:902986/151‑1 (MQ=255)
                                                                      aTTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTT‑TTTAAAAAAAAT‑CATTTGACT‑C‑A‑AAACAAAAAAAAAATTCCAGA‑TTTTTTTTTTAACAAAAGAATAAAAAAATTaaa                                                                                >  1:818732/1‑150 (MQ=255)
                                                                              ttAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTT‑TTTAAAACACCC‑CATTTAACT‑C‑C‑AAACAAAAACAAACCTCCCGA‑TTTATTATTTTATCCAAGAATAAAAAACTTAAAAAAtata                                                                         >  1:435049/1‑149 (MQ=255)
                                                                                 aaGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGAATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGtt                                                                    <  2:490069/151‑1 (MQ=255)
                                                                                   gTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TTTAAAAACACT‑CATTAACCC‑C‑A‑AAAAAAAAACAAAATTCCAGA‑TTTTTTTATTTACCAAAAAAAAAAAAACTTaaaaaaaaaaatttttt                                                                  >  2:891373/1‑140 (MQ=255)
                                                                                     aTGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTT‑TGTTAAAAAAAC‑CATTATCCT‑C‑ACAAAAAAAAAAAAAATCCCAAA‑TTTTTTATTTAA‑AAATAAATAAAAAAATTaaaaaaaaaaatattttt                                                                 >  1:1183035/1‑138 (MQ=255)
                                                                                     aTGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTc                                                                <  1:821572/151‑1 (MQ=255)
                                                                                         aCTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTaaa                                                            <  2:919783/151‑1 (MQ=255)
                                                                                           tGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTTAAAACACC‑CCTTATCCC‑CTA‑AAAAAAAAAAAAACTT‑CAAC‑TTTTTTTTTTTACCAAAGAATAAAAAAATTAAAAAAAAAAATGTtttt                                                                 >  2:1025060/1‑144 (MQ=255)
                                                                                            ggATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAg                                                           <  1:273870/149‑1 (MQ=255)
                                                                                            ggATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAg                                                           <  1:295981/149‑1 (MQ=255)
                                                                                            ggATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGa                                                          <  2:879707/150‑1 (MQ=255)
                                                                                             gATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑ATA‑TTTAAAAAAACT‑CCCTATCCTCC‑A‑AAAAAAAAAAAAAATTCC‑GA‑TTTTTTATTTTAACAAAGAAAAAAAAAATTAAAAAAAATAATGTtttttttaaaaa                                                         >  2:601923/1‑142 (MQ=255)
                                                                                              atatCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TTTAAAAACACT‑CATTAACCT‑C‑A‑AAAAAAAAACAAAAATCCAAA‑TTTTTTATTTTAACAATAAAAAAAAAAATTAAAAAAAAAACTTTTATTTTTAAGaaa                                                        >  2:524703/1‑150 (MQ=255)
                                                                                                atCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTT‑TTTAAAAAAACA‑CATTTCCTT‑C‑A‑AAA‑AAAAACAAAAATCCACA‑GTTTTTTTTTTACCAAAAAATAAAAAACTTAAAAAAAAAAATGTTATTTTTaaaaaaaa                                                      >  1:1197616/1‑150 (MQ=255)
                                                                                                atCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGaaaaa                                                      <  2:1117644/150‑1 (MQ=255)
                                                                                                   tttATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTT‑TGTTTAAAAACC‑CATTATACC‑C‑A‑CAAAAAAAAAAAAATTCCCGA‑TTTTTTACATTAACAAAAAATAACAAACCTAAAAAAAAAAATTTTTTTTTTAaaaaaaaa                                                     >  1:920812/1‑148 (MQ=255)
                                                                                                    ttATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TTTAAAAACACT‑CCTTATCCC‑C‑C‑AAAAAAAAAAAAAAATCCAGA‑TTTTTTTATTTAACAAAGAAAAAAAAACTTAAAAAAAAAAATTTTTTTTTAAAAAAAAAAATa                                                  >  2:738080/1‑150 (MQ=255)
                                                                                                       ttCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAAc                                                <  2:124028/149‑1 (MQ=255)
                                                                                                        tCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAg                                              <  1:961060/150‑1 (MQ=255)
                                                                                                            aaGAACTATTCATCAACAGCGTATAGAGGCGG‑ATA‑TTTAAAAACACA‑CATTTTTCC‑C‑A‑AAAAAAAAAAAAACAACCTTA‑TTTTTTAATTATCCAAAGAAAAAAAAAACTAAAAAAAATTATTTTATTTTTAaaaaaaaaaaaaaaagca                                           >  1:831288/1‑139 (MQ=255)
                                                                                                             aGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑G‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTa                                        <  1:533765/151‑1 (MQ=255)
                                                                                                              gAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TTTTAAAAAACTCCATTA‑TCT‑C‑A‑AAAAAAAAAAAAAAATTCAGA‑TTTTTTTTTTTAACAAATAAAAAAAAAAATAAAAAAAAAAATGTTTTTTTTAAGAAAAAAATaaa                                                >  1:918183/1‑141 (MQ=255)
                                                                                                              gAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTaa                                       <  1:81907/151‑1 (MQ=255)
                                                                                                               aaCTATTCATCAACAGCGTATAGAGGCGG‑TTT‑TTTAAAAACAAT‑CCCTATCCT‑C‑C‑AAAAAAAAACAAAAATCCAAATTTTTTTATTTAACAAAT‑AAAAAAAAAATTAAAAAAAAAAATGTTTTTTTTAAAAAAAAAATAAAAAGAATTaa                                       >  2:646684/1‑150 (MQ=255)
                                                                                                               aaCTATTCATCAACAGCGTATAGAGGCGG‑TTT‑TGTAAAAACACC‑CATTAGCCC‑C‑C‑AAAAAAAAAAAAAAATCCCAA‑TTTTTTTCTTTAACAAAGAAAAAAAAACTTAAAAAAAATAATTTTTTTTTTAAAAAAAAAATAAAAAGGATTaa                                       >  2:548083/1‑150 (MQ=255)
                                                                                                               aaCTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAAt                                      <  2:238250/151‑1 (MQ=255)
                                                                                                                    ttCATCAACAGCGTATAGAGGCGG‑TTT‑TTTAAAAACACC‑CCTTAAACT‑C‑C‑AAAAAAAAAAAAAAATCCCCA‑TTTTTTTCTTTAACAAAGAATAAAAAAATTTAAAAAAATACTGTTTTTTTAAAAAAAAAAATAAAAAGCATTaa                                       >  2:1049/1‑145 (MQ=255)
                                                                                                                      catcaACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAata                                <  1:192730/150‑1 (MQ=255)
                                                                                                                         caACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaa                            <  2:638627/151‑1 (MQ=255)
                                                                                                                         caACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTa                             <  1:1037498/150‑1 (MQ=255)
                                                                                                                         caACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTa                             <  1:1483/150‑1 (MQ=255)
                                                                                                                         caACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTa                             <  2:859331/150‑1 (MQ=255)
                                                                                                                             aGCGTATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaaaaaa                        <  1:253421/151‑1 (MQ=255)
                                                                                                                             aGCGTATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaaaaaa                        <  1:987429/151‑1 (MQ=255)
                                                                                                                                gTATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACATCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAaata                     <  1:868093/151‑1 (MQ=255)
                                                                                                                                  ataGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAaataa                    <  1:928932/150‑1 (MQ=255)
                                                                                                                                          gg‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATAGTAAAAAAATAATAGATTaaa           <  1:614657/151‑1 (MQ=255)
                                                                                                                                           g‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTaaa           >  2:875738/1‑150 (MQ=255)
                                                                                                                                              tA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAAtt        >  2:175054/1‑151 (MQ=255)
                                                                                                                                                   tAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAAAATCCAGA‑GTTATTACTTAACCAATGAATAACAAGATTAAAAAATAAACTGTAATTATAAAGAAAAAATTAACAGGCAGTAATTCAATATTAAAAAAATAATAGATTAAAATTTCt     <  1:210914/150‑1 (MQ=255)
                                                                                                                                                    aaaaCCACT‑CATTAGCCT‑C‑A‑AAACAAAGATAAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCtt    <  2:667626/150‑1 (MQ=255)
                                                                                                                                                     aaaCCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTaa  <  2:877801/151‑1 (MQ=255)
                                                                                                                                                     aaaCCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTaa  <  2:1021269/151‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                        
GTTTTTACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCACT‑CATTAGCCT‑C‑A‑AAACAAAGACAAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAA  >  NC_000913/2903542‑2903839

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