Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,434,8610GT43.2% 129.4 / 42.4 81P259Q (CCA→CAA) truAtRNA pseudouridine(38‑40) synthase
Reads supporting (aligned to +/- strand):  ref base G (12/34);  new base T (35/0);  total (47/34)
Fisher's exact test for biased strand distribution p-value = 5.23e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAG‑CCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTG  >  NC_000913/2434713‑2435010
                                                                                                                                                    |                                                                                                                                                      
aaaTCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGa                                                                                                                                                      <  2:152760/151‑1 (MQ=255)
    cGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGa                                                                                                                                                   <  1:977110/150‑1 (MQ=255)
    cGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGa                                                                                                                                                   <  2:1001565/150‑1 (MQ=255)
    cGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGa                                                                                                                                                   <  2:1061162/150‑1 (MQ=255)
    cGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGa                                                                                                                                                   <  1:797164/150‑1 (MQ=255)
       tGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCa                                                                                                                                                >  2:577799/1‑150 (MQ=255)
            aaataaatCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGAACATACCgg                                                                                                                                          >  2:982042/1‑151 (MQ=255)
            aaataaatCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCgg                                                                                                                                          <  2:227366/151‑1 (MQ=255)
             aataaatCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGt                                                                                                                                         <  2:863426/151‑1 (MQ=255)
               taaatCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCgg                                                                                                                                          >  2:595240/1‑148 (MQ=255)
                 aatCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAgg                                                                                                                                     >  2:792797/1‑151 (MQ=255)
                   tCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATAACGGTCAggg                                                                                                                                    >  1:633889/1‑150 (MQ=255)
                       gTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATc                                                                                                                               <  2:1254031/151‑1 (MQ=255)
                         ccATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATcc                                                                                                                              <  1:867907/150‑1 (MQ=255)
                          cATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCCAACCGGAGAGGGGAATACa                                                                                                                             >  2:135813/1‑150 (MQ=255)
                          cATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCAc                                                                                                                            <  1:210109/151‑1 (MQ=255)
                          cATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCAc                                                                                                                            <  1:453026/151‑1 (MQ=255)
                             atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGGAAaccacc                                                                                                                           >  2:135540/1‑149 (MQ=255)
                               ttGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACcgcg                                                                                                                        <  1:42765/150‑1 (MQ=255)
                                     aTGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAACCGTCAGGGGAAACCACCGCGGCCAGa                                                                                                                  >  1:925903/1‑150 (MQ=255)
                                      tGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGCTTTTGGAGAACATACCGGGGGGGGTAATACACCCCCAcc                                                                                                                     >  2:1143748/1‑146 (MQ=255)
                                        ccTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCCGTCAGGGGAATCCACCCCGACCAGATa                                                                                                                >  1:886905/1‑149 (MQ=255)
                                        ccTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACa                                                                                                              <  2:246578/151‑1 (MQ=255)
                                        ccTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACa                                                                                                              <  2:676524/151‑1 (MQ=255)
                                        ccTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACa                                                                                                              <  2:882602/151‑1 (MQ=255)
                                        ccTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTTGGAGATCAAACCGGGCAGGGGAAACCACCCCGCCCCGATa                                                                                                                >  2:786166/1‑149 (MQ=255)
                                        ccTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTTGAAAATCCTACCCGGCAGGGGAAaccacccccaac                                                                                                                     >  2:581900/1‑138 (MQ=255)
                                          ttGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCAGGGGAATCCACCCCGGCCAGATa                                                                                                                >  1:957801/1‑147 (MQ=255)
                                           tGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAG‑cc                                                                                                          >  1:311182/1‑151 (MQ=255)
                                                tCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGAACATAACGGGCAGGGGAATTCACCCCGCCCCGAAAAAG‑CCCTTcc                                                                                                     >  2:502315/1‑151 (MQ=255)
                                                    cATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGAATATAACGGTCAGGGGAATCCACCCCGACCAGATACAA‑CCCTTCCCCtt                                                                                                 >  2:336369/1‑151 (MQ=255)
                                                    cATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGGCAGGGGAATCCACCGCGGCCCGATAAAA‑CCCCTCCGCtt                                                                                                 >  1:684953/1‑151 (MQ=255)
                                                         ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACATACCCGGCAGGGGAATCCACCGCGACCAGATACAG‑CCCTTCCCCtttt                                                                                               >  2:437647/1‑148 (MQ=255)
                                                         ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCCGTCCGGGGAATCCACCCCGACCAGATACAGCCCCTTCCCCtttt                                                                                               >  2:670393/1‑149 (MQ=255)
                                                         ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATAATAACGGGCAGGGGAAACAACCCAGGCCAGAAACAA‑CCCTTccccct                                                                                                 >  1:1261093/1‑142 (MQ=255)
                                                           cAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAG‑CCCTTCCGCTTTTGCCg                                                                                           <  1:117233/150‑1 (MQ=255)
                                                            aGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATTAAACCCGACAGGGGAATCTACCGCCAACCGATATAA‑CCCCTCCGCTTTTTccct                                                                                          >  2:607474/1‑148 (MQ=255)
                                                              gCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAG‑CCCTTCCGCTTTTGCCGTtgc                                                                                       <  2:1038198/151‑1 (MQ=255)
                                                                 gTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAG‑CCCTTCCGCTTTTGCCGTtgct                                                                                      <  1:162350/149‑1 (MQ=255)
                                                                   aaCTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGTCAGGGGAATCCACCCCGCCCAGAAACAG‑CCCTTCCCCTTTTTCCCTTTCTgccc                                                                                  >  1:422482/1‑150 (MQ=255)
                                                                    aCTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAACCGGCAGGGGAAACCACCGCGACCAGATACAG‑CCCCTCCGCTTTTTCCCTTGCTgccgc                                                                                 >  1:469806/1‑151 (MQ=255)
                                                                      tACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCAGGGGAAACCACCGCGACCCGATACAG‑CCCTTCCCCTTTTTCCGGTGCTgccgcc                                                                                >  2:1147780/1‑150 (MQ=255)
                                                                          ttGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAACGGGCAGGGGAATTCACCCCGACCAGATACAG‑ACCTTCCCCTTTTTCCGTTGCTtccccc                                                                                >  2:48357/1‑146 (MQ=255)
                                                                           tGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGTCAGGGGAAATCACCCCGACCAGATAAAG‑CCCTTCCCCTTTTGCCGCTGCTGCCGCCCGCGt                                                                           >  2:1073298/1‑150 (MQ=255)
                                                                             gCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCCGGCAGGGGAATCCCCCGCGGCCCGATACAG‑CCCCTCCGCCTTTGCCGGTGATGCCGCCCGCGGTCt                                                                        >  1:595633/1‑151 (MQ=255)
                                                                                  aGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGGAATCCACCCCGACCAGATACAG‑CCCCTCCGCTTTTTCCGTTGCTGCCGCCAGCGTTCTGTcc                                                                    >  1:669583/1‑150 (MQ=255)
                                                                                   gCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAG‑CCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCtt                                                                  <  1:502319/151‑1 (MQ=255)
                                                                                    ccTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAG‑CCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCttt                                                                 <  2:1067250/151‑1 (MQ=255)
                                                                                    ccTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAG‑CCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCttt                                                                 <  1:177593/151‑1 (MQ=255)
                                                                                        aaaTGTTCCGCGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCCTTCCGCCCAGGGGAATACCACGCGACCCACCACAG‑CCCTTTCCCTCTTTCCGGTCGTTCCGCCAGCGGTCTGGCttttttc                                                              >  2:780480/1‑147 (MQ=255)
                                                                                         aaTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGATGTTATACCCGTCAGGGGAATCCATCGCGGGCCGAAACAG‑CCCCTCCTCTTTTTCCGGTGGTGCCCCCAGCGGTCTGTCCttttgcgc                                                            >  2:1210324/1‑146 (MQ=255)
                                                                                             ttCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAACGGGCAGGGGAATCCACCCCGACCAGATACAG‑CCCTTCCCCTTTTTCCGTTGCTGCCCCCAGCGTTCTTTCCTTTtccgcc                                                           >  2:980184/1‑148 (MQ=255)
                                                                                               ccGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAG‑CCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCaa                                                      <  1:1196898/151‑1 (MQ=255)
                                                                                               ccGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAG‑CCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCa                                                       <  1:1176593/150‑1 (MQ=255)
                                                                                               ccGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAG‑CCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCa                                                       <  1:982053/150‑1 (MQ=255)
                                                                                                  aGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATCCACCGCGCCCAGATAAAG‑GCCTTCCGCTTTTTCCCGTTCTGCCGCCCGCGTTCTGTCCTTTTCCGCCCGCAActc                                                   >  2:679122/1‑151 (MQ=255)
                                                                                                  aGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAG‑CCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAActc                                                   <  2:217036/151‑1 (MQ=255)
                                                                                                    cATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAAACCACCGCGACCAGATACAA‑CCCCTCCGCTTTTGCCGTTGCTGCCGCCCGCGTTCTGTCCTTTTCCGCCAGCAActc                                                   >  1:330638/1‑149 (MQ=255)
                                                                                                     attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATAACCGGCAGGGGAATCCACCGCGACCCGATACAG‑CCCTTTCGCTTTTTCCGTTGCTGCCCCCCGCGGTCTGTCCTTTTCCGCCCGCCAccc                                                   >  2:856350/1‑146 (MQ=255)
                                                                                                       tattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAG‑CCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCt                                               <  2:941989/150‑1 (MQ=255)
                                                                                                        attCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATGATACCGGTCAGGGGAATCCACCGCGGCCAGATACAA‑CCCCTCCGCCTTTTCCGTTGCTGCCGCCCGCGGTCTGTCCTTTGCCGCCCGCAACTCTTCTa                                              >  2:596348/1‑150 (MQ=255)
                                                                                                          tctcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAG‑CCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATc                                            <  1:24549/150‑1 (MQ=255)
                                                                                                             cGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAA‑CCCTTCCGCTTTTTCCGTTGCTGCCGCCCGCGGTCTGCCCTTTTCCGCCCGCAACCCTGCTATCCCGc                                        >  2:191647/1‑151 (MQ=255)
                                                                                                                 aGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAG‑CCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGctct                                     <  2:28209/150‑1 (MQ=255)
                                                                                                                 aGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAG‑CCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCtctc                                    <  1:38282/151‑1 (MQ=255)
                                                                                                                   tCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCAAACCCGTCAGGGGAATCCACCGCGGCCCCATACAG‑CCCTTTCGCCTTTTCCGTTGCTGCCGCCAGGGTTTTGTCCTTTTTCGCCCGCAACGCTGCTCTCCAGCTCTccc                                  >  1:926675/1‑150 (MQ=255)
                                                                                                                     cGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCCGGCAGGGGAATCCACCGCGACCAGATACAG‑CCCCTCCGCTTTTGCCGTTGTTGCCGCCCGCGTTCTGTCCTTTTCCGCCCGCAACTCTTCTATCCAGCTCTcccg                                 >  1:325162/1‑148 (MQ=255)
                                                                                                                     cGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGAACATACCCGTCAGGGGAATCCACCCCGACCAGATACAA‑CCCTTCCCCTTTTTCCGTTGCTGCCGCCAGCGGTCTGTCCTTTTTCGCCCGCAACTCTGCTATTCAGCTCTcccg                                 >  1:812272/1‑148 (MQ=255)
                                                                                                                     cGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAG‑CCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGc                                <  1:733702/151‑1 (MQ=255)
                                                                                                                      gCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGTCAGGGGAATCCACCCCGACCCGATAAAG‑CCCTTCCCCTTTTTCCGGTGCTGCCCCCCGCGTTCTGTCCTTTTCCGCCAGCAACTCTGCTATCCAGCTCTCCGGc                                >  1:1256595/1‑150 (MQ=255)
                                                                                                                      gCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAG‑CCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCt                               <  1:123001/151‑1 (MQ=255)
                                                                                                                        cAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGAACATACCCGGCAGGGGAAACCACCCCGACCAGATAAAG‑CCCCTCCGCCTTTTCCGCTGCTGCCCCCAGCGGTCTGTCCTTTTCCCCCCGCAAATCTGCTATCCAGCTCTcccg                                 >  2:194160/1‑145 (MQ=255)
                                                                                                                           aaaTAGCGGGCCCATTGGCGGGTTTGGAAAATCATACCCGGCAAGGGAAATCACCCCCACCAGATACAG‑CCCTTCCCCTTTTGCCCTTGCTGCCCCCCGCGGTCTGTCCCTTTCCGCCCGCAACACTGCTATTCAGCTCTCCGGCTgg                             >  2:285034/1‑148 (MQ=255)
                                                                                                                               aGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAG‑CCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTgg                       <  2:330635/150‑1 (MQ=255)
                                                                                                                                gCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCCGGGGAATCCACCCCGACCAGATACAG‑CCCTTCCGCTTTTTCCGGTGGTGCCCCCCGCGGTCTTTCCTTTTCCGCCCGCAACTCTGCTATCCAGCTCTCCCGCTGGTTGTggg                      >  2:806993/1‑150 (MQ=255)
                                                                                                                                  gggCCCATTGGCGGTTTTTGGAGATCATAACGGGCAGGGGAATCCACCGCGACCAGATACAG‑CCCTTCCGCTTTTGCCGTTGCTGCCGCCCGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCa                    >  2:629311/1‑150 (MQ=255)
                                                                                                                                      ccATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAG‑CCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGAc               <  2:44661/151‑1 (MQ=255)
                                                                                                                                             cGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAG‑CCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATc        >  1:170616/1‑151 (MQ=255)
                                                                                                                                              ggTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAG‑CCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCa       <  2:929849/151‑1 (MQ=255)
                                                                                                                                                 tttGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAG‑CCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAgg     >  1:890742/1‑150 (MQ=255)
                                                                                                                                                   tGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAG‑CCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTg  >  1:1192458/1‑151 (MQ=255)
                                                                                                                                                    ggAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAG‑CCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTg  <  1:50095/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
AAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAG‑CCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTG  >  NC_000913/2434713‑2435010

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: ‑