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

AAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTG  >  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)
                                           tGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAGcc                                                                                                          >  1:311182/1‑151 (MQ=255)
                                                tCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGAACATAACGGGCAGGGGAATTCACCCCGCCCCGAAAAAGCCCTTcc                                                                                                     >  2:502315/1‑151 (MQ=255)
                                                    cATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGAATATAACGGTCAGGGGAATCCACCCCGACCAGATACAACCCTTCCCCtt                                                                                                 >  2:336369/1‑151 (MQ=255)
                                                    cATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGGCAGGGGAATCCACCGCGGCCCGATAAAACCCCTCCGCtt                                                                                                 >  1:684953/1‑151 (MQ=255)
                                                         ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACATACCCGGCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCCCtttt                                                                                               >  2:437647/1‑148 (MQ=255)
                                                         ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCCGTCCGGGGAATCCACCCCGACCAGATACAGCCCCTTCCCCtttt                                                                                               >  2:670393/1‑149 (MQ=255)
                                                         ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATAATAACGGGCAGGGGAAACAACCCAGGCCAGAAACAACCCTTccccct                                                                                                 >  1:1261093/1‑142 (MQ=255)
                                                           cAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCg                                                                                           <  1:117233/150‑1 (MQ=255)
                                                            aGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATTAAACCCGACAGGGGAATCTACCGCCAACCGATATAACCCCTCCGCTTTTTccct                                                                                          >  2:607474/1‑148 (MQ=255)
                                                              gCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTtgc                                                                                       <  2:1038198/151‑1 (MQ=255)
                                                                 gTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTtgct                                                                                      <  1:162350/149‑1 (MQ=255)
                                                                   aaCTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGTCAGGGGAATCCACCCCGCCCAGAAACAGCCCTTCCCCTTTTTCCCTTTCTgccc                                                                                  >  1:422482/1‑150 (MQ=255)
                                                                    aCTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAACCGGCAGGGGAAACCACCGCGACCAGATACAGCCCCTCCGCTTTTTCCCTTGCTgccgc                                                                                 >  1:469806/1‑151 (MQ=255)
                                                                      tACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCAGGGGAAACCACCGCGACCCGATACAGCCCTTCCCCTTTTTCCGGTGCTgccgcc                                                                                >  2:1147780/1‑150 (MQ=255)
                                                                          ttGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAACGGGCAGGGGAATTCACCCCGACCAGATACAGACCTTCCCCTTTTTCCGTTGCTtccccc                                                                                >  2:48357/1‑146 (MQ=255)
                                                                           tGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGTCAGGGGAAATCACCCCGACCAGATAAAGCCCTTCCCCTTTTGCCGCTGCTGCCGCCCGCGt                                                                           >  2:1073298/1‑150 (MQ=255)
                                                                             gCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCCGGCAGGGGAATCCCCCGCGGCCCGATACAGCCCCTCCGCCTTTGCCGGTGATGCCGCCCGCGGTCt                                                                        >  1:595633/1‑151 (MQ=255)
                                                                                  aGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGGAATCCACCCCGACCAGATACAGCCCCTCCGCTTTTTCCGTTGCTGCCGCCAGCGTTCTGTcc                                                                    >  1:669583/1‑150 (MQ=255)
                                                                                   gCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCtt                                                                  <  1:502319/151‑1 (MQ=255)
                                                                                    ccTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCttt                                                                 <  2:1067250/151‑1 (MQ=255)
                                                                                    ccTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCttt                                                                 <  1:177593/151‑1 (MQ=255)
                                                                                        aaaTGTTCCGCGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCCTTCCGCCCAGGGGAATACCACGCGACCCACCACAGCCCTTTCCCTCTTTCCGGTCGTTCCGCCAGCGGTCTGGCttttttc                                                              >  2:780480/1‑147 (MQ=255)
                                                                                         aaTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGATGTTATACCCGTCAGGGGAATCCATCGCGGGCCGAAACAGCCCCTCCTCTTTTTCCGGTGGTGCCCCCAGCGGTCTGTCCttttgcgc                                                            >  2:1210324/1‑146 (MQ=255)
                                                                                             ttCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAACGGGCAGGGGAATCCACCCCGACCAGATACAGCCCTTCCCCTTTTTCCGTTGCTGCCCCCAGCGTTCTTTCCTTTtccgcc                                                           >  2:980184/1‑148 (MQ=255)
                                                                                               ccGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCaa                                                      <  1:1196898/151‑1 (MQ=255)
                                                                                               ccGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCa                                                       <  1:1176593/150‑1 (MQ=255)
                                                                                               ccGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCa                                                       <  1:982053/150‑1 (MQ=255)
                                                                                                  aGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATCCACCGCGCCCAGATAAAGGCCTTCCGCTTTTTCCCGTTCTGCCGCCCGCGTTCTGTCCTTTTCCGCCCGCAActc                                                   >  2:679122/1‑151 (MQ=255)
                                                                                                  aGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAActc                                                   <  2:217036/151‑1 (MQ=255)
                                                                                                    cATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAAACCACCGCGACCAGATACAACCCCTCCGCTTTTGCCGTTGCTGCCGCCCGCGTTCTGTCCTTTTCCGCCAGCAActc                                                   >  1:330638/1‑149 (MQ=255)
                                                                                                     attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATAACCGGCAGGGGAATCCACCGCGACCCGATACAGCCCTTTCGCTTTTTCCGTTGCTGCCCCCCGCGGTCTGTCCTTTTCCGCCCGCCAccc                                                   >  2:856350/1‑146 (MQ=255)
                                                                                                       tattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCt                                               <  2:941989/150‑1 (MQ=255)
                                                                                                        attCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATGATACCGGTCAGGGGAATCCACCGCGGCCAGATACAACCCCTCCGCCTTTTCCGTTGCTGCCGCCCGCGGTCTGTCCTTTGCCGCCCGCAACTCTTCTa                                              >  2:596348/1‑150 (MQ=255)
                                                                                                          tctcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATc                                            <  1:24549/150‑1 (MQ=255)
                                                                                                             cGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAACCCTTCCGCTTTTTCCGTTGCTGCCGCCCGCGGTCTGCCCTTTTCCGCCCGCAACCCTGCTATCCCGc                                        >  2:191647/1‑151 (MQ=255)
                                                                                                                 aGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGctct                                     <  2:28209/150‑1 (MQ=255)
                                                                                                                 aGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCtctc                                    <  1:38282/151‑1 (MQ=255)
                                                                                                                   tCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCAAACCCGTCAGGGGAATCCACCGCGGCCCCATACAGCCCTTTCGCCTTTTCCGTTGCTGCCGCCAGGGTTTTGTCCTTTTTCGCCCGCAACGCTGCTCTCCAGCTCTccc                                  >  1:926675/1‑150 (MQ=255)
                                                                                                                     cGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCCGGCAGGGGAATCCACCGCGACCAGATACAGCCCCTCCGCTTTTGCCGTTGTTGCCGCCCGCGTTCTGTCCTTTTCCGCCCGCAACTCTTCTATCCAGCTCTcccg                                 >  1:325162/1‑148 (MQ=255)
                                                                                                                     cGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGAACATACCCGTCAGGGGAATCCACCCCGACCAGATACAACCCTTCCCCTTTTTCCGTTGCTGCCGCCAGCGGTCTGTCCTTTTTCGCCCGCAACTCTGCTATTCAGCTCTcccg                                 >  1:812272/1‑148 (MQ=255)
                                                                                                                     cGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGc                                <  1:733702/151‑1 (MQ=255)
                                                                                                                      gCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGTCAGGGGAATCCACCCCGACCCGATAAAGCCCTTCCCCTTTTTCCGGTGCTGCCCCCCGCGTTCTGTCCTTTTCCGCCAGCAACTCTGCTATCCAGCTCTCCGGc                                >  1:1256595/1‑150 (MQ=255)
                                                                                                                      gCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCt                               <  1:123001/151‑1 (MQ=255)
                                                                                                                        cAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGAACATACCCGGCAGGGGAAACCACCCCGACCAGATAAAGCCCCTCCGCCTTTTCCGCTGCTGCCCCCAGCGGTCTGTCCTTTTCCCCCCGCAAATCTGCTATCCAGCTCTcccg                                 >  2:194160/1‑145 (MQ=255)
                                                                                                                           aaaTAGCGGGCCCATTGGCGGGTTTGGAAAATCATACCCGGCAAGGGAAATCACCCCCACCAGATACAGCCCTTCCCCTTTTGCCCTTGCTGCCCCCCGCGGTCTGTCCCTTTCCGCCCGCAACACTGCTATTCAGCTCTCCGGCTgg                             >  2:285034/1‑148 (MQ=255)
                                                                                                                               aGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTgg                       <  2:330635/150‑1 (MQ=255)
                                                                                                                                gCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCCGGGGAATCCACCCCGACCAGATACAGCCCTTCCGCTTTTTCCGGTGGTGCCCCCCGCGGTCTTTCCTTTTCCGCCCGCAACTCTGCTATCCAGCTCTCCCGCTGGTTGTggg                      >  2:806993/1‑150 (MQ=255)
                                                                                                                                  gggCCCATTGGCGGTTTTTGGAGATCATAACGGGCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCCGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCa                    >  2:629311/1‑150 (MQ=255)
                                                                                                                                      ccATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGAc               <  2:44661/151‑1 (MQ=255)
                                                                                                                                             cGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATc        >  1:170616/1‑151 (MQ=255)
                                                                                                                                              ggTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCa       <  2:929849/151‑1 (MQ=255)
                                                                                                                                                 tttGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAgg     >  1:890742/1‑150 (MQ=255)
                                                                                                                                                   tGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTg  >  1:1192458/1‑151 (MQ=255)
                                                                                                                                                    ggAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTg  <  1:50095/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
AAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTG  >  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: