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

AAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCA  >  NC_000913/2434713‑2435005
                                                                                                                                                    |                                                                                                                                                 
aaaTCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGa                                                                                                                                                 <  2:669470/151‑1 (MQ=255)
 aaTCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGaa                                                                                                                                                <  1:388386/151‑1 (MQ=255)
     gATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGa                                                                                                                                              >  1:216555/1‑149 (MQ=255)
     gATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATc                                                                                                                                            <  2:1043052/151‑1 (MQ=255)
      aTGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTgg                                                                                                                                                  >  1:697700/1‑144 (MQ=255)
      aTGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCa                                                                                                                                           <  1:463033/151‑1 (MQ=255)
       tGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGGTCa                                                                                                                                           >  2:911544/1‑150 (MQ=255)
        gAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCAt                                                                                                                                          <  1:416949/150‑1 (MQ=255)
           gAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCg                                                                                                                                     <  1:975568/151‑1 (MQ=255)
           gAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCg                                                                                                                                     <  1:579012/151‑1 (MQ=255)
           gAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCg                                                                                                                                     <  1:551375/151‑1 (MQ=255)
             aataaatCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCgg                                                                                                                                    >  2:500603/1‑150 (MQ=255)
               taaatCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTc                                                                                                                                  >  1:321181/1‑150 (MQ=255)
                  atCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCCGTCAggg                                                                                                                              >  2:218351/1‑151 (MQ=255)
                   tCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACATAACCGTCCggg                                                                                                                              >  2:1058420/1‑150 (MQ=255)
                    cAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGGCAgggg                                                                                                                             >  1:387253/1‑149 (MQ=255)
                     aGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAggg                                                                                                                              >  2:925536/1‑148 (MQ=255)
                        tCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGTAATGATACCCGCCAGGGtgcggc                                                                                                                        >  1:101944/1‑146 (MQ=255)
                            tatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAAaccacc                                                                                                                     >  2:786712/1‑150 (MQ=255)
                             atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGTAAAGCATCACCGCCAGGGTAATACAccc                                                                                                                    >  1:195658/1‑149 (MQ=255)
                               ttGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACcgcg                                                                                                                  <  1:669793/150‑1 (MQ=255)
                                     aTGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACATACCGGGCAGGGGAATACACCCCCAcc                                                                                                               >  1:877439/1‑147 (MQ=255)
                                     aTGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAAATCATACCCGTCAGGGAAAACCaaccccacc                                                                                                               >  1:484077/1‑139 (MQ=255)
                                        ccTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATTATAACCGGCCGGGGGATCCACCCCGACCa                                                                                                              >  2:928584/1‑145 (MQ=255)
                                            gTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATAACGGTCCGGGGAATCCACCGCGACCAGATACAAccc                                                                                                    >  1:989904/1‑151 (MQ=255)
                                            gTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGccc                                                                                                    <  1:1088114/151‑1 (MQ=255)
                                             tACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGAACATACCGGGCAGGGGAAACCACCCCGACCAGATAAAAccc                                                                                                    >  1:1070613/1‑150 (MQ=255)
                                                         ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTg                                                                                         <  1:1072629/149‑1 (MQ=255)
                                                           cAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCCTACCGGGCAGGGGAATCCACCGCGCCCAGATACAGCCCTTCCGCTTTTGcc                                                                                       >  1:262949/1‑149 (MQ=255)
                                                           cAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAACCCTTCCGCTTTTTcc                                                                                       >  1:209932/1‑149 (MQ=255)
                                                            aGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGtt                                                                                    <  1:137280/151‑1 (MQ=255)
                                                             ggCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTtg                                                                                   <  1:219096/151‑1 (MQ=255)
                                                              gCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTtgc                                                                                  <  2:1047362/151‑1 (MQ=255)
                                                                 gTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTtgct                                                                                 <  1:128393/149‑1 (MQ=255)
                                                                      tACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAAACCACCGCGACCAGATACAGCCCTTCCGCTTTTTCCGGTGCTtcccccc                                                                          >  2:107522/1‑150 (MQ=255)
                                                                         gTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATCCACCGCGACCCGATACAGCCCTTCCGCCTTTTCCGGTGCTgccgcc                                                                           >  1:398545/1‑147 (MQ=255)
                                                                         gTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCCGTCAGGGGAATACACCGCGGCCCGAAACAACCCCTCCGCTTTTTTCGTTGCTgccgcc                                                                           >  1:579402/1‑147 (MQ=255)
                                                                            ggCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATAACGGGCAGGGGAAACCACCCCGACCAGAAACAGCCCTTCCGCTTTTTCCGGTGCTGCCGCCAGCggt                                                                     >  1:1052757/1‑148 (MQ=255)
                                                                            ggCTACAGCCTGAAATGTTCCGAACATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATAACGGGCAGGGGAATCCACCGCGACCAGATACAACCCTTCCCCTTTTTCCGTTGCTgccccc                                                                           >  1:783795/1‑144 (MQ=255)
                                                                                 cAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATTCACCCCGACCAGATACAGCCCTTCCGCTTTTTCCCTTGCTGCCGCCAGCGTTCTGTcc                                                               >  1:569620/1‑151 (MQ=255)
                                                                                 cAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTcc                                                               <  1:501845/151‑1 (MQ=255)
                                                                                  aGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCt                                                              <  1:500608/151‑1 (MQ=255)
                                                                                   gCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCtt                                                             <  2:580994/151‑1 (MQ=255)
                                                                                       gAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATCCACCCCGACCAGATACAGCCCTTCCGCCTTTGTCGTTTCTGCCGCCCGCGTTTTGTCCctt                                                            >  2:61844/1‑148 (MQ=255)
                                                                                        aaaTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAATCCACCCCGACCCGATACAGCCCTTCCGCCTTTTCCGTTGCTGCCGCCCGCGGTCTTTCCTTTtcc                                                         >  2:1067511/1‑150 (MQ=255)
                                                                                            gTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAg                                                    <  1:150433/151‑1 (MQ=255)
                                                                                                cGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGGAATCCACCGCGGCCAGATACAGCCCCTCCGCTTTTTCCGGTGGTGCCGCCAGCGGTCTGTTCTTTGCCGCCCGCaa                                                 >  2:261229/1‑150 (MQ=255)
                                                                                                  aGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCCGGCAGGGGAATCCACCGCGACCAGATACAACCCTTCCGCTTTTTCCGTTTCTGCCGCCAGCGGTCTTTCCTTTTCCGCCAGCaa                                                 >  1:966566/1‑148 (MQ=255)
                                                                                                   gCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGTCAGGGGAAACCACCCCGGCCCGATACAACCCTTCCGCCTTTTCCGTTGCTTCCGCCCGCGGTCTGTCCTTTTCCGCCAGCaa                                                 >  1:422000/1‑147 (MQ=255)
                                                                                                   gCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAActct                                             <  2:195657/151‑1 (MQ=255)
                                                                                                    cATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCCGGCAGGGGAATCCACCCCGACCAGATAAAGCCCCTCCCCTTTTTCCCTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCCActc                                              >  1:1076379/1‑149 (MQ=255)
                                                                                                        attCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCCGGGGAATCCACCCCGACCAGATAAAGCCCCTTCCCTTTTTCCGTTTCTTCCCCCCGCGGTCTTTCCTTTTCCGCCCGCaa                                                 >  2:730262/1‑142 (MQ=255)
                                                                                                          tctcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGAGCATACCCGGGCGGGGGATCCCACCCGACCAGATACAGCCCTTCCTCTTTTGTCGGTGCTGGCGCCAGCGTTCTGTCCTTTTCCGCCCCCCACCCTGCTATcc                                      >  2:459739/1‑151 (MQ=255)
                                                                                                          tctcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGGCAGGGGAATCCCCCGCGACCAGATACAGCCCTTCCCCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTTCCGCCCGCAACTCTGCTATcc                                      >  1:117950/1‑151 (MQ=255)
                                                                                                           ctcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCCGGCAGGGGAATCCACCGCGGCCAGATACAGCCCTTTCCCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTTCCGCCCGCAACTCTGCTATcc                                      >  1:211653/1‑150 (MQ=255)
                                                                                                           ctcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGTCAGGGTAATCCACCCCGACCAGATACAGCCCTTCCCCTTTTTCCGTTTCTGCCCCCAGCGTTCTGTCCTTTGCCGCCCGCAACCCTGCTATcc                                      >  1:1017887/1‑150 (MQ=255)
                                                                                                             cGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATcc                                      <  2:1065785/148‑1 (MQ=255)
                                                                                                              gTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGct                                  <  1:837598/151‑1 (MQ=255)
                                                                                                                tAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGGCAGGGGAAACCACCCCCACCCGAAACAGCCCCTCCCCTTTTTCCGGTGCTTCCGCCCGCGGTCTTTCCCTTTCCGCCCGCAACTCTGCTTTTCAGctc                                 >  2:92069/1‑150 (MQ=255)
                                                                                                                 aGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCCGGGTAATCCACCCCGGCCCGATACAGCCCCTCCCCCTTTTGCGGTGCTGCCCCCCGCCGTCTTTCCCTTTCCGCCAGCCAATCTGCTATTCAGc                                   >  2:1131033/1‑147 (MQ=255)
                                                                                                                    ccGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCgg                            <  2:216554/151‑1 (MQ=255)
                                                                                                                    ccGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCg                             <  2:946914/150‑1 (MQ=255)
                                                                                                                     cGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAATCCACCGCGCCCAGATACAGCCCCTCCCCTTTTTCCCTTTCTTCCCCCCGCCTTCTGTCCTTTGCCGCCCGCCACTCTGCTATCCACCTCTcc                              >  2:149584/1‑148 (MQ=255)
                                                                                                                     cGCCAGAAATAGCGGGCCCATTGGCGGGTTTTGGAGATCATACCGGGCAGGGGAATCCACCGCGACCAGAGACAGCCCCTCCCCTTTTTCCGTTGCTGCCCCCAGCCTTCTGTCCCTTTCCGCCCACAAATCTTCTATCCACCcctc                               >  2:808885/1‑147 (MQ=255)
                                                                                                                      gCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATAACGGTCAGGGGAATCCACCCCGACCAGATACAGCCCCTCCGCTTTTTCCGTTTCTGCCGCCCCCGGTCTTTCCTTTTCCCCCAGCCACTCTTCTTTCCAGCTCTcc                              >  2:54010/1‑147 (MQ=255)
                                                                                                                       ccAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTg                         <  1:503201/151‑1 (MQ=255)
                                                                                                                         aGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCCTCCGCCTTTTCCGTTGCTGCCCCCAGCGTTCTGTCCTTTTCCGCCCACAACTCTTCTATCCCGCTCTcc                              >  2:216814/1‑144 (MQ=255)
                                                                                                                          gAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGGCAGGGGAAACCACCCCGACCCGATACAGCCCTTCCGCTTTTTCCGGTGCTGCCGCCCGCGGTCTGTCCTTTTCCGCCCGCAACTCTGCTATCCAGCCCTCCCGCTgggt                      >  1:622329/1‑149 (MQ=255)
                                                                                                                             aTAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATCCCCCGCGAACAGATACAACCCTTCCCCCTTTTCCCTTGCTGCCCCCAGCGTTCTGTCCTTTTCCGCCAGCCAATCTCCTATTCAGCTCTCCCGCTGGgtgt                    >  1:644724/1‑150 (MQ=255)
                                                                                                                              tAGCGGGCCCATTGGCGGTTTTTGGAGAATATAACCGGCAGGGGAATCCACCGCGGCCAGATACAACCCCTCCGCTTTTTCCGGTGCTGCCGCCCGCGGTCTGTCCCTTTCCGCCCGCCACTCTGCTAATCCGCTCTCCGGCTGGGTGTgg                  >  2:418505/1‑151 (MQ=255)
                                                                                                                              tAGCGGGCCCATTGGCGGTTTTTGAAGATCAAACCCGGCAGGGTAATCCACCGCGACCCGATACAGCCCCTCCGCCTTTTCCGTTGCTGCCGCCAGCGCTCTTTCCTTTTCCGCCAGCAACTCCGCTATCCCGCTCTCCCGCTGGGTGTgg                  >  2:915292/1‑151 (MQ=255)
                                                                                                                               aGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTggg                 <  2:759851/151‑1 (MQ=255)
                                                                                                                                 cGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCa               <  2:375637/151‑1 (MQ=255)
                                                                                                                                  gggCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCAc              <  2:440065/151‑1 (MQ=255)
                                                                                                                                      ccATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGAc          <  2:426903/151‑1 (MQ=255)
                                                                                                                                            gCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCAt    >  2:993311/1‑151 (MQ=255)
                                                                                                                                              ggTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCa  <  2:590859/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCA  >  NC_000913/2434713‑2435005

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: