Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,434,8610GT48.8% 75.1 / 56.9 80P259Q (CCA→CAA) truAtRNA pseudouridine(38‑40) synthase
Reads supporting (aligned to +/- strand):  ref base G (14/27);  new base T (39/0);  total (53/27)
Fisher's exact test for biased strand distribution p-value = 2.55e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTG  >  NC_000913/2434712‑2435010
                                                                                                                                                     |                                                                                                                                                       
aaaaTCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTg                                                                                                                                                         <  1:353277/150‑1 (MQ=255)
  aaTCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGaa                                                                                                                                                      <  1:243585/151‑1 (MQ=255)
     cGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGAt                                                                                                                                                   <  2:913103/151‑1 (MQ=255)
       aTGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGaa                                                                                                                                                      >  1:164360/1‑146 (MQ=255)
             aaataaatCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATAcccg                                                                                                                                           >  2:600188/1‑149 (MQ=255)
             aaataaatCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCCTAccc                                                                                                                                            >  2:318627/1‑149 (MQ=255)
              aataaatCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCgg                                                                                                                                           >  1:1197692/1‑150 (MQ=255)
              aataaatCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATAcc                                                                                                                                             >  2:578833/1‑148 (MQ=255)
                   atCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAggg                                                                                                                                     >  1:1135379/1‑151 (MQ=255)
                       ggTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAggg                                                                                                                                     >  1:1143992/1‑147 (MQ=255)
                         tCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGGAATcc                                                                                                                               >  2:346016/1‑151 (MQ=255)
                          ccATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGGAATcc                                                                                                                               >  2:810199/1‑150 (MQ=255)
                          ccATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCa                                                                                                                              <  2:1190844/151‑1 (MQ=255)
                            atatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATAACCGTCAGGGGAACCCAc                                                                                                                             >  2:226546/1‑150 (MQ=255)
                             tatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCAATGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATCCAcc                                                                                                                            >  1:188106/1‑150 (MQ=255)
                                ttGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACcgcg                                                                                                                         <  1:653043/150‑1 (MQ=255)
                                    cTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGTAATCCACCGCGAcc                                                                                                                      >  2:771600/1‑149 (MQ=255)
                                     tATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAACGGTCAGGGGAATTCACCGCGAcc                                                                                                                      >  1:787130/1‑148 (MQ=255)
                                     tATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCCGGCAGGGGAATCCACCCCGGcc                                                                                                                      >  2:688712/1‑148 (MQ=255)
                                        gCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATAc                                                                                                                <  2:380089/151‑1 (MQ=255)
                                         ccTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCCGTCAGGGGAATCCACCGCGACCCGATa                                                                                                                 >  1:1140594/1‑149 (MQ=255)
                                         ccTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGGCAGGGGAATCCACCCCGGCCAGaaa                                                                                                                 >  1:391029/1‑147 (MQ=255)
                                          cTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCCGTCAGGGGAATCCACCGCGACCCGATa                                                                                                                 >  1:267342/1‑148 (MQ=255)
                                              tACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGGCCGGGGAATCCACCGCGACCCGATACAGccc                                                                                                           >  1:133363/1‑150 (MQ=255)
                                              tACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCAAAACGGTCAGGGGAAACCACCGCGACCAGATACAGccc                                                                                                           >  2:497029/1‑150 (MQ=255)
                                               aCTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCtt                                                                                                         <  1:626676/151‑1 (MQ=255)
                                                   gTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACATAACCGTCAGGGGAAACCACCCCGGCCAGAAACAACCCCTcc                                                                                                      >  2:1193164/1‑149 (MQ=255)
                                                     cATCTTTCAGGCTGTAACTACGTTGGCTCCAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACATACCAGCCAAGGGAATCCTCCCCGACCCGATa                                                                                                                 >  1:847374/1‑137 (MQ=255)
                                                          ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCtttt                                                                                                <  2:706931/148‑1 (MQ=255)
                                                            cAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCAGGGGAATCCACCCCGACCAGATACAGCCCTTCCGCTTTTTCCGt                                                                                          >  2:248615/1‑151 (MQ=255)
                                                             aGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAAATCATACCGGGCAGGGTAATCCACCCCGACCCGATACAGCCCTTCCGCCTTTGCCg                                                                                           >  1:600932/1‑149 (MQ=255)
                                                             aGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGtt                                                                                         <  1:1168987/151‑1 (MQ=255)
                                                             aGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGtt                                                                                         <  1:628468/151‑1 (MQ=255)
                                                              ggCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTtg                                                                                        <  2:812158/151‑1 (MQ=255)
                                                                cTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTtg                                                                                        >  1:1001104/1‑149 (MQ=255)
                                                                      cTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCAGGGGAATCCACCGCCACCAGATAAAGCCCTTCCGCTTTTTCCCTtgctgc                                                                                    >  2:91/1‑147 (MQ=255)
                                                                      cTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCCGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTTCCGTTGCTgccgcc                                                                                >  2:1062270/1‑151 (MQ=255)
                                                                       tACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCAGGGGAAACCACCGCGACCAGATACAACCCCTCCGCTTTTTCCGTTGCTgccccc                                                                                >  2:992760/1‑150 (MQ=255)
                                                                            tGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCCTACCCGTCAGGGGAAACCACCCCGGCCAGATACAACCCCTTCCCTTTTTGCGTTTCCGCCCCCcgcggt                                                                          >  2:1124253/1‑149 (MQ=255)
                                                                             ggCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGtt                                                                          <  2:542703/150‑1 (MQ=255)
                                                                               cTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATTATACCGGGCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGTCGGTGGTGGCGCCAGCGGTCt                                                                        >  2:641870/1‑150 (MQ=255)
                                                                                 aCAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCCGGCAGGGTAAATCCCCGCGACCCGATACCGCCCTTCCGCTTTTGCCGTTTTTGCCCCCCGCGGTCTGt                                                                      >  1:1070312/1‑150 (MQ=255)
                                                                                  cAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGTAGATCATACCGGTCAGGGGAATTCACCGCCACCAGATACAACCCCTCCCCTTTTTCCGTTGCTGCCGCCCGCGTTCTTTcc                                                                    >  1:405973/1‑151 (MQ=255)
                                                                                    gCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATTATACCGGGCAGGGGAATCCACCCCGACCAGATAAAACCCCTCCCCTTTTGCCGTTGCCGCCCCCCGCGGTTTGTCCtt                                                                  >  1:76716/1‑151 (MQ=255)
                                                                                    gCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCtt                                                                  <  1:668041/151‑1 (MQ=255)
                                                                                    gCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCtt                                                                  <  2:1197816/151‑1 (MQ=255)
                                                                                     ccTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCCGGGTAATCCACCGCGACCAGATACAGCCCTTCCCCTTTTGCCGTTGCTGCCGCCCGCGCTCTGTCCttt                                                                 >  1:1061058/1‑151 (MQ=255)
                                                                                     ccTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCttt                                                                 <  2:600986/151‑1 (MQ=255)
                                                                                      cTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATCCACCCCGACCAGATACAACCCTTCCCCTTTTGCCGTTGCTGCCGCCAGCGTTCTTTCCttt                                                                 >  2:116128/1‑150 (MQ=255)
                                                                                      cTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAAACCACCCCGACCAGATACAGCCCTTCCGCTTTTTCCGTTGGTGCCGCCCGCGGTCTGTCCtttt                                                                >  2:88401/1‑150 (MQ=255)
                                                                                        gAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGTAGATTATACCCGGGAGGGTAATACACCCCGCACAGATACACCCCCTCCGCTTTTTCCGCTGATTCCCCCAGCGCTCTGTCCTTTgc                                                               >  2:1110688/1‑150 (MQ=255)
                                                                                        gAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCTTACCGGTCCGGGTAATGGAACCCCGGCAGATACAGCCCTTCCCATTTTTGCCGTGTTGCCCCCAGCGGTTTGTCCttt                                                                 >  1:112266/1‑148 (MQ=255)
                                                                                        gAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAACGGGCAGGGGAATCCCCCCCGGCCAGATACAGCCCTTCCCCTTTTTCCGGTTCTGCCCCCCGCGGTCTGTTCTTTgcc                                                              >  1:1119030/1‑151 (MQ=255)
                                                                                        gAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTg                                                                <  1:23555/149‑1 (MQ=255)
                                                                                         aaaTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGGCAGGGGAATCCAACCCGGCCAGATAAAACCCTTCCGCCTTTTCCGGTTCTGCCGCCCGCGTTCTTTCCTTTtcc                                                              >  1:788842/1‑150 (MQ=255)
                                                                                          aaTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGTTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGTCCAGATACAGCCCTTCCGCGTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTgccgc                                                            <  2:847456/151‑1 (MQ=255)
                                                                                              ttCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACATAACTGTCCGGGGAATCCACCGCGGCCAGATAAAGCCCTTCCGCCTTTGTCCGTTCTGGCGCCAGCGGTTTGTCCTTTtccgcc                                                           >  2:311878/1‑148 (MQ=255)
                                                                                                ccGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGTCAGGGGAATCCACCGCGACCCGATACAGCCCTTCCGCTTTTTCCGTTGCTGCCGCCAGCGGTCTTTCCTTTtcccccc                                                          >  1:76634/1‑146 (MQ=255)
                                                                                                      attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGc                                                <  1:1254101/151‑1 (MQ=255)
                                                                                                          ttCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATTATACCGGTCAGGGGAAACCACCCCGACCCGATACAGCCCCTCCGCTTTTTCCCTTGCTGCCCCCAGCGTTTTTTCCTTTTCCCCCCGCAACTCTCCTAt                                             >  1:781442/1‑150 (MQ=255)
                                                                                                          ttCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGGTCTGTCCTTTGCCGCCCGCAACTCTGCTAt                                             >  1:410260/1‑150 (MQ=255)
                                                                                                              agTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGGCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCCCCTTTTCCCTTGGTTCCCCCAGCGTTCTGTCCTTTGCCGCCCGCCACTCTGCTATCCCGc                                        >  2:664632/2‑151 (MQ=255)
                                                                                                                ttAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGGCAGGGGAATCCACCCCGACCAGATACAACCCTTCCGCTTTTTCCGTTGCTGCCGCCAGCGGTCTGTCCTTTTCCGCCAGCAACTCTGCTATCCAGctc                                      >  2:320762/1‑151 (MQ=255)
                                                                                                                ttAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGACAGAGCATACCGGTCAGGGTAACCACCCGCAACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCGGTCCTTTGCCGCCAGCAACTCTGCTATCCAGc                                        <  1:279776/149‑1 (MQ=255)
                                                                                                                   gTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTcc                                   <  2:370087/151‑1 (MQ=255)
                                                                                                                     ccGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCCGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTTCCGCCCGCAACTCTGCTATCCAGCTCTcc                                   >  2:114386/1‑149 (MQ=255)
                                                                                                                       gCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCt                               <  2:311293/151‑1 (MQ=255)
                                                                                                                        ccAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGGCGGGGAATCCACCGCCACCCGAAACAGCCCCTCCGCTTTTTCCGTTGCTGCCGCCAGCGTTCTTTCCCTTTCCGCCAACAACTCTGCTATCCAACCCTCCCGCTg                              >  2:222988/1‑151 (MQ=255)
                                                                                                                             aaTAGCGGGCCCATTGGCGGTTTTTGAAGATCATAACCGGCAGGGGAATCCACCCCGGCCAGATACAGCCCTTCCGCTTTTTCCGTTGCTGCCCCCCGCGTTCTTTCCTTTTCCGCCCGCCAATCCGCTATCCAGCTCTCCCGCTGGtt                           >  1:985414/1‑149 (MQ=255)
                                                                                                                             aaTAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTtgt                         <  1:157562/151‑1 (MQ=255)
                                                                                                                             aaTAGCGGGCCCATTGGCGGTTTTGGAAGATCATAACGGTCAGGGGAATCCACCGCGACCAGATACAACCCCTCCGCCTTTTCCCTTGCTTCCCCCCGCGGTCTTTCCTTTTCCGCCCGCCACTCTGCTATCCAGCTCTCCCGCTGGTtg                          >  2:719637/1‑150 (MQ=255)
                                                                                                                              aTAGCGGGCCCAATGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTtgtg                        <  2:188126/151‑1 (MQ=255)
                                                                                                                                    ggCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGTAATCCACCCCGACCAGATACAGCCCTTCCGCTTTTGCCCGTGCTGCCCCCAGCGGTCTGTCCTTTTCCGCCAGCCACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCAcc                  >  1:988638/1‑151 (MQ=255)
                                                                                                                                     gCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCg                 <  2:902342/151‑1 (MQ=255)
                                                                                                                                           tGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACtt             <  1:275703/149‑1 (MQ=255)
                                                                                                                                           tGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTc            <  1:791569/150‑1 (MQ=255)
                                                                                                                                             gCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCAt         <  2:308426/151‑1 (MQ=255)
                                                                                                                                               ggTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATc        >  1:893225/1‑150 (MQ=255)
                                                                                                                                               ggTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATc        >  1:593421/1‑150 (MQ=255)
                                                                                                                                               ggTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCa       <  2:1140712/151‑1 (MQ=255)
                                                                                                                                                gTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATc        >  2:729698/1‑149 (MQ=255)
                                                                                                                                                     ggAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTg  <  2:604740/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
AAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTG  >  NC_000913/2434712‑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: