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

GTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑AGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACC  >  NC_000913/2434736‑2435028
                                                                                                                                                   |                                                                                                                                                    
gTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATA‑CCGGTCAGGGGaa                                                                                                                                                    >  2:537859/1‑149 (MQ=255)
  ccATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATcc                                                                                                                                                 <  2:585652/150‑1 (MQ=255)
  ccATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCa                                                                                                                                                <  2:370124/151‑1 (MQ=255)
      atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACcgc                                                                                                                                            <  2:532648/151‑1 (MQ=255)
      atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACcgc                                                                                                                                            <  2:364876/151‑1 (MQ=255)
      atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACcg                                                                                                                                             <  2:695017/150‑1 (MQ=255)
         tGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATA‑CCGGGCAGGGGAATCCAcc                                                                                                                                              >  1:545962/1‑146 (MQ=255)
          gCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATA‑CCCGGCAGGGGAATCCACCCCGAc                                                                                                                                         >  2:399753/1‑150 (MQ=255)
          gCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATA‑ACGGGCAGGGGAATCCACCGCGAcc                                                                                                                                        >  1:108867/1‑151 (MQ=255)
           ccTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATA‑CCGGGCAGGGGAAATCACCGCGAcc                                                                                                                                        >  2:849263/1‑150 (MQ=255)
             tATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATA‑CCGGGCAGGGGAATCCACCCCGAcc                                                                                                                                        >  1:914069/1‑148 (MQ=255)
             tATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATA‑ACCGGCAGGGGAATCCACCGCGAcc                                                                                                                                        >  2:823007/1‑148 (MQ=255)
                gCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACATA‑CCGGTCAGGGGAATCCACCCCGGCCAGATa                                                                                                                                   >  2:751595/1‑150 (MQ=255)
                     gTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACATA‑CCCGGCAGGGGAATACACCCCGCCCAGATAA‑AAccc                                                                                                                            >  1:959885/1‑151 (MQ=255)
                       aCTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑AGCCCtt                                                                                                                          <  2:449858/151‑1 (MQ=255)
                       aCTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑AGCCCt                                                                                                                           <  2:329031/150‑1 (MQ=255)
                       aCTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTGGAAGACCATA‑CCGGCCAGGGTAATACCCCCCGACCAGGTAC‑AAccc                                                                                                                            >  1:57086/1‑149 (MQ=255)
                          cGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑AGCCCTT‑CCg                                                                                                                      <  1:267439/151‑1 (MQ=255)
                          cGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGGAATCCACCGCGACCAGATAC‑AGCCCTT‑cc                                                                                                                       >  1:202425/1‑150 (MQ=255)
                            tCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATAATA‑ACGGGCAGGGGAAACCCCCGCGACCAGaaaaaaccc                                                                                                                             >  1:350001/1‑136 (MQ=255)
                            tCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGATTTGGTAGATAATA‑TCCGGCAGGGGAATCCACCGCGACCAGATAC‑AGCACCT‑CCGcc                                                                                                                    >  2:154842/1‑150 (MQ=255)
                              aTCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATA‑CCGGTCAGGGGAAACCACCGCGACCAGATAC‑AGCCCCT‑CCGCttt                                                                                                                  >  2:919022/1‑151 (MQ=255)
                               tCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATAATA‑ACGGGCAGGGGAATCCACCGCGACCCGATAC‑CGCCCTT‑CCCCcttt                                                                                                                 >  1:62052/1‑151 (MQ=255)
                               tCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑AGCCCTC‑CCGCt                                                                                                                    <  2:869955/148‑1 (MQ=255)
                                  ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATA‑CCGGCTAGGGGAATCCACCCCGGACAGATACAAcccctt‑cccctt                                                                                                                   >  2:84805/1‑147 (MQ=255)
                                  ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑AGCCCTT‑CCGCtttt                                                                                                                 <  1:200784/148‑1 (MQ=255)
                                        cTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGAACATA‑CCCGGCCGGGGAATCCCCCCCGCCCCGAAACAAGCCCTC‑CCCTTTTTCCCGtg                                                                                                           >  1:989916/1‑148 (MQ=255)
                                            aaCTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATA‑ACCGGCAGGGGAAACCACCCCGACCCGATAC‑AGCCCTTCCCCCTTTTGCCGt                                                                                                            >  1:259936/1‑144 (MQ=255)
                                                 cGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑AGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCa                                                                                                 <  1:768563/149‑1 (MQ=255)
                                                  gTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATA‑CCGGTCAGGGGAATCCACCGCGACCAGAAAC‑AGCCCCT‑CCGCTTTTTCCGTTGCTGCCGCCcgcg                                                                                              >  1:236601/1‑151 (MQ=255)
                                                               tGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATTATA‑CCGGTCAGGGGAATACACCGCGACCAGATAC‑AGCCCCT‑CCGCCTTTTCCCGTGCTGCCCCCCGCGGTTTGTTCttt                                                                                   >  2:685401/1‑149 (MQ=255)
                                                                gAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATA‑CCGGTCAGGGGAATACACCCCGACCAGATAC‑AACCCCT‑CCCCCTTTTCCCTTGCTGCCGCCCGCGGTCTGTCCttt                                                                                   >  2:225502/1‑148 (MQ=255)
                                                                   aTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACAAA‑CCCGGCAGGGGAAAACACCGCGACCAGGTAC‑AGCCCCT‑CCGCTTTTTCCGTTGGTGGCGCCCGCGttttttgccttgtcaac                                                                             >  1:857119/1‑136 (MQ=255)
                                                                     gTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑AGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCa                                                                            <  1:8121/150‑1 (MQ=255)
                                                                      ttCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATA‑CCGGGCAGGGGAATCCACCCCGACCAGATAA‑AGCCCCT‑CCGCCTTTGCCGTTTCTGCCCCCAGCGGTCTGTCCTTTgccgcc                                                                             >  2:858792/1‑148 (MQ=255)
                                                                         cGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATA‑CCTGGCAGGGGAATCCCCCGCGACCAGATAC‑AGCCCTT‑TCGCTTTTGCCGTTGTGGCCGCCCACGATCTGTCCTTTTCCGCCCGCaa                                                                        >  1:493242/1‑150 (MQ=255)
                                                                         cGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATA‑CCCGGCCGGGGAATCCACCCCGACCCGATAC‑AGCCCTT‑CCGCCTTTTCCCTTGGTGCCCCCCGCGGTCTTTCCCTTtccgccc                                                                            >  2:689082/1‑145 (MQ=255)
                                                                           aGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑AGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAActc                                                                     <  1:369231/151‑1 (MQ=255)
                                                                             cATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATA‑CCGGGCAGGGGAATCCACCGCGACCCGATAC‑AGCCCCT‑CCGCTTTTTCCGGTGCTGCCGCCCGCGTTCTGTCCTTTTCCGCCCGCCAccc                                                                     >  2:930707/1‑147 (MQ=255)
                                                                             cATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑AGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTg                                                                   <  1:399774/151‑1 (MQ=255)
                                                                                 attCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATA‑ACGGGCCGGGGAATCCACCGCGACCAGATAC‑AGCCCTT‑CCGCTTTTGCCGTTGCTGCCCCCCGCGTTCTGTCCTTTGCCGCCAGCAAATCTTCTAt                                                               >  1:952215/1‑151 (MQ=255)
                                                                                     tcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATAATA‑CCGGGCAGGGGAATCCACCGCGGCCAGATAC‑AGCCCTT‑CCGCTTTTGCCGTTGCTTCCCCCCGCGTTCTGTCCCTTTCCGCCCGCAACCCTGCTATcc                                                             >  1:789752/1‑149 (MQ=255)
                                                                                       gTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑AGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGct                                                         <  1:683835/151‑1 (MQ=255)
                                                                                           gTCCGCCTGAAAAAGCGGGCCCATTGGCGGTTTAGGAAGTTCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑AGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCATctct                                                       <  2:746476/149‑1 (MQ=255)
                                                                                             ccGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑AGCCCTT‑CCGCTTTTGCCGTTGCTGCCACCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCgg                                                   <  1:931205/151‑1 (MQ=255)
                                                                                                ccAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑AGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTg                                                <  2:79886/151‑1 (MQ=255)
                                                                                                 cAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATAAAA‑ACGGTCAGGGGAATACACCCCGACCCGAAAC‑AACCCTT‑CCGCTTTTTCCGTTGCTGCCGCCAGCGGTCTGTCCTTTTCCGCCAGCAACTCCGTTTTTCAGCTTCCCgggttg                                               >  1:1018406/1‑147 (MQ=255)
                                                                                                 cAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACATAGCGGGCCAGGGAAAT‑CACCCCCACCAGATAC‑AGGCCCT‑CCCCCTTTGCCGTTTCTGCCGCCCGCGTTCGGGCCTTTGCCGCCAGCTAATCTGCGCTCCAGCTCTccggccgg                                               >  2:279845/1‑151 (MQ=255)
                                                                                                   gAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑AGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGtt                                             <  2:1046202/151‑1 (MQ=255)
                                                                                                    aaaTAGCGGGCCCATTGGCGGTTTTTGGAAAACATA‑CCGGGCAGGGGAATTCACCCCGACCAGATAC‑AACCCCT‑CCGCTTTTTCCGTTGCTTCCGCCCGCGTTCTGTCCTTTTCCGCCCGCAACTCTTCTATCCAGCTCTCCCGCTGGtt                                             >  2:548353/1‑150 (MQ=255)
                                                                                                     aaTAGCGGGCCCATTGGCGGTTTTTGGAGAAAATA‑ACGGGGAGGGGAAATCACCGCGACCAGATAA‑AGCCCCT‑CCGCTTTTTCCGGTGCTGCCGCCCGCGGTCTGTCCTTTGCCGCCAGCAACTCTGCTATTCAGCTCTCCGGCTgggt                                             >  1:742160/1‑147 (MQ=255)
                                                                                                     aaTAGCGGGCCCATTGGCGGTTTTGGAAGGACAAA‑CCGGTCAGGGGAATCCACCCCGACCCGATAC‑AGCCCCT‑TCGCCTTTTCCGGTGCTGCCCCCCGCGTTCTTCCCTTTGCCCCCAGCCAATCTTCTTTCCCGCTCTCCCGCttggt                                             >  2:905847/1‑145 (MQ=255)
                                                                                                        aGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑AGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTgg                                         <  2:824693/150‑1 (MQ=255)
                                                                                                               ccATTGGCGGTTTTGGAAGATCATA‑ACGGGCCGGGGAAACCACCGCGACCAGATAA‑AGCACCT‑CCGCCTTTTCCGTTGGCGCCCCCCGCGGTCTGTCCTTTTTCCCCCGCAACTCTGCTTTCCCGCCCTCCGGCTGGTTTGGGGCACCGa                                  >  2:550376/1‑150 (MQ=255)
                                                                                                                      cGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑AGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATc                          >  1:107603/1‑151 (MQ=255)
                                                                                                                                  aTCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑AGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACa              <  1:537888/151‑1 (MQ=255)
                                                                                                                                   tCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑AGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACa              >  2:363746/1‑150 (MQ=255)
                                                                                                                                     aTA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑AGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACCat            >  1:464115/1‑150 (MQ=255)
                                                                                                                                      tA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑AGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAatat          >  2:762030/1‑151 (MQ=255)
                                                                                                                                         ccGGTCAGGGTAATCCACCGCGACCAGATAC‑AGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTc        <  1:804402/151‑1 (MQ=255)
                                                                                                                                           ggTCAGGGTAATCCACCGCGACCAGATAC‑AGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCt      <  2:925965/151‑1 (MQ=255)
                                                                                                                                               aGGGTAATCCACCGCGACCAGATAC‑AGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAcc  >  2:353522/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑AGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACC  >  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: ‑