Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,434,8820TG43.9% 107.2 / 18.5 66Y252S (TAC→TCC) truAtRNA pseudouridine(38‑40) synthase
Reads supporting (aligned to +/- strand):  ref base T (10/27);  new base G (29/0);  total (39/27)
Fisher's exact test for biased strand distribution p-value = 1.49e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTG‑CCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCAT  >  NC_000913/2434733‑2435030
                                                                                                                                                      |                                                                                                                                                     
cAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGt                                                                                                                                                       <  2:44388/150‑1 (MQ=255)
cAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGTa                                                                                                                                                      <  2:693821/151‑1 (MQ=255)
     ccATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGTAATcc                                                                                                                                                  <  1:172045/150‑1 (MQ=255)
          tttGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAAATCATAC‑CGGGCAGGGGAATCCAcc                                                                                                                                               >  2:628806/1‑148 (MQ=255)
           ttGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGATTTGGATGAACCTAA‑CAGTCCGGGGAATCCa                                                                                                                                                 >  2:924956/1‑145 (MQ=255)
             gCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAAATCATAC‑CGGTCAGGGGAATCCACCGCGAcc                                                                                                                                         >  2:415722/1‑151 (MQ=255)
              ccTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGTAATCCACCGCGACCa                                                                                                                                        <  1:839640/151‑1 (MQ=255)
                 aTGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGTAATCCACCGCGACCAGAt                                                                                                                                     <  1:107712/151‑1 (MQ=255)
                  tGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGGCAGGGGAATCCACCCCGACCAGAt                                                                                                                                     >  1:688253/1‑150 (MQ=255)
                       tGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAC‑CCGTCCGGGGAAACCACCGCGACCAGATACAGcc                                                                                                                               >  1:827751/1‑151 (MQ=255)
                        gTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGTAATCCACCGCGACCAGATACAGccc                                                                                                                              <  1:82328/151‑1 (MQ=255)
                          aCTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGAGCATAA‑CGGGCAGGGGAAACCACCCCGCCCAGGTAAAGccc                                                                                                                              >  1:688240/1‑149 (MQ=255)
                          aCTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGGAATCCACCGCGACCAGATACAGCCCtt                                                                                                                            >  2:373257/1‑151 (MQ=255)
                             cGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTcc                                                                                                                          <  1:1013141/150‑1 (MQ=255)
                             cGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTcc                                                                                                                          <  2:137976/150‑1 (MQ=255)
                                cATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATAC‑CGGGCAGGGGAATCCACCGCGACCCGATAAAGCCCTTCCGCtt                                                                                                                      >  1:1050458/1‑151 (MQ=255)
                                     ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTg                                                                                                                   <  2:978276/149‑1 (MQ=255)
                                      tCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTg                                                                                                                   <  1:269831/148‑1 (MQ=255)
                                        aGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATAC‑CCGGCAGGGGAATCCACCCCGGCCAGATACAGCCCCTCCGCtttt                                                                                                                    >  1:821692/1‑145 (MQ=255)
                                        aGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGtt                                                                                                              <  1:721873/151‑1 (MQ=255)
                                           cTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTTCCg                                                                                                                >  2:248600/1‑146 (MQ=255)
                                                   aCGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGGCAGGGGAATCCACCGCGACCCGATACAGCCCCTCCGCTTTTGCCGTTGCTg‑ccccc                                                                                                    >  2:155083/1‑149 (MQ=255)
                                                      ttGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGGAATCCACCCCGACCCGATACAGCCCTTCCCCTTTTTCCGTTGCTG‑CCGCCAGCg                                                                                                >  2:970905/1‑150 (MQ=255)
                                                            aCAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTG‑CCGCCAGCGTTCTGTc                                                                                         <  2:245406/151‑1 (MQ=255)
                                                              aGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAC‑CGGGCAGGGGAATCCACCCGGGCCAGATAAAGCCCCTTCGCTTTTGGCGTTGCTG‑CCCCCCGCGtttttttc                                                                                        >  1:206362/1‑144 (MQ=255)
                                                              aGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAAATCAATA‑CAGTGCGGGGAATACACCGCCGCCAGACACAGCCCTTCCGCCTTTGCCGTTTCTGCCCCCCAGC‑TTCTGTCCt                                                                                       >  2:99783/1‑151 (MQ=255)
                                                                ccTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATAC‑CGGGCAGGGGAAACCACCGCGACCAGATACAGCCCTTCCGCTTTTTCCGTTGCTT‑CCCCCCGCGTTTTGTCCttt                                                                                     >  1:897020/1‑151 (MQ=255)
                                                                ccTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTG‑CCGCCAGCGTTCTGTCCttt                                                                                     <  2:224578/151‑1 (MQ=255)
                                                                   gAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTG‑CCGCCAGCGTTCTGTCCTTTgcc                                                                                  <  2:407129/151‑1 (MQ=255)
                                                                   gAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGGAATCCACCCCGACCAGATACAGCCCTTCCGCTTTTTCCGTTGCTG‑CCGCCAGCGTTCTGTCCCTTgcc                                                                                  >  1:439880/1‑151 (MQ=255)
                                                                   gAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTTGACGATCATAC‑CCGTCAGGGGAATCCACCGCGACCAGATAAAGCCCCTCCCCCTTTTCCGCTGCTG‑CCCCCAACGTTCCGTCCTTTtcc                                                                                  >  2:117959/1‑151 (MQ=255)
                                                                       tGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATAC‑CCGGCAGGGGAAACCACCCCGGCCAGATACAGCCCCTCCGCTTTTGCCGTTGCTG‑CCGCCAGCGGTCTGTTCTTTtccgcc                                                                               >  1:438118/1‑150 (MQ=255)
                                                                        gTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATAC‑CCGGCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTTCCGTTGCTG‑CCGCCAGCGTTCTGTCCTTTtccgcc                                                                               >  2:824120/1‑149 (MQ=255)
                                                                          tCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGCTTTTGAAGATCATAA‑CCGTCAGGGGAATCCACCGCGGGCAGATACAGGCCTTTCCCTTTTGCCGTTGCTG‑CCGCCCGCGTTCTGTCCTTTtccgcc                                                                               >  1:326855/1‑147 (MQ=255)
                                                                             gAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATAC‑CGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTTCCGTTGCTG‑CCGCCAGCGTTCTGTCCTTTTCCGCCCGCaa                                                                          >  1:674669/1‑149 (MQ=255)
                                                                                cATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTG‑CCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTg                                                                     <  1:411779/151‑1 (MQ=255)
                                                                                cATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGGAATCCACCGCGACCAGAAACAGCCCCTCCGCTTTTGCCGTTGCTG‑CCGCCAGCGTTCTGTCCTTTGCCGCCAGCAActct                                                                      >  2:509431/1‑150 (MQ=255)
                                                                                 attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGGCAGGGGAATCCACCCCGACCAGATACAGCCCTTCCCCTTTTGCCGTTGCTG‑CCCCCCGCGTTCTGTCCTTTTCCGCCCGCa                                                                           >  2:846489/1‑144 (MQ=255)
                                                                                   tattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCCTACTCGGT‑AGGGTAATCCATCGCGACCCGATAAAGCCCTTCCGTTTTTGACCGTGCTG‑CCTCCAGCGGTCTTTCCTTTGCCGCCAGCACCtcttct                                                                   >  2:525077/1‑150 (MQ=255)
                                                                                    attCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTG‑CCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTAt                                                                 <  2:434691/151‑1 (MQ=255)
                                                                                          gTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGCAGATCATAC‑CCGGCAAGGGAATCCACCCCGACCAGGTACAGCCCCTCCGCTTTTGCCCCTTCTG‑CCGCCAGCGCTTTGTTCCTTTCCGCCCGCAAATCTGCTATcc                                                               >  2:828229/1‑147 (MQ=255)
                                                                                          gTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGGAATCCACCCCGACCAGATACAGCCCTTCCCCTTTTGCCGTTGCTG‑CCCCCCGCGTTCTGTCCTTTTCCGCCAGCAACTCTGCTATcc                                                               >  2:210630/1‑147 (MQ=255)
                                                                                            tAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTG‑CCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGctc                                                          <  1:517793/150‑1 (MQ=255)
                                                                                               tCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTG‑CCGCCCGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTTTCCAGCTCTcc                                                       >  2:322438/1‑150 (MQ=255)
                                                                                                ccGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTG‑CCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCg                                                      <  2:369227/150‑1 (MQ=255)
                                                                                                ccGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTG‑CCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCg                                                      <  2:201828/150‑1 (MQ=255)
                                                                                                 cGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTG‑CCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGc                                                    <  1:190013/151‑1 (MQ=255)
                                                                                                        aaTAGCGGGCCCATTGGCGGTTTTTGGAGATCATAC‑CGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTTCTG‑CCCCCAGCGTTCTGTCCTTTTCCGCCCGCAACTCTGCTATCCAGCTCTCCGGCTgggt                                               >  2:374393/1‑147 (MQ=255)
                                                                                                        aaTAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTG‑CCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTtgt                                             <  2:368505/151‑1 (MQ=255)
                                                                                                           aGCGGGCCCATTGGCGGTTTTTGAAGATCATAC‑CGGTCAGGGGAATCCACCCCGACCAGATACAGCCCTTCCGCTTTTTCCGTTTCTG‑CCGCCAGCGTTCTGTCCTTTTCCCCCCGCAACTCTGCTATCCAGCTCTCCCGCTGGTTGTggg                                          >  1:669698/1‑151 (MQ=255)
                                                                                                           aGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTG‑CCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTgg                                           <  2:821709/150‑1 (MQ=255)
                                                                                                            gCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTG‑CCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGc                                         <  1:665000/151‑1 (MQ=255)
                                                                                                                gCCCATTGGCGGTTTTGGAAGATCATAC‑CGGGCAGGGGAATCCACCCCGACCAGATACAGCCCTTCCGCTTTTGCCGGTGCTG‑CCGCCAGCGGTCTGGCCTTTGCCGCCCGCAACTCTGCTATCCAGCTCTCCCGCTGGTTGTGGGCAcc                                      >  1:41384/1‑150 (MQ=255)
                                                                                                                   cATTGGCGGTTTTTGGAGATCAAAC‑CGGTCAGGGTAATCCACCGCGACCCGATACAGCCCCTCCGCTTTTTCCGTTGCTG‑CCGCCCGCGTTCTGTCCTTTGCCGCCCGCAACTCCGCTATCCAGCTCTCCGGCTGGTTTTGGGCAcc                                      >  1:584633/1‑147 (MQ=255)
                                                                                                                       ggCGGTTTTGGAAGATCATAC‑CGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTG‑CCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTcc                               <  1:146541/150‑1 (MQ=255)
                                                                                                                        gCGGTTTTGGAAGATCATAC‑CGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTG‑CCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCAt                             >  2:399714/1‑151 (MQ=255)
                                                                                                                           gTTTTGGAAGATCATAC‑CGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTG‑CCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATc                            >  1:290212/1‑149 (MQ=255)
                                                                                                                               tGGAAGATCATAC‑CGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTG‑CCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTg                      >  1:106194/1‑151 (MQ=255)
                                                                                                                                ggAAGATCATAC‑CGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTG‑CCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTg                      <  2:307487/150‑1 (MQ=255)
                                                                                                                                   aGATCATAC‑CGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTG‑CCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGa                  <  1:628794/151‑1 (MQ=255)
                                                                                                                                      tCATAC‑CGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTG‑CCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACa                >  2:867879/1‑150 (MQ=255)
                                                                                                                                       cATAC‑CGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTG‑CCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACaa               <  2:603961/150‑1 (MQ=255)
                                                                                                                                       cATAC‑CGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTG‑CCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAat              >  2:310776/1‑151 (MQ=255)
                                                                                                                                           c‑cGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTG‑CCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTc          >  2:792733/1‑151 (MQ=255)
                                                                                                                                                tCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTG‑CCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTg       >  1:140096/1‑150 (MQ=255)
                                                                                                                                                  aGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTG‑CCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAcc    >  1:643430/1‑151 (MQ=255)
                                                                                                                                                    ggTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTG‑CCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCat  <  1:209926/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAC‑CGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTG‑CCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCAT  >  NC_000913/2434733‑2435030

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