Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,903,7150GA47.4% 36.9 / 84.8 78intergenic (+341/+1032)ygcE/queEputative kinase/7‑carboxy‑7‑deazaguanine synthase; queosine biosynthesis
Reads supporting (aligned to +/- strand):  ref base G (5/36);  new base A (37/0);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 3.58e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTAC‑TTAA‑CCAATGAATAACAAGCTTAAAAAATATACTGTTA‑TTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAACGATTTAAGAATCATACAAATA  >  NC_000913/2903567‑2903861
                                                                                                                                                             |                                                                                                                                                      
cTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGAc                                                                                                                                                      <  2:1082916/151‑1 (MQ=255)
      gTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTT‑TTTTAAAAAA‑C‑A‑CATTT‑AACT‑Ca‑aaaa‑aaaaaaaaaaat                                                                                                                                                >  1:303323/1‑148 (MQ=255)
        cTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAAAAA‑C‑T‑CAATT‑ACCT‑Ca‑aaaa‑aaaaaaaaaaatc                                                                                                                                               >  1:592224/1‑146 (MQ=255)
              aaTAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑t                                                                                                                                         <  2:944010/149‑1 (MQ=255)
              aaTAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑t                                                                                                                                         <  2:825706/149‑1 (MQ=255)
               aTAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTT‑TTTTAAACAA‑C‑T‑CATTA‑TTCT‑CC‑CAAC‑AAAAAAAAACCTcccacattt                                                                                                                                       >  1:1170243/1‑144 (MQ=255)
                  aaaGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑Tttatta                                                                                                                                   <  1:334758/151‑1 (MQ=255)
                  aaaGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑Tttatta                                                                                                                                   <  1:1758014/151‑1 (MQ=255)
                  aaaGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑Tttatta                                                                                                                                   <  1:448560/151‑1 (MQ=255)
                  aaaGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑Tttatta                                                                                                                                   <  2:1107258/151‑1 (MQ=255)
                   aaGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTAc                                                                                                                                  <  1:1499733/151‑1 (MQ=255)
                   aaGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTAc                                                                                                                                  <  2:774653/151‑1 (MQ=255)
                    aGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TTTAAAAACA‑C‑C‑CATTA‑GCCT‑CC‑AAAA‑AAAAAAAAAAATCCAGA‑Tttattttt                                                                                                                                 >  2:1252358/1‑148 (MQ=255)
                      aCAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTAccac                                                                                                                               <  2:272160/151‑4 (MQ=255)
                      aCAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTAC‑tt                                                                                                                               <  2:1661964/150‑1 (MQ=255)
                      aCAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTAC‑tt                                                                                                                               <  2:254099/150‑1 (MQ=255)
                       cAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTT‑TTTAAAACAA‑C‑C‑CTTTA‑‑CCT‑CC‑AAAC‑AAAAAAAAAAATTCAGA‑ATTTTTAA‑ttta                                                                                                                             >  1:1746807/1‑149 (MQ=255)
                       cAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTT‑TTTAAAACAA‑C‑C‑CATTA‑‑CCT‑CA‑AAAC‑AAAAAAAAAAATTCAAA‑TTTTTTTC‑ttta                                                                                                                             >  1:72985/1‑149 (MQ=255)
                         aTACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTT‑TTTTAAAACA‑C‑C‑CCTTT‑ACCT‑CC‑AAAA‑AAAAAAAAAATTCCCGA‑TTTTTTTC‑TTTA‑c                                                                                                                           >  2:712160/1‑150 (MQ=255)
                         aTACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTT‑TGTAAAAAAC‑C‑C‑CATTA‑GTCT‑CC‑AAAA‑AAAAAAAAAATTCCAGA‑TTTTTTTA‑TTaaaa                                                                                                                           >  2:693899/1‑149 (MQ=255)
                         aTACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTAC‑TTAA‑cc                                                                                                                          <  2:187047/151‑1 (MQ=255)
                         aTACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTAC‑TTAA‑cc                                                                                                                          <  2:790150/151‑1 (MQ=255)
                         aTACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAAAAA‑A‑T‑CATTA‑AACC‑CA‑AAAA‑AAAAAAAAAAATCCCGA‑ATTTTTAA‑TTAA‑c                                                                                                                           >  1:1034897/1‑150 (MQ=255)
                         aTACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAACGTATAGAGGCGG‑TTA‑TGTAAAGCCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTAC‑TTAA‑cc                                                                                                                          <  2:6438/151‑1 (MQ=255)
                             gAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTT‑TTTAAAAACA‑C‑A‑CATTA‑ACCT‑CA‑AAAC‑AAAAACAAAAATCCCGA‑TTTTTTAC‑TTAA‑CCAAta                                                                                                                      >  1:496908/1‑150 (MQ=255)
                                   cACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTAC‑TTAA‑CCAATGAATAAc                                                                                                                <  2:230763/151‑1 (MQ=255)
                                     cTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTT‑TTTTAAAACA‑C‑T‑CATTA‑TACT‑CC‑AAAA‑AAAAAAAAAAATTCACA‑TTTTTTAC‑TTTA‑ACAAATAATAaaaa                                                                                                              >  1:1372005/1‑151 (MQ=255)
                                       ttATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTAC‑TTAA‑CCAATGAATAACAAGc                                                                                                            <  1:394669/151‑1 (MQ=255)
                                                   aCTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTTAAAAAA‑A‑C‑CACTA‑CACT‑CA‑AAAA‑AAAAACAAAATTTCAGA‑ATTATTAA‑TTTA‑ACAAAGAAAAACAAACCTAAAAAAaata                                                                                                >  1:1717537/1‑151 (MQ=255)
                                                    cTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTT‑TTTTAAAAAA‑A‑A‑CATTATACCC‑CA‑AAAA‑AAAAAAAAATTTCC‑GA‑TTTTTTAA‑TTTA‑ACAATGAAAAAAAAAACTAAAAAAaataa                                                                                               >  1:321897/1‑150 (MQ=255)
                                                    cTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTAC‑TTAA‑CCAATGAATAACAAGCTTAAAAAATATAc                                                                                               <  1:1286676/151‑1 (MQ=255)
                                                       aaaGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTTC‑TTAA‑CCAATGAATAACAAGCTTAAAAAATATACTGt                                                                                            <  1:1265239/151‑1 (MQ=255)
                                                          gTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTT‑TGTAAAAACA‑C‑T‑CATTA‑ACCT‑CA‑AAAC‑AAAAAAAAAAATCCAGA‑TTTATTAC‑TTAA‑CCAAAGAAAAAAAAGCTTAAAAAATATACTGTTA‑t                                                                                        >  1:81604/1‑151 (MQ=255)
                                                          gTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TTTAAAAAAA‑C‑C‑CATTA‑ACCT‑CC‑AAAA‑AAAAAAAAAATTCCACA‑TTTTTTTCTTTAA‑CAAATAAAAAAAAAACTTaaaaaaaaaatttttt                                                                                          >  2:446428/1‑141 (MQ=255)
                                                            aTGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTAC‑TTAA‑CCAATGAATAACAAGCTTAAAAAATATACTGTTA‑TTc                                                                                      <  2:114832/151‑1 (MQ=255)
                                                              gAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAAATGCCAGA‑TTTATTAC‑TTAA‑CCAATGAATAACAAGCTTAAAAAATATACTGTTA‑TTc                                                                                      <  2:1480015/149‑1 (MQ=255)
                                                                   ggATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑GTT‑TTTTAAAACA‑C‑TCCATTT‑CCCC‑CA‑AAA‑‑AAAAACAAACATCCCGA‑TTTATTTA‑TTAA‑CCAAATAAAAAAAAACTTAAAAAAAATTATGTTA‑tttttaaaaa                                                                               >  1:468651/1‑143 (MQ=255)
                                                                        tCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTT‑TTTTAAAACA‑CAC‑CATTG‑ACTC‑CA‑AAAA‑AAAAAAAAAATTCC‑GA‑ATTTTTAA‑TTTA‑ACAAAAAATAAAAAACATTAAAAAAATAATGTTT‑TTTTTAaaaaaaaaa                                                                          >  1:861817/1‑150 (MQ=255)
                                                                        tCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTT‑TTTAAAAACA‑C‑C‑CATTT‑ACCTCCA‑AAAA‑AAAAAAAAATTTCCAAT‑TTTTTTAT‑TTAA‑CAAATAAAAAAAAAACTAAAAAAAAAAAATTTTT‑TTTTaaaaaaaaaa                                                                           >  2:738791/1‑151 (MQ=255)
                                                                        tCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTT‑TGTTAAAAAA‑C‑A‑CCTTA‑ACCC‑CA‑AAAA‑AAAAACAAACATCCAGA‑TTTTTTAC‑TTAA‑CCAACAAAAAAAAAACTTAAAAAAAATACTATTT‑TTTTTAaaaaaaaaa                                                                          >  1:1144150/1‑150 (MQ=255)
                                                                        tCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTT‑TGTAAAAACA‑C‑C‑CCTTT‑TCCC‑CC‑AAAA‑AAAAACCAAATTCCCGA‑TTTTTTTT‑TTTA‑ACAAATAAAAAAAAACCTAAAAAAAAAAATTTTT‑TTTTaaaaaaaaaaa                                                                          >  2:17448/1‑150 (MQ=255)
                                                                        tCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAAACC‑C‑A‑CCTTA‑ACCC‑CC‑AAAA‑AAAAAAAAACATCAAGA‑ATTTTTAA‑TTTA‑ACAACGACGAACAAAGTTTAAAAATAAACTGTTT‑TTTTATAGaaaaaaa                                                                          >  2:1082793/1‑150 (MQ=255)
                                                                            tATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTAC‑TTAA‑CCAATGAATAACAAGCTTAAAAAATATACTGTTA‑TTCTAAAGAAAAAATTAAc                                                                      <  1:868652/151‑1 (MQ=255)
                                                                                cATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTT‑TTTAAAAACA‑C‑C‑CATTT‑TTCT‑CA‑AAAA‑AAAAAAAAAAATCCAGA‑TTTATTAC‑TTTA‑ACAAAGAAAAAAAAACCTTAAAAAAATAATTTTT‑TTTTTAAAAAAAAAATaaaaaga                                                                  >  2:640727/1‑146 (MQ=255)
                                                                                cATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTT‑TGTTAAAAAC‑C‑T‑CATTT‑GCCT‑CC‑AAAC‑AAAAAAAAAACTCCACATTTTTTTTT‑TTAT‑CAAATTAAAAAAAACCTTAAAAAA‑AAAATCTTT‑TTTTTAAAAAAAAATTaaaaaga                                                                  >  2:237965/1‑146 (MQ=255)
                                                                                cATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTT‑TGTAAAAACA‑C‑A‑CAATA‑TCCT‑CA‑AAAC‑AAAAAAAAAAATTCAGA‑TTTATTTC‑TTTA‑ACAAAGAAAAAAAAAATTTAAAAAAATAATGTTA‑TTTTTTAAAAAAAAATaaaaaga                                                                  >  1:315582/1‑146 (MQ=255)
                                                                                cATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TTTAAAAACA‑C‑C‑CATTT‑TCCT‑CC‑AAAA‑AAAAAAAAAATTCCAGA‑TTTTTTTA‑TTTA‑ACAAAGAATAAAAAACTTAAAAAAAATAATTTTA‑TTTTTAAAAAAAAAATaaaaag                                                                   >  2:209326/1‑146 (MQ=255)
                                                                                cATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TTTAAAAACA‑C‑C‑CATTA‑GCCT‑CA‑AAAC‑AAAAACAAACATCCCGA‑TTTTTTAT‑TTTA‑ACAATAAAAAAAAAACTTAAAAAAAAAAATTTTT‑TTTTTAaaaaaaaaaaaaaaaaa                                                                  >  2:1626856/1‑142 (MQ=255)
                                                                                cATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑GTATTTTAAAAAAA‑C‑C‑ATTTT‑ACCC‑CA‑AAAA‑AAAACAAAAATTCCAAA‑TTT‑TTAA‑TTTA‑CCAAAGAAAAAAAAAAATTAAAAAAATAATTTTATTTTTTAAAAAAAAATTAACAgg                                                                   >  1:1405694/1‑151 (MQ=255)
                                                                                  tAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTT‑TGTAAAAACA‑C‑C‑CCTTA‑TCCT‑CA‑AAAA‑AAAAAAAAAAATCCCGA‑TTTTTTAA‑TTAA‑ACCAAGAATAACAAACTTAAAAAAAAAAATGTTA‑TTTTAAAAAAAAAAATaaaaagaaa                                                                >  1:221487/1‑144 (MQ=255)
                                                                                        cTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTAC‑TTAA‑CCAATGAATAACAAGCTTAAAAAATATACTGTTA‑TTCTAAAGAAAAAATTAACAGGCATTAATTc                                                          <  2:1641956/151‑1 (MQ=255)
                                                                                        cTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTAC‑TTAA‑CCAATGAATAACAAGCTTAAAAAATATACTGTTA‑TTCTAAAGAAAAAATTAACAGGCATTAATTc                                                          <  1:1295084/151‑1 (MQ=255)
                                                                                        cTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTAC‑TTAA‑CCAATGAATAACAAGCTTAAAAAATATACTGTTA‑TTCTAAAGAAAAAATTAACAGGCATTAATTc                                                          <  1:1267432/151‑1 (MQ=255)
                                                                                        cTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTAC‑TTAA‑CCAATGAATAACAAGCTTAAAAAATATACTGTTA‑TTCTAAAGAAAAAATTAACAGGCATTAATTc                                                          <  1:306594/151‑1 (MQ=255)
                                                                                        cTATTCATCAACAGCGTATAGAGGCGG‑CTT‑TGTAAAACCACC‑C‑CTTTT‑ACTT‑CA‑AAAA‑AAAAAAAAAATTCACAT‑TTTTTTAT‑TTAA‑CAAATAAAAAAAAAAC‑TAAAAAAAATAATGTTA‑TTTTTAAAAAAAAAATAAAAAGGAATAAttt                                                          >  1:899630/1‑150 (MQ=255)
                                                                                          aTTCATCAACAGCGTATAGAGGCGG‑TTT‑TTTAAAAACA‑CAC‑CATTT‑TTTT‑AA‑AAAC‑AAAAAAAAACTACCAGT‑TTTTTTAC‑TTAA‑CAAATAATAAAAAAAATTAAAAAAATAACATTTT‑TTTTTATAAAAAAATAAACAGCCTTTAAttaa                                                         >  1:1533892/1‑149 (MQ=255)
                                                                                          aTTCATCAACAGCGTATAGAGGCGG‑TTT‑TTTAAAAACA‑C‑C‑CATTA‑GCCT‑CA‑AAAA‑AAAAAAAAAATTCCCGA‑TTTATTAT‑TTAA‑ACAATAAAAAAAAAACTTAAAAAAAATAATTTTA‑TTTTAAAAAAAAAAATAAAAAGCATTAATTCaa                                                        >  2:1505977/1‑151 (MQ=255)
                                                                                          aTTCATCAACAGCGTATAGAGGCGG‑TTT‑TGTTAAAAAA‑C‑C‑CCTTT‑TCCT‑CC‑AAACAAAAAAAAAAATTCCAAA‑TTTTTTTT‑TTAA‑CAAAAAAAAAAAAACATTAAAAAA‑ATAACTTTT‑TTTTAAATAAAAAAATAAAAAGCACTaaattaa                                                        >  1:1444783/1‑146 (MQ=255)
                                                                                          aTTCATCAACAGCGTATAGAGGCGG‑TTA‑TTTAAAAACA‑C‑C‑CATTA‑ACCT‑CA‑AAAA‑AAAAAAAAAAATCCCGA‑TTTATTAA‑TTAA‑ACAAAAAAAAAAAAAATTAAAAAAAAAAATTTTT‑TTTTTAAAAAAAAAATAAAAAGAATTaaattaa                                                        >  2:381851/1‑146 (MQ=255)
                                                                                          aTTCATCAACAGCGTATAGAGGCGG‑TTA‑TTTAAAAACA‑A‑C‑CATTA‑ACCT‑CC‑AAAA‑AAAAAAAAAATTCCAGA‑TTTTTTTT‑TTAA‑CCAATGAATAAAAAACTTAAAAAAAATACTGTTT‑TTTTAAAGAAAAAAAAAAAAAGAATTAATTTaa                                                        >  1:1339031/1‑151 (MQ=255)
                                                                                          aTTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAAACA‑A‑C‑CCTTA‑ACAC‑CA‑AAAC‑AAAAAAAAACATCCAGATTTTTTTAT‑TAAA‑CCAATAATAAAAAACCTTAAAAAAAAAAAGGTTT‑TTTTAAA‑AAAAAATTAAAAAGCATTAATTCaa                                                        >  1:882540/1‑151 (MQ=255)
                                                                                             catcaACAGCGTATAGAGGCGG‑TTA‑TTTAAAAACA‑C‑C‑CATTA‑GACTCCA‑AAAA‑AAAAAAAAACTCCCAAT‑TTTTTTTT‑TTAA‑CAAAAAAAAAAAAAATTTAAAAAA‑AAAATGTTT‑TTTTTAAAAAAAAAATAAAAAGGATTaaatttaaat                                                     >  2:105781/1‑143 (MQ=255)
                                                                                             catcaACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTAC‑TTAA‑CCAATGAATAACAAGCTTAAAAAATATACTGTTA‑TTCTAAAGAAAAAATTAACAGGCATTAATTCAatat                                                     <  2:195412/151‑1 (MQ=255)
                                                                                             catcaACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTAC‑TTAA‑CCAATGAATAACAAGCTTAAAAAATATACTGTTA‑TTCTAAAGAAAAAATTAACAGGCATTAATTCAatat                                                     <  2:1734275/151‑1 (MQ=255)
                                                                                                     gCGTATAGAGGCGG‑TTA‑TGTAAAAACA‑C‑C‑CATTA‑ACCC‑CA‑AAAA‑AAAAACAAAAATTCAGA‑TTTATTAC‑TTAA‑CCAAAGAATAAAAAACTTAAAAAAAATACTGTTA‑TTTTTAAGAAAAAAATAACAAGCATTAAATCAAAATTataaaaa                                             >  1:946137/1‑151 (MQ=255)
                                                                                                       gTATAGAGGCGG‑TTA‑TGTAAAACCA‑C‑T‑ATCTA‑ACCT‑CC‑AAAC‑AAAAAAAAAATTCCCGA‑TTTTTTTT‑TTAATCCAATAAAAAAAAACATTAAAAAAAAAATTTTTT‑TTTTAAA‑AAAAAATTAAAAGGAATTAATTTAATAATaaaaaaaaa                                           >  1:1511848/1‑149 (MQ=255)
                                                                                                       gTATAGAGGCGG‑TTA‑TGTAAAAACA‑C‑C‑CATTA‑TCCT‑CACAAAA‑AAAAAAAAAATCCCAAA‑TTTATTAT‑TTAA‑‑CAAAAAATAAAAAAATTTAAAAAAAAATTGTTA‑TTTTTAAAAAAAAAATAAAAAGCATTAATTTAAAATTtaaaaaaaa                                           >  1:814537/1‑149 (MQ=255)
                                                                                                             agGCGG‑TTA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTAC‑TTAA‑CCAATGAATAACAAGCTTAAAAAATATACTGTTA‑TTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAg                                       <  1:573007/149‑1 (MQ=255)
                                                                                                                    ttA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTAC‑TTAA‑CCAATGAATAACAAGCTTAAAAAATATACTGTTA‑TTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAAt                               <  1:840708/151‑1 (MQ=255)
                                                                                                                    ttA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTAC‑TTAA‑CCAATGAATAACAAGCTTAAAAAATATACTGTTA‑TTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAAt                               <  1:738124/151‑1 (MQ=255)
                                                                                                                      a‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTAC‑TTAA‑CCAATGAATAACAAGCTTAAAAAATATACTGTTA‑TTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAAttt                             >  1:1048682/1‑151 (MQ=255)
                                                                                                                      a‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTAC‑TTAA‑CCAATGAATAACAAGCTTAAAAAATATACTGTTA‑TTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAAttt                             >  1:814485/1‑151 (MQ=255)
                                                                                                                           aaaaCCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTAC‑TTAA‑CCAATGAATAACAAGCTTAAAAAATATACTGTTA‑TTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTa                         <  2:1319857/151‑1 (MQ=255)
                                                                                                                            aaaCCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTAC‑TTAA‑CCAATGAATAACAAGCTTAAAAAATATACTGTTA‑TTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTaa                        >  1:397547/1‑151 (MQ=255)
                                                                                                                                   c‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTAC‑TTAA‑CCAATGAATAACAAGCTTAAAAAATATACTGTTA‑TTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAACGAttt                  >  1:1506283/1‑151 (MQ=255)
                                                                                                                                              ccT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTAC‑TTAA‑CCAATGAATAACAAGCTTAAAAAATATACTGTTA‑TTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAACGATTTAAGAAt            <  1:950297/149‑1 (MQ=255)
                                                                                                                                                  cA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTAC‑TTAA‑CCAATGAATAACAAGCTTAAAAAATATACTGTTA‑TTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAACGATTTAAGAATCATAc       <  1:1662282/151‑1 (MQ=255)
                                                                                                                                                     aaaC‑AAAGACAAACTTCCAGA‑TTTATTAC‑TTAA‑CCAATGAATAACAAGCTTAAAAAATATACTGTTA‑TTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAACGATTTAAGAATCATACaa     >  2:858268/1‑151 (MQ=255)
                                                                                                                                                       aC‑AAAGACAAACTTCCAGA‑TTTATTAC‑TTAA‑CCAATGAATAACAAGCTTAAAAAATATACTGTTA‑TTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAACGATTTAAGAATCATACAAAt   <  2:763211/151‑1 (MQ=255)
                                                                                                                                                        c‑AAAGACAAACTTCCAGA‑TTTATTAC‑TTAA‑CCAATGAATAACAAGCTTAAAAAATATACTGTTA‑TTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAACGATTTAAGAATCATACAAATa  <  2:49563/151‑1 (MQ=255)
                                                                                                                                                             |                                                                                                                                                      
CTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTA‑TGTAAAACCA‑C‑T‑CATTA‑GCCT‑CA‑AAAC‑AAAGACAAACTTCCAGA‑TTTATTAC‑TTAA‑CCAATGAATAACAAGCTTAAAAAATATACTGTTA‑TTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAACGATTTAAGAATCATACAAATA  >  NC_000913/2903567‑2903861

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