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

CAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCAT  >  NC_000913/2434733‑2435030
                                                                                                                                                      |                                                                                                                                                     
cAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGt                                                                                                                                                       <  2:44388/150‑1 (MQ=255)
cAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTa                                                                                                                                                      <  2:693821/151‑1 (MQ=255)
     ccATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATcc                                                                                                                                                  <  1:172045/150‑1 (MQ=255)
          tttGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAAATCATACCGGGCAGGGGAATCCAcc                                                                                                                                               >  2:628806/1‑148 (MQ=255)
           ttGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGATTTGGATGAACCTAACAGTCCGGGGAATCCa                                                                                                                                                 >  2:924956/1‑145 (MQ=255)
             gCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAAATCATACCGGTCAGGGGAATCCACCGCGAcc                                                                                                                                         >  2:415722/1‑151 (MQ=255)
              ccTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCa                                                                                                                                        <  1:839640/151‑1 (MQ=255)
                 aTGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGAt                                                                                                                                     <  1:107712/151‑1 (MQ=255)
                  tGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGGCAGGGGAATCCACCCCGACCAGAt                                                                                                                                     >  1:688253/1‑150 (MQ=255)
                       tGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGTCCGGGGAAACCACCGCGACCAGATACAGcc                                                                                                                               >  1:827751/1‑151 (MQ=255)
                        gTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGccc                                                                                                                              <  1:82328/151‑1 (MQ=255)
                          aCTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGAGCATAACGGGCAGGGGAAACCACCCCGCCCAGGTAAAGccc                                                                                                                              >  1:688240/1‑149 (MQ=255)
                          aCTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCtt                                                                                                                            >  2:373257/1‑151 (MQ=255)
                             cGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTcc                                                                                                                          <  1:1013141/150‑1 (MQ=255)
                             cGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTcc                                                                                                                          <  2:137976/150‑1 (MQ=255)
                                cATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATCCACCGCGACCCGATAAAGCCCTTCCGCtt                                                                                                                      >  1:1050458/1‑151 (MQ=255)
                                     ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTg                                                                                                                   <  2:978276/149‑1 (MQ=255)
                                      tCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTg                                                                                                                   <  1:269831/148‑1 (MQ=255)
                                        aGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCCGGCAGGGGAATCCACCCCGGCCAGATACAGCCCCTCCGCtttt                                                                                                                    >  1:821692/1‑145 (MQ=255)
                                        aGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGtt                                                                                                              <  1:721873/151‑1 (MQ=255)
                                           cTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTTCCg                                                                                                                >  2:248600/1‑146 (MQ=255)
                                                   aCGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGGCAGGGGAATCCACCGCGACCCGATACAGCCCCTCCGCTTTTGCCGTTGCTgccccc                                                                                                    >  2:155083/1‑149 (MQ=255)
                                                      ttGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGGAATCCACCCCGACCCGATACAGCCCTTCCCCTTTTTCCGTTGCTGCCGCCAGCg                                                                                                >  2:970905/1‑150 (MQ=255)
                                                            aCAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTc                                                                                         <  2:245406/151‑1 (MQ=255)
                                                              aGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAATCCACCCGGGCCAGATAAAGCCCCTTCGCTTTTGGCGTTGCTGCCCCCCGCGtttttttc                                                                                        >  1:206362/1‑144 (MQ=255)
                                                              aGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAAATCAATACAGTGCGGGGAATACACCGCCGCCAGACACAGCCCTTCCGCCTTTGCCGTTTCTGCCCCCCAGCTTCTGTCCt                                                                                       >  2:99783/1‑151 (MQ=255)
                                                                ccTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAAACCACCGCGACCAGATACAGCCCTTCCGCTTTTTCCGTTGCTTCCCCCCGCGTTTTGTCCttt                                                                                     >  1:897020/1‑151 (MQ=255)
                                                                ccTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCttt                                                                                     <  2:224578/151‑1 (MQ=255)
                                                                   gAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTgcc                                                                                  <  2:407129/151‑1 (MQ=255)
                                                                   gAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGGAATCCACCCCGACCAGATACAGCCCTTCCGCTTTTTCCGTTGCTGCCGCCAGCGTTCTGTCCCTTgcc                                                                                  >  1:439880/1‑151 (MQ=255)
                                                                   gAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTTGACGATCATACCCGTCAGGGGAATCCACCGCGACCAGATAAAGCCCCTCCCCCTTTTCCGCTGCTGCCCCCAACGTTCCGTCCTTTtcc                                                                                  >  2:117959/1‑151 (MQ=255)
                                                                       tGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCCGGCAGGGGAAACCACCCCGGCCAGATACAGCCCCTCCGCTTTTGCCGTTGCTGCCGCCAGCGGTCTGTTCTTTtccgcc                                                                               >  1:438118/1‑150 (MQ=255)
                                                                        gTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCCGGCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTTCCGTTGCTGCCGCCAGCGTTCTGTCCTTTtccgcc                                                                               >  2:824120/1‑149 (MQ=255)
                                                                          tCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGCTTTTGAAGATCATAACCGTCAGGGGAATCCACCGCGGGCAGATACAGGCCTTTCCCTTTTGCCGTTGCTGCCGCCCGCGTTCTGTCCTTTtccgcc                                                                               >  1:326855/1‑147 (MQ=255)
                                                                             gAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTTCCGTTGCTGCCGCCAGCGTTCTGTCCTTTTCCGCCCGCaa                                                                          >  1:674669/1‑149 (MQ=255)
                                                                                cATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTg                                                                     <  1:411779/151‑1 (MQ=255)
                                                                                cATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGGAATCCACCGCGACCAGAAACAGCCCCTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAActct                                                                      >  2:509431/1‑150 (MQ=255)
                                                                                 attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGGCAGGGGAATCCACCCCGACCAGATACAGCCCTTCCCCTTTTGCCGTTGCTGCCCCCCGCGTTCTGTCCTTTTCCGCCCGCa                                                                           >  2:846489/1‑144 (MQ=255)
                                                                                   tattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCCTACTCGGTAGGGTAATCCATCGCGACCCGATAAAGCCCTTCCGTTTTTGACCGTGCTGCCTCCAGCGGTCTTTCCTTTGCCGCCAGCACCtcttct                                                                   >  2:525077/1‑150 (MQ=255)
                                                                                    attCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTAt                                                                 <  2:434691/151‑1 (MQ=255)
                                                                                          gTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGCAGATCATACCCGGCAAGGGAATCCACCCCGACCAGGTACAGCCCCTCCGCTTTTGCCCCTTCTGCCGCCAGCGCTTTGTTCCTTTCCGCCCGCAAATCTGCTATcc                                                               >  2:828229/1‑147 (MQ=255)
                                                                                          gTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGGAATCCACCCCGACCAGATACAGCCCTTCCCCTTTTGCCGTTGCTGCCCCCCGCGTTCTGTCCTTTTCCGCCAGCAACTCTGCTATcc                                                               >  2:210630/1‑147 (MQ=255)
                                                                                            tAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGctc                                                          <  1:517793/150‑1 (MQ=255)
                                                                                               tCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCCGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTTTCCAGCTCTcc                                                       >  2:322438/1‑150 (MQ=255)
                                                                                                ccGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCg                                                      <  2:369227/150‑1 (MQ=255)
                                                                                                ccGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCg                                                      <  2:201828/150‑1 (MQ=255)
                                                                                                 cGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGc                                                    <  1:190013/151‑1 (MQ=255)
                                                                                                        aaTAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTTCTGCCCCCAGCGTTCTGTCCTTTTCCGCCCGCAACTCTGCTATCCAGCTCTCCGGCTgggt                                               >  2:374393/1‑147 (MQ=255)
                                                                                                        aaTAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTtgt                                             <  2:368505/151‑1 (MQ=255)
                                                                                                           aGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGGAATCCACCCCGACCAGATACAGCCCTTCCGCTTTTTCCGTTTCTGCCGCCAGCGTTCTGTCCTTTTCCCCCCGCAACTCTGCTATCCAGCTCTCCCGCTGGTTGTggg                                          >  1:669698/1‑151 (MQ=255)
                                                                                                           aGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTgg                                           <  2:821709/150‑1 (MQ=255)
                                                                                                            gCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGc                                         <  1:665000/151‑1 (MQ=255)
                                                                                                                gCCCATTGGCGGTTTTGGAAGATCATACCGGGCAGGGGAATCCACCCCGACCAGATACAGCCCTTCCGCTTTTGCCGGTGCTGCCGCCAGCGGTCTGGCCTTTGCCGCCCGCAACTCTGCTATCCAGCTCTCCCGCTGGTTGTGGGCAcc                                      >  1:41384/1‑150 (MQ=255)
                                                                                                                   cATTGGCGGTTTTTGGAGATCAAACCGGTCAGGGTAATCCACCGCGACCCGATACAGCCCCTCCGCTTTTTCCGTTGCTGCCGCCCGCGTTCTGTCCTTTGCCGCCCGCAACTCCGCTATCCAGCTCTCCGGCTGGTTTTGGGCAcc                                      >  1:584633/1‑147 (MQ=255)
                                                                                                                       ggCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTcc                               <  1:146541/150‑1 (MQ=255)
                                                                                                                        gCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCAt                             >  2:399714/1‑151 (MQ=255)
                                                                                                                           gTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATc                            >  1:290212/1‑149 (MQ=255)
                                                                                                                               tGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTg                      >  1:106194/1‑151 (MQ=255)
                                                                                                                                ggAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTg                      <  2:307487/150‑1 (MQ=255)
                                                                                                                                   aGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGa                  <  1:628794/151‑1 (MQ=255)
                                                                                                                                      tCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACa                >  2:867879/1‑150 (MQ=255)
                                                                                                                                       cATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACaa               <  2:603961/150‑1 (MQ=255)
                                                                                                                                       cATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAat              >  2:310776/1‑151 (MQ=255)
                                                                                                                                           ccGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTc          >  2:792733/1‑151 (MQ=255)
                                                                                                                                                tCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTg       >  1:140096/1‑150 (MQ=255)
                                                                                                                                                  aGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAcc    >  1:643430/1‑151 (MQ=255)
                                                                                                                                                    ggTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCat  <  1:209926/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCAT  >  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: