Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,434,8610GT44.6% 67.4 / 37.3 66P259Q (CCA→CAA) truAtRNA pseudouridine(38‑40) synthase
Reads supporting (aligned to +/- strand):  ref base G (9/27);  new base T (29/0);  total (39/27)
Fisher's exact test for biased strand distribution p-value = 7.54e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGG‑TAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCC‑GTTGCTGCCGCCAGCGTTCTG‑TCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCT  >  NC_000913/2434712‑2435009
                                                                                                                                                     |                                                                                                                                                       
aaaaTCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTgg                                                                                                                                                        <  2:578486/151‑1 (MQ=255)
       aTGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCa                                                                                                                                                 <  1:972493/151‑1 (MQ=255)
         gAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCAt                                                                                                                                                <  1:369901/150‑1 (MQ=255)
            gAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAcc                                                                                                                                             >  2:146690/1‑150 (MQ=255)
                taaatCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAACgg                                                                                                                                           >  2:850474/1‑148 (MQ=255)
                 aaatCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCgg                                                                                                                                           >  2:555662/1‑147 (MQ=255)
                     cAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCAAAACGGGCAggg                                                                                                                                     >  2:875122/1‑149 (MQ=255)
                      aGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGG‑Taa                                                                                                                                 <  2:293900/151‑1 (MQ=255)
                        gTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGG‑TAATc                                                                                                                               <  2:796212/151‑1 (MQ=255)
                          ccATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGCTTTTGGAGAGGATACCTGTGAGGG‑Taaatc                                                                                                                              >  1:397209/1‑147 (MQ=255)
                             tatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTAGGAGATCATAACGGTCA‑GG‑TAAaccacccc                                                                                                                         >  2:670534/1‑149 (MQ=255)
                              atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGG‑TAATCCACcgc                                                                                                                         <  2:376156/151‑1 (MQ=255)
                              atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGG‑TAATCCACcgc                                                                                                                         <  1:146086/151‑1 (MQ=255)
                                  gCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGG‑GAATTCACCCCGGcc                                                                                                                     >  1:928900/1‑151 (MQ=255)
                                    cTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAATATAACGGGCAGGG‑GAATACAcccc                                                                                                                         >  2:530771/1‑143 (MQ=255)
                                     tATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATACTAACCGGCCGGG‑GAATCCAcccccaac                                                                                                                     >  1:500814/1‑142 (MQ=255)
                                     tATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGAACATACCGGGCAGGG‑GAATTCCCCGGGAACa                                                                                                                    >  2:262075/1‑149 (MQ=255)
                                      aTGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGG‑TAATCCACCGCGACCAGa                                                                                                                  <  2:373942/150‑1 (MQ=255)
                                          cTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGG‑TAAACCACCGCGACCAGATACa                                                                                                              >  1:1192062/1‑150 (MQ=255)
                                            tGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGTAAGAATAAACCGGGCAGGG‑GAATTCACCCCGACCCGAGACaaac                                                                                                           >  1:225619/1‑148 (MQ=255)
                                              tACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCCTACCGGTCAGGG‑GAATCCACCCccaccaaaaaaaaccc                                                                                                          >  2:682860/1‑140 (MQ=255)
                                              tACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGGCAGGG‑GAATCCACCCCGACCAGATACAAccc                                                                                                          >  1:1058565/1‑150 (MQ=255)
                                               aCTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGG‑TAATCCACCGCGACCAGATACAGCCCt                                                                                                         <  2:115479/150‑1 (MQ=255)
                                                  cGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGG‑TAATCCACCGCGACCAGATACAGCCCTTcc                                                                                                      <  1:933100/150‑1 (MQ=255)
                                                      aTCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCCGGCTGGG‑GAATACACCGCGACCAGAGAAAGCCCCTCCGCttt                                                                                                 >  2:1195841/1‑151 (MQ=255)
                                                          ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAGGATTATAACGGGCCGGG‑GAATCCCCCCCGACCAGATAAAGCCCCTCCCCcttt                                                                                                >  1:1010270/1‑148 (MQ=255)
                                                          ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGG‑TAATCCACCGCGACCAGATACAGCCCTTCCGCtttt                                                                                                <  2:1154884/148‑1 (MQ=255)
                                                          ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGG‑TAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGc                                                                                              <  2:262184/150‑1 (MQ=255)
                                                           tCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAGGATTATATCCGTCAGGG‑TGATTCCCCCCGACCCGAAATAGCCCCTCCGCCTTTTCC‑g                                                                                           >  1:500564/1‑151 (MQ=255)
                                                             aGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGG‑TAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCC‑Gt                                                                                          <  2:554217/150‑1 (MQ=255)
                                                                        aCGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCAGCCAGAAATAGCGGGCCCATTGGCGGGTTTTCAACATTATATCGGGCAGGG‑GAATCCCCCCCGGCCCAACAAAACCCTTCCTCATTTTCC‑Gtttttgt                                                                                    >  1:1088570/1‑140 (MQ=255)
                                                                          gTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACATACCCGTCAGGG‑GAATCCCCCCCGACCAGATACAGCCCCTCCCCTTTTTCC‑CTTGCTTGCCCCAGc                                                                             >  1:584931/1‑150 (MQ=255)
                                                                                tACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCCTACCCGGCAGGG‑GAAACCACCCCGACCCGATACAGCCCTTCCCCTTTTTCC‑GTTGCTGCCGCCCGCGctct                                                                        >  1:383363/1‑149 (MQ=255)
                                                                                 aCAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGG‑TAATCCGCCGCGACCATATACAGCCCTTCCGCTTTTGCC‑GTTGCTGCCGCCAGCGTTCTg                                                                       <  1:295252/149‑1 (MQ=255)
                                                                                 aCAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGG‑TAAACCACCGCGACCAGATACAGCCCTTCCGCTTTTGCC‑GTTGCTGCCGCCAGCGTTCTg                                                                       <  1:796579/149‑1 (MQ=255)
                                                                                   aGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGG‑TAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCC‑GTTGCTGCCGCCAGCGTTCTG‑TCCt                                                                  <  1:225565/151‑1 (MQ=255)
                                                                                    gCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGG‑GAATCCACCCCGACCAGAAAAAGCCCCTTCGCCTTTTCC‑GTTTCTGCCCCCCGCGTTCTG‑TCCtt                                                                 >  2:706146/1‑151 (MQ=255)
                                                                                     ccTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGATTTTGGAGAACATAACGGCCGAGG‑GAAACCACCCCGACCAGCCACACCCCCTCCGCTTTTTCC‑GTTGCTGCCGCCAACGTTCTG‑TTCttt                                                                >  2:469385/1‑151 (MQ=255)
                                                                                      cTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTTGAAGATCATAACGGTCCGGG‑GAAACCACCCCGACCCGACACAGCCCCTCCGCTTTTTCC‑GGTGCTtccccccaccgtccgtcc                                                                    >  1:857949/1‑135 (MQ=255)
                                                                                        gAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGG‑TAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCC‑GTTGCTGCCGCCAGCGTTCTG‑TCCTTTgc                                                              <  2:149421/150‑1 (MQ=255)
                                                                                         aaaTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGAACAAAACCGTCAGGG‑GAAACCAACGCGGCCAGAAACAACCCTTTCGCTTTTGCC‑GTTGCTGCCGCCAACGGTctt‑tccctttcc                                                             >  2:181113/1‑150 (MQ=255)
                                                                                           aTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACATACCCGGCAGGG‑GAAACCACCCCGGCCCGGTACAGCCCTTCCGCTTTTTTC‑GCTGCTCCCCCCCGCGTTCTG‑TTCTTTtccgcc                                                          >  1:184654/1‑151 (MQ=255)
                                                                                             gTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCAGGG‑GAATCCACCCCGGCCAGATACAACCCCTCCGCTTTTTCC‑GTTGCTTCCCCCACCGGTCTG‑TCCTTTttcgccc                                                         >  2:493438/1‑149 (MQ=255)
                                                                                             gTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGAATATACCGGTCAGGG‑GAAATCACCCCGCACAGAAACAACCCTTCCGCCTTTGCC‑GTTGCTGCCCCCCGCGCTCTG‑TCCTTTtccgcc                                                          >  2:1163015/1‑149 (MQ=255)
                                                                                                    gCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAACCGGCAGGG‑GAAACCACCGCGACCAGAAACAGCCCTTCCGCCTTTTTC‑GTTGGTGCCGCCAGCGTTCTG‑TCCCTTGCCGCCAGCAAcac                                                  >  1:1051342/1‑148 (MQ=255)
                                                                                                    gCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGG‑TAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCC‑GTTGCTGCCGCCAGCGTTCTG‑TCCTTTGCCGCCAGCAActct                                                 <  1:520423/151‑1 (MQ=255)
                                                                                                       ttattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGTCAGGG‑GAATCCACCCCCACCAGATACCGCCCCTCCCCTTTTTCC‑CTTGCTTCCCCCAGCGCTCTG‑TCCCTTGCCCCCCGCaa                                                     >  1:835616/1‑144 (MQ=255)
                                                                                                       ttattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAATGGTGAGGG‑GGATCCCCCCCGCCCAGATAAAGCCCCTCCCCTTTTTCC‑GTTGCTGCCCCCAGCGGTCTG‑TTCTTTGCCGCCCGCCActc                                                  >  1:1020210/1‑147 (MQ=255)
                                                                                                        tattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGATTTTGAACAATAATCCGCGAAGGG‑TAATTCAACCCCACCCAAGATAACCCTTTCTGTTTTTTC‑TCTGCTGGCGCCACCGTTCTT‑TTCTTTTCCGCCaatgaccctgc                                               >  2:54767/1‑139 (MQ=255)
                                                                                                         attCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGG‑TAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCC‑GTTGCTGCCGCCAGCGTTCTG‑TCCTTTGCCGCCAGCAACTCTGCTAt                                            <  2:114265/151‑1 (MQ=255)
                                                                                                         attCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGAACATACCGGGCAGGG‑GAAACCACCGCGACCAGATAAAACCCTTCCGCCTTTGCC‑GTTGCTGCCCCCCGCGGTCTG‑TCCTTTTCCGCCAGCAACTCTg                                                >  2:972644/1‑147 (MQ=255)
                                                                                                           tctcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAAGTCATACCGGGGAGGG‑GGAACCCCCACGCACCGAGAAACCCCCTCCGCCTTTGGC‑GTTTCTGCCGCCCACGATCTGTTCCTTTgcccccgcaaaacttctttcc                                           >  1:472808/1‑136 (MQ=255)
                                                                                                            ctcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCCGTCAGGG‑GAATCCCCCGCGACCAGATACAGCCCCTCCGCCTTTTCC‑GTTGCTGCCGCCAGCGTTCTG‑TCCCTTTCCCCCAGCAACTCTGCTAttc                                          >  2:462396/1‑148 (MQ=255)
                                                                                                             tcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGGAAGCATACCGGTCAGGGTTAATCC‑CCGCGCCCCGATACACCCCCTCCTCTTTTTTC‑GTTGCTGCCGCCCGCGGTCTC‑TCCTTTGTCCCCCGCAACTCCGccctcc                                          >  1:942058/1‑144 (MQ=255)
                                                                                                              cGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGG‑TAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCC‑GTTGCTGCCGCCAGCGTTCTG‑TCCTTTGCCGCCAGCAACTCTGCTATcc                                          <  1:146686/148‑1 (MQ=255)
                                                                                                                  aGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGG‑TAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCC‑GTTGCTGCCGCCAGCGTTCTG‑TCCTTTGCCGCCAGCAACTCTGCTATCCAGctct                                    <  1:469372/150‑1 (MQ=255)
                                                                                                                     ccGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCAAAACGGGCAGGG‑GAAACCACCGCGACCCGATACAGCCCTTCCGCTTTTTTC‑GTTTGTGCCGCCCGCGGTTTG‑TCCTTTTCCGCCCGCAACCCTTCTATTCAGCTcccccg                                >  1:907452/1‑149 (MQ=255)
                                                                                                                     ccGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGG‑TAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCC‑GTTGCTGCCGCCAGCGTTCTG‑TCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCg                                 <  2:1063617/150‑1 (MQ=255)
                                                                                                                       gCCAGAAATAGCGGGCCCATTGGCGGATTTTGGAGGTCAAACCAGCCCGGG‑GAAAACACCCCCAACCGATAAAATCCCTCCCCTTTTTCCGGCTGTCGCCG‑CAGCCCTCTG‑TCCCTTTCCGCCAGCAACCCTGCTCTCCCCCTCTcccg                                >  1:153930/1‑147 (MQ=255)
                                                                                                                            aaaTAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGG‑TAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCC‑GTTGCTGCCGCCAGCGTTCTG‑TCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTtg                         <  2:831043/151‑1 (MQ=255)
                                                                                                                             gaTAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGG‑TAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCC‑GTTGCTGCCGCCAGCGTTCTG‑TCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTtgt                        <  1:535246/150‑1 (MQ=255)
                                                                                                                                aGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGG‑TAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCC‑GTTGCTGCCGCCAGCGTTCTG‑TCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTggg                     <  1:1080185/151‑1 (MQ=255)
                                                                                                                                      cccATTGGCGGTTTTGGAAGATCATACCGGTCAGGG‑TAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCC‑GTTGCTGCCGCCAGCGTTCTG‑TCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGa               <  2:667944/151‑1 (MQ=255)
                                                                                                                                             gCGGTTTTGGAAGATCATACCGGTCAGGG‑TAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCC‑GTTGCTGCCGCCAGCGTTCTG‑TCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCAt        >  2:750097/1‑151 (MQ=255)
                                                                                                                                               ggTTTTGGAAGATCATACCGGTCAGGG‑TAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCC‑GTTGCTGCCGCCAGCGTTCTG‑TCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCa      <  2:884870/151‑1 (MQ=255)
                                                                                                                                                 ttttGGAAGATCATACCGGTCAGGG‑TAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCC‑GTTGCTGCCGCCAGCGTTCTG‑TCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAgg    >  1:839859/1‑151 (MQ=255)
                                                                                                                                                    tGGAAGATCATACCGGTCAGGG‑TAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCC‑GTTGCTGCCGCCAGCGTTCTG‑TCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCt  <  2:607581/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
AAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGG‑TAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCC‑GTTGCTGCCGCCAGCGTTCTG‑TCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCT  >  NC_000913/2434712‑2435009

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