Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,434,8820TG52.5% 68.9 / 99.6 62Y252S (TAC→TCC) truAtRNA pseudouridine(38‑40) synthase
Reads supporting (aligned to +/- strand):  ref base T (6/23);  new base G (32/0);  total (39/23)
Fisher's exact test for biased strand distribution p-value = 1.27e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.00e-11

GTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACC  >  NC_000913/2434736‑2435028
                                                                                                                                                   |                                                                                                                                                    
gTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGGaa                                                                                                                                                    >  2:537859/1‑149 (MQ=255)
  ccATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATcc                                                                                                                                                 <  2:585652/150‑1 (MQ=255)
  ccATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCa                                                                                                                                                <  2:370124/151‑1 (MQ=255)
      atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACcgc                                                                                                                                            <  2:532648/151‑1 (MQ=255)
      atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACcgc                                                                                                                                            <  2:364876/151‑1 (MQ=255)
      atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACcg                                                                                                                                             <  2:695017/150‑1 (MQ=255)
         tGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAATCCAcc                                                                                                                                              >  1:545962/1‑146 (MQ=255)
          gCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGGCAGGGGAATCCACCCCGAc                                                                                                                                         >  2:399753/1‑150 (MQ=255)
          gCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATAACGGGCAGGGGAATCCACCGCGAcc                                                                                                                                        >  1:108867/1‑151 (MQ=255)
           ccTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAAATCACCGCGAcc                                                                                                                                        >  2:849263/1‑150 (MQ=255)
             tATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAATCCACCCCGAcc                                                                                                                                        >  1:914069/1‑148 (MQ=255)
             tATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATAACCGGCAGGGGAATCCACCGCGAcc                                                                                                                                        >  2:823007/1‑148 (MQ=255)
                gCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACATACCGGTCAGGGGAATCCACCCCGGCCAGATa                                                                                                                                   >  2:751595/1‑150 (MQ=255)
                     gTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACATACCCGGCAGGGGAATACACCCCGCCCAGATAAAAccc                                                                                                                            >  1:959885/1‑151 (MQ=255)
                       aCTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCtt                                                                                                                          <  2:449858/151‑1 (MQ=255)
                       aCTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCt                                                                                                                           <  2:329031/150‑1 (MQ=255)
                       aCTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTGGAAGACCATACCGGCCAGGGTAATACCCCCCGACCAGGTACAAccc                                                                                                                            >  1:57086/1‑149 (MQ=255)
                          cGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCg                                                                                                                      <  1:267439/151‑1 (MQ=255)
                          cGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTcc                                                                                                                       >  1:202425/1‑150 (MQ=255)
                            tCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATAATAACGGGCAGGGGAAACCCCCGCGACCAGaaaaaaccc                                                                                                                             >  1:350001/1‑136 (MQ=255)
                            tCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGATTTGGTAGATAATATCCGGCAGGGGAATCCACCGCGACCAGATACAGCACCTCCGcc                                                                                                                    >  2:154842/1‑150 (MQ=255)
                              aTCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCAGGGGAAACCACCGCGACCAGATACAGCCCCTCCGCttt                                                                                                                  >  2:919022/1‑151 (MQ=255)
                               tCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATAATAACGGGCAGGGGAATCCACCGCGACCCGATACCGCCCTTCCCCcttt                                                                                                                 >  1:62052/1‑151 (MQ=255)
                               tCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTCCCGCt                                                                                                                    <  2:869955/148‑1 (MQ=255)
                                  ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGCTAGGGGAATCCACCCCGGACAGATACAAccccttcccctt                                                                                                                   >  2:84805/1‑147 (MQ=255)
                                  ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCtttt                                                                                                                 <  1:200784/148‑1 (MQ=255)
                                        cTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGAACATACCCGGCCGGGGAATCCCCCCCGCCCCGAAACAAGCCCTCCCCTTTTTCCCGtg                                                                                                           >  1:989916/1‑148 (MQ=255)
                                            aaCTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATAACCGGCAGGGGAAACCACCCCGACCCGATACAGCCCTTCCCCCTTTTGCCGt                                                                                                            >  1:259936/1‑144 (MQ=255)
                                                 cGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCa                                                                                                 <  1:768563/149‑1 (MQ=255)
                                                  gTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGTCAGGGGAATCCACCGCGACCAGAAACAGCCCCTCCGCTTTTTCCGTTGCTGCCGCCcgcg                                                                                              >  1:236601/1‑151 (MQ=255)
                                                               tGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATTATACCGGTCAGGGGAATACACCGCGACCAGATACAGCCCCTCCGCCTTTTCCCGTGCTGCCCCCCGCGGTTTGTTCttt                                                                                   >  2:685401/1‑149 (MQ=255)
                                                                gAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCAGGGGAATACACCCCGACCAGATACAACCCCTCCCCCTTTTCCCTTGCTGCCGCCCGCGGTCTGTCCttt                                                                                   >  2:225502/1‑148 (MQ=255)
                                                                   aTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACAAACCCGGCAGGGGAAAACACCGCGACCAGGTACAGCCCCTCCGCTTTTTCCGTTGGTGGCGCCCGCGttttttgccttgtcaac                                                                             >  1:857119/1‑136 (MQ=255)
                                                                     gTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCa                                                                            <  1:8121/150‑1 (MQ=255)
                                                                      ttCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGGCAGGGGAATCCACCCCGACCAGATAAAGCCCCTCCGCCTTTGCCGTTTCTGCCCCCAGCGGTCTGTCCTTTgccgcc                                                                             >  2:858792/1‑148 (MQ=255)
                                                                         cGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCTGGCAGGGGAATCCCCCGCGACCAGATACAGCCCTTTCGCTTTTGCCGTTGTGGCCGCCCACGATCTGTCCTTTTCCGCCCGCaa                                                                        >  1:493242/1‑150 (MQ=255)
                                                                         cGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCCGGCCGGGGAATCCACCCCGACCCGATACAGCCCTTCCGCCTTTTCCCTTGGTGCCCCCCGCGGTCTTTCCCTTtccgccc                                                                            >  2:689082/1‑145 (MQ=255)
                                                                           aGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAActc                                                                     <  1:369231/151‑1 (MQ=255)
                                                                             cATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAATCCACCGCGACCCGATACAGCCCCTCCGCTTTTTCCGGTGCTGCCGCCCGCGTTCTGTCCTTTTCCGCCCGCCAccc                                                                     >  2:930707/1‑147 (MQ=255)
                                                                             cATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTg                                                                   <  1:399774/151‑1 (MQ=255)
                                                                                 attCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATAACGGGCCGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCCCCCGCGTTCTGTCCTTTGCCGCCAGCAAATCTTCTAt                                                               >  1:952215/1‑151 (MQ=255)
                                                                                     tcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATAATACCGGGCAGGGGAATCCACCGCGGCCAGATACAGCCCTTCCGCTTTTGCCGTTGCTTCCCCCCGCGTTCTGTCCCTTTCCGCCCGCAACCCTGCTATcc                                                             >  1:789752/1‑149 (MQ=255)
                                                                                       gTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGct                                                         <  1:683835/151‑1 (MQ=255)
                                                                                           gTCCGCCTGAAAAAGCGGGCCCATTGGCGGTTTAGGAAGTTCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCATctct                                                       <  2:746476/149‑1 (MQ=255)
                                                                                             ccGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCACCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCgg                                                   <  1:931205/151‑1 (MQ=255)
                                                                                                ccAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTg                                                <  2:79886/151‑1 (MQ=255)
                                                                                                 cAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATAAAAACGGTCAGGGGAATACACCCCGACCCGAAACAACCCTTCCGCTTTTTCCGTTGCTGCCGCCAGCGGTCTGTCCTTTTCCGCCAGCAACTCCGTTTTTCAGCTTCCCgggttg                                               >  1:1018406/1‑147 (MQ=255)
                                                                                                 cAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACATAGCGGGCCAGGGAAATCACCCCCACCAGATACAGGCCCTCCCCCTTTGCCGTTTCTGCCGCCCGCGTTCGGGCCTTTGCCGCCAGCTAATCTGCGCTCCAGCTCTccggccgg                                               >  2:279845/1‑151 (MQ=255)
                                                                                                   gAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGtt                                             <  2:1046202/151‑1 (MQ=255)
                                                                                                    aaaTAGCGGGCCCATTGGCGGTTTTTGGAAAACATACCGGGCAGGGGAATTCACCCCGACCAGATACAACCCCTCCGCTTTTTCCGTTGCTTCCGCCCGCGTTCTGTCCTTTTCCGCCCGCAACTCTTCTATCCAGCTCTCCCGCTGGtt                                             >  2:548353/1‑150 (MQ=255)
                                                                                                     aaTAGCGGGCCCATTGGCGGTTTTTGGAGAAAATAACGGGGAGGGGAAATCACCGCGACCAGATAAAGCCCCTCCGCTTTTTCCGGTGCTGCCGCCCGCGGTCTGTCCTTTGCCGCCAGCAACTCTGCTATTCAGCTCTCCGGCTgggt                                             >  1:742160/1‑147 (MQ=255)
                                                                                                     aaTAGCGGGCCCATTGGCGGTTTTGGAAGGACAAACCGGTCAGGGGAATCCACCCCGACCCGATACAGCCCCTTCGCCTTTTCCGGTGCTGCCCCCCGCGTTCTTCCCTTTGCCCCCAGCCAATCTTCTTTCCCGCTCTCCCGCttggt                                             >  2:905847/1‑145 (MQ=255)
                                                                                                        aGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTgg                                         <  2:824693/150‑1 (MQ=255)
                                                                                                               ccATTGGCGGTTTTGGAAGATCATAACGGGCCGGGGAAACCACCGCGACCAGATAAAGCACCTCCGCCTTTTCCGTTGGCGCCCCCCGCGGTCTGTCCTTTTTCCCCCGCAACTCTGCTTTCCCGCCCTCCGGCTGGTTTGGGGCACCGa                                  >  2:550376/1‑150 (MQ=255)
                                                                                                                      cGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATc                          >  1:107603/1‑151 (MQ=255)
                                                                                                                                  aTCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACa              <  1:537888/151‑1 (MQ=255)
                                                                                                                                   tCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACa              >  2:363746/1‑150 (MQ=255)
                                                                                                                                     aTACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACCat            >  1:464115/1‑150 (MQ=255)
                                                                                                                                      tACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAatat          >  2:762030/1‑151 (MQ=255)
                                                                                                                                         ccGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTc        <  1:804402/151‑1 (MQ=255)
                                                                                                                                           ggTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCt      <  2:925965/151‑1 (MQ=255)
                                                                                                                                               aGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAcc  >  2:353522/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACC  >  NC_000913/2434736‑2435028

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: