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

TAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTG  >  NC_000913/2434711‑2435010
                                                                                                                                                      |                                                                                                                                                     
tAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGATTTg                                                                                                                                                       >  1:693343/1‑151 (MQ=255)
 aaaaTCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTgg                                                                                                                                                      <  1:1319908/151‑1 (MQ=255)
    aTCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAg                                                                                                                                                   <  2:815460/151‑1 (MQ=255)
          gAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATa                                                                                                                                             <  1:1256500/151‑1 (MQ=255)
            ggAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAcc                                                                                                                                           <  2:912125/151‑1 (MQ=255)
                ataaatCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGt                                                                                                                                        >  2:1030380/1‑150 (MQ=255)
                  aaatCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCa                                                                                                                                      <  2:10818/150‑1 (MQ=255)
                  aaatCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCa                                                                                                                                      <  2:1153633/150‑1 (MQ=255)
                  aaatCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAg                                                                                                                                     <  1:224362/151‑1 (MQ=255)
                   aatCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATATCGGTCAgg                                                                                                                                    >  2:130018/1‑151 (MQ=255)
                    atCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAggg                                                                                                                                   <  2:202727/151‑1 (MQ=255)
                     tCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGTGCATAACGGTCAcgg                                                                                                                                   >  2:1127085/1‑150 (MQ=255)
                      cAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTTTAAGAGCATACCCGGCAggg                                                                                                                                   >  1:235059/1‑149 (MQ=255)
                            cATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATTCa                                                                                                                            >  1:994132/1‑150 (MQ=255)
                             atatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCAAACCCGTCAGGGGAAaccacc                                                                                                                          >  1:120619/1‑151 (MQ=255)
                             atatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTGGAAGATCATAACCGGCAGGGGGATTCAcc                                                                                                                          >  2:1163266/1‑151 (MQ=255)
                               atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGTAGGTCTTAACGGTCCGGGGGATCcccccc                                                                                                                        >  1:282106/1‑149 (MQ=255)
                               atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACcgc                                                                                                                        <  2:75186/151‑1 (MQ=255)
                                 ttGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCCGGGGAAaccacccc                                                                                                                        >  1:214928/1‑147 (MQ=255)
                                 ttGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACATACCGGGCCGGGGAAAACACcg                                                                                                                         >  1:930457/1‑148 (MQ=255)
                                    ccTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCa                                                                                                                   <  2:202134/151‑1 (MQ=255)
                                     cTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATTAAACCGGGCAGGGGAATCCAccc                                                                                                                         >  1:332668/1‑143 (MQ=255)
                                      tATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAAATCATAACGGGCCGGGGaaaacaacccgac                                                                                                                     >  2:682347/1‑136 (MQ=255)
                                      tATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTTGAAGATCAAATCGGGCCGGGGAATCCACCGCGAcc                                                                                                                    >  2:902936/1‑148 (MQ=255)
                                       aTGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATTATAACGGTCAGGGGAATCCACCCCGccc                                                                                                                    >  2:1101290/1‑147 (MQ=255)
                                       aTGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGAt                                                                                                                <  2:233178/151‑1 (MQ=255)
                                         gCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTTGAGAACATTACTGTCAGGGGAAAACAACGCGACCAg                                                                                                                  >  1:473558/1‑147 (MQ=255)
                                            ttGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGTCAGGGGAATCCACCGCGGCCCGATACaac                                                                                                           >  2:52275/1‑149 (MQ=255)
                                            ttGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCCGGGGAAACCACCCCGACCAGATACa                                                                                                             >  2:659777/1‑149 (MQ=255)
                                               tACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCAAAACCGGCAGGGGAATCCACCCCGAccccaaacagacc                                                                                                         >  1:785035/1‑139 (MQ=255)
                                                aCTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCt                                                                                                        <  2:71524/150‑1 (MQ=255)
                                                 cTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCCTACCGGGCAGGGTGATCCACCCCGACCCCCTACACCCCTTc                                                                                                      >  1:379509/1‑151 (MQ=255)
                                                 cTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGGCAGGGGAATTCACCCCGGCCAGATACAAccc                                                                                                         >  2:239540/1‑148 (MQ=255)
                                                     tCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATAACGGTCAGGGGAATCCACCGCGACCCGAAACAGCCCCTcccc                                                                                                   >  1:1083987/1‑148 (MQ=255)
                                                          tttCAGGCTGCAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCAGGGGAATCCACCCCGACCAGAAACAGCCCTTCCCCtttt                                                                                               >  2:838154/1‑149 (MQ=255)
                                                           ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCAcc                                                                                                                          <  2:705738/121‑1 (MQ=255)
                                                            tCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGGGCGGGGAATCCACCCCGACCCGATACAACCCCTCCGCtttt                                                                                               >  1:227221/1‑147 (MQ=255)
                                                                   gTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTtgctg                                                                                     <  2:49747/150‑1 (MQ=255)
                                                                      aCTACGTTGGCTACAGCCTGAAATGTTCCGCGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTgccgc                                                                                 <  2:738420/151‑1 (MQ=255)
                                                                         aCGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAATCCACCGCGCACCGATAAAACCCTTCCCCTTTTTCCGTTTCTgccgcc                                                                                >  2:1237675/1‑149 (MQ=255)
                                                                         aCGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCa                                                                               <  1:959827/150‑1 (MQ=255)
                                                                         aCGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAg                                                                              <  2:695936/151‑1 (MQ=255)
                                                                          cGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGc                                                                             <  1:599859/151‑1 (MQ=255)
                                                                           gTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATAATAACGGGCAGGGGAAATCACCCCGCCCAGATAAAGCCCCTTCCCTTTTTCCGTTGCTTCCCCCcgcg                                                                            >  1:1188466/1‑151 (MQ=255)
                                                                           gTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGc                                                                             <  2:97101/150‑1 (MQ=255)
                                                                           gTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCg                                                                            <  2:1255653/151‑1 (MQ=255)
                                                                                  aCAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGTAATCATAACCGTCGGGGGTATCCCCCCCGACCAACTTAACCCCCTTCGCTTTTTCCCTTGCTTCTGCCCGCGGTTTGt                                                                      >  2:514533/1‑150 (MQ=255)
                                                                                    aGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCt                                                                   <  2:1094303/151‑1 (MQ=255)
                                                                                       cTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAAATCACCCCGGCCAGAAAAAACCCCTCCCCCTTTTCCGGTTCTgccgcgaaccggagagcac                                                                   >  2:639252/1‑134 (MQ=255)
                                                                                       cTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTgccg                                                                                  <  1:639233/133‑1 (MQ=255)
                                                                                          aaaTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGTGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTgccg                                                             <  2:1098613/151‑1 (MQ=255)
                                                                                                 ccGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGAGCAAAACGGGCAGGGGAAACCACCGCGACCAGATACAACCCCTCCGCTTTTTCCGTTGCTGCCGCCCGCGGTCTGTCCTTTtccgcc                                                           >  2:854490/1‑146 (MQ=255)
                                                                                                   gAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAACCGGCAGGGGAATCCACCGCCACCAGAAACAACCCTTCCGCCTTTTCCGTTGCTGCCCCCCGCGTTTTGTCCTTTGCCGCCCGCAAcc                                                    >  2:833375/1‑150 (MQ=255)
                                                                                                     gCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTTGAAGATCCTACCCGTCAGGGTAATCCACCGCGACCCAATACAGCCCTTCCCCCTTTTCCCTTTCTGCCGCCCACGTTCTGTCCCTTTCCGCCCGCaaaac                                                   >  1:439754/1‑147 (MQ=255)
                                                                                                      cATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCAAACCGGGCCGGGGAATACACCGCGGCCAGATACAGGCCTTCCGCTTTTGCCCTTGCTGCCGCCCGCGGTTTGTCCTTTGCCGCCCGCAACTCTg                                                 >  2:269530/1‑151 (MQ=255)
                                                                                                       attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAACCGTCAGGGGAATCCACCCCGACCCGATACAGCCCTTCCGCTTTTTCCGTTGCTGCCGCCAGCGTTTTGTCCTTTGCCGCCCGCCAccc                                                   >  2:1144745/1‑146 (MQ=255)
                                                                                                       attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAATATACCCGTCAGGGGAATCCCCCGCGCCCAGATACAACCCCTCCGCCTTTTCCGTTTCTGCCGCCCCCGTTCTGTCCTTTtcc                                                              >  1:877748/1‑137 (MQ=255)
                                                                                                            tctcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAACCCTTCCGCTTTTTCCGTTGCTGCCGCCAGCGTTCTGTCCTTTTCCGCCCGCAACCCTGCTATcc                                           >  2:901555/1‑151 (MQ=255)
                                                                                                             ctcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTTGTAGATTATACCCGGCGCGGGAATACCCCCCGACCCGATACAACCCTTCCGCCTTTTTCGGTGTTTCCCCCACCGTTCTTTCCTCTTTCGCCCGCAACTCTTCTATcc                                           >  1:969580/1‑150 (MQ=255)
                                                                                                              tcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATAACGGGCAGGGGAATCCACCGCGACCAGATACAACCCTTCCGCTTTTTCCGTTGCTGCCGCCAGCGTTCTGTCCTTTTTCGCCCGCCACTCTTCTAt                                             >  2:575986/1‑147 (MQ=255)
                                                                                                                  tAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGTAGAACATACCGGTTAGGGGGATACACCCCGACCACATACCGACCTTCCGCTTTTTCCGTTGGTGCCCCCCGCGGTCTGTCCTTTTGCGCCCGCCAATCTTCTAttccgc                                        >  1:938886/1‑143 (MQ=255)
                                                                                                                  tAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCAAAACGGTCAGGGTAAACCACCCCGGCCAGATACAACCCTTCCGCTTTTTCCGGTGCTGCCCCCCGCGGTCTGTCCCTTTCCCCCAGCAACCCCTCTATCCAGccc                                      >  2:775541/1‑148 (MQ=255)
                                                                                                                  tAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATAACCGTTAGGGTAATCCACCCCGACCCGATACAGCCCCTCCGCTTTTGGCGTTGCTGCCCCCAGCGGTCTGTCCTTTGCCCCCCGCAACTCTGCTTTCCa                                          >  2:1019248/1‑146 (MQ=255)
                                                                                                                  tAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGTAGAACATAACGGGGAGGGGAATTCCCCCCGACCAGATACAGCCCCTCCGCCTTTTGCGTTGCTGCCCCCCGCGGTCTGGCCTTTTCCGCCCGCAACCCCGCTTTCCAGct                                       >  2:472445/1‑149 (MQ=255)
                                                                                                                     tCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGGCAGGGGAAACCACCCCGACCCGAAAAAACCCCTCCCCTTTTTCCCTTGCTGCCGCCCGCGGTCTGTTCTTTTCCCCCAGCAACTCTGCTATCCAGCTCTcc                                   >  1:1150116/1‑150 (MQ=255)
                                                                                                                      ccGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAAATCACACCCGGCAGGGGAATTCACCCCGGCCCGATAAAAACCCTCCGCTTTTTTCGCTTCTGCCGCCCGCGTTTTGGCCCTTTTCCCCCGCAACTCTGCTTTCCCGctct                                     >  1:1224907/1‑147 (MQ=255)
                                                                                                                        gCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCt                               <  1:1113890/151‑1 (MQ=255)
                                                                                                                          cAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGGCAGGGGAATCCACCCCGGCCAGATACAACCCTTCCGCTTTTTCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAAATCTGCTATCCAGCTCTCCCGCTgg                             >  2:999153/1‑151 (MQ=255)
                                                                                                                           aGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCCCTTTTTCCGTTTCTGCCGCCAGCGTTCTGTCCTTTTCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTgg                             >  2:2336/1‑150 (MQ=255)
                                                                                                                            gAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGTCAGGGGAAACCCCCCCGGCCCGATAAAGCCCCTCCGCTTTTTCCGCTGCTGCCGCCCCCGTTCTGTCCTTTTCCGCCCGCAACTCTGCTATCCAGCCCTCCGGCTgg                             >  1:67661/1‑149 (MQ=255)
                                                                                                                               aTAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTtgtg                        <  1:47007/151‑1 (MQ=255)
                                                                                                                                        ccATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGAc               <  2:263824/151‑1 (MQ=255)
                                                                                                                                             ggCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTcc           <  2:760427/150‑1 (MQ=255)
                                                                                                                                               cGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATc        >  2:833501/1‑151 (MQ=255)
                                                                                                                                                ggTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCa       <  2:235156/151‑1 (MQ=255)
                                                                                                                                                 gTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAg      <  2:422636/151‑1 (MQ=255)
                                                                                                                                                    ttGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGc    <  1:891594/150‑1 (MQ=255)
                                                                                                                                                    ttGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGc    <  1:641217/150‑1 (MQ=255)
                                                                                                                                                      ggAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTg  <  2:973227/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTG  >  NC_000913/2434711‑2435010

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: