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

CGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTG  >  NC_000913/2434717‑2435010
                                                                                                                                                 |                                                                                                                                                       
cGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGa                                                                                                                                                    <  2:43484/150‑1 (MQ=255)
 gATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGaa                                                                                                                                                      >  2:637607/1‑147 (MQ=255)
   tGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCa                                                                                                                                                 >  1:409166/1‑150 (MQ=255)
    gAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCAt                                                                                                                                                >  1:676263/1‑150 (MQ=255)
     aGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGgagaatcta                                                                                                                                               >  1:160681/1‑145 (MQ=255)
     aGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAc                                                                                                                                              <  1:91208/151‑1 (MQ=255)
       gAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAAgttatacc                                                                                                                                              >  1:26508/1‑143 (MQ=255)
       gAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATaac                                                                                                                                             >  1:649041/1‑148 (MQ=255)
          ataaatCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTc                                                                                                                                         >  2:1052897/1‑151 (MQ=255)
           taaatCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAGCATACCCGt                                                                                                                                          >  1:537025/1‑149 (MQ=255)
            aaatCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGgagaatattaccg                                                                                                                                           >  1:605320/1‑138 (MQ=255)
            aaatCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCa                                                                                                                                        <  1:645986/150‑1 (MQ=255)
              atCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAggg                                                                                                                                     <  1:286845/151‑1 (MQ=255)
                cAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAggg                                                                                                                                     >  1:981583/1‑149 (MQ=255)
                 aGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGAATATACCCGTCAggg                                                                                                                                     >  2:1210177/1‑148 (MQ=255)
                   gTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCCGGGGaa                                                                                                                                  >  1:359593/1‑149 (MQ=255)
                   gTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTGGAAGATCATACCGGTCAGGGGaa                                                                                                                                  >  2:875738/1‑149 (MQ=255)
                        tatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGGTTCTAACGGGCCCGGTAATTCAcc                                                                                                                            >  1:631182/1‑150 (MQ=255)
                        tatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGGCAGGGGAAACCAc                                                                                                                             >  1:1037875/1‑149 (MQ=255)
                        tatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATAATAACCGGCAGGGGAAaccacc                                                                                                                            >  1:238918/1‑150 (MQ=255)
                         atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGGAATCCACcgc                                                                                                                          >  2:1223738/1‑151 (MQ=255)
                            tGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATAACGGTCAGGGTAATCCAcccc                                                                                                                          >  1:582459/1‑146 (MQ=255)
                             gCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGa                                                                                                                       <  1:497427/149‑1 (MQ=255)
                             gCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGAc                                                                                                                      <  1:551063/150‑1 (MQ=255)
                                 aTGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAg                                                                                                                   <  1:1007028/149‑1 (MQ=255)
                                  tGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCCGGGGAATCCACCGCGAcc                                                                                                                     >  1:1149038/1‑146 (MQ=255)
                                    ccTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACa                                                                                                              <  1:753573/151‑1 (MQ=255)
                                       tGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAg                                                                                                             <  2:837186/149‑1 (MQ=255)
                                         tACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGGAATCCACCCCGACCAGATACAGccc                                                                                                         >  2:157118/1‑150 (MQ=255)
                                         tACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCAAACCGGGCAGGGGAATCCACCGCGACCAGATACAGccc                                                                                                         >  1:895472/1‑150 (MQ=255)
                                          aCTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCCGGGGAATCCACCGCGACCAAAAACAGccc                                                                                                         >  2:502360/1‑149 (MQ=255)
                                          aCTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGCTTTGGAAAAACAAACCGGTTAGGGGAATTCACCCCGAACAAATACAGCCCt                                                                                                        >  2:629583/1‑150 (MQ=255)
                                              gTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCg                                                                                                    >  2:658734/1‑150 (MQ=255)
                                                   cTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCtttt                                                                                               >  2:501488/1‑150 (MQ=255)
                                                     ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCtttt                                                                                               <  2:896687/148‑1 (MQ=255)
                                                       cAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTTTAAGATCATACCGGGCCAGGTAATCCACCGCGGCCCGAAACAGACCCTCCCCCttttc                                                                                              >  1:723945/1‑147 (MQ=255)
                                                        aGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGTGCGGGGAATACACCCCGACCAGAAACAACCCCTCCGCTTTTTcc                                                                                            >  2:435522/1‑148 (MQ=255)
                                                        aGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGGCAGGGGAATCCACCCCGACCAGATACAGCCCCTCCCCCTTTTCCGtt                                                                                         >  1:161755/1‑151 (MQ=255)
                                                              tAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAAATCACCCCGACCAGATACAACCCTTCCCCTTTTTCCGCTGCTgcc                                                                                   >  2:941457/1‑151 (MQ=255)
                                                               aaCTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTgccg                                                                                  <  1:1126039/151‑1 (MQ=255)
                                                                 cTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATACACCCCGACCCGATACACGCCTTCCGCCTTTTCCGTTGCTgccgcc                                                                                >  1:547133/1‑151 (MQ=255)
                                                                   aCGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCAGGGGAATTCACCCCGACCAGATAAAACCCTTCCCCTTTTTCCGTTGCTgccgcc                                                                                >  1:791646/1‑149 (MQ=255)
                                                                   aCGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCa                                                                               <  2:537026/150‑1 (MQ=255)
                                                                      ttGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGCTTTGGGAAATTATACCCGACAGGGGAAATCACCCCCGCCAGAATCAACCCCTCCCCTTTTGCCGTTTCCGCCCCCAGCGt                                                                           >  2:232743/1‑151 (MQ=255)
                                                                        ggCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGTAGATTATACCGGTCCGGGGAAAACACCCCGACAAGATACAACGCTTCCCGTTTTTTCGTTGCTgccgcccccg                                                                            >  1:787800/1‑144 (MQ=255)
                                                                        ggCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTc                                                                         <  1:942662/151‑1 (MQ=255)
                                                                          cTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGTCAGGGGAATCCACCCCGAACCGATACAACCCTTCCTCTTTTGCCGGTGCTGCCGCCAGccgtc                                                                         >  1:227637/1‑145 (MQ=255)
                                                                            aCAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGt                                                                      <  1:1093823/150‑1 (MQ=255)
                                                                                ccTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCCGTCAGGGGAATCCACCGCGACCAGATACAACCCTTCCGCTTTTTCCGTTGCTGCCGCCAGCGGTCTGTCCttt                                                                 >  2:269452/1‑151 (MQ=255)
                                                                                ccTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCttt                                                                 <  2:738498/151‑1 (MQ=255)
                                                                                  tGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAAACCACCGCGACCAGATACAGCCCCTCCGCCTTTGCCGTTGCTGCCGCCAGCGGTCTGTCCttt                                                                 >  1:552207/1‑149 (MQ=255)
                                                                                  tGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAACCGACCCGGGAATCCACCCCCACCAGATACAGCCCTTCCGCCTTTTCCCCTGCTTCCGCCACCGTTCTGTCCctt                                                                 >  2:685653/1‑149 (MQ=255)
                                                                                   gAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGGCCGGGGAATCCACCCCGGCCAGATACAACCCCTCCCCTTTTTCCGTTGCTGCCGCCAGCGTTCTGTCCttt                                                                 >  1:836649/1‑148 (MQ=255)
                                                                                         ttCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCCGTCAGGGGAATCCACCGCGACCAGATACAGCCCCTCCGCTTTTGCCGTTGCTGCCGCCAGCGGTCTTTCCTTTGccgcccgc                                                        >  1:1068074/1‑151 (MQ=255)
                                                                                         ttCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCCGGCAGGGGAATTCACCCCGACCAGATACAGCCCTTCCGCTTTTTCCGTTGCTGCCGCCCGCGGTTTTTCCCTTgccccc                                                           >  1:652501/1‑148 (MQ=255)
                                                                                           ccGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCCGTCAGGGGAATCCACCCCGCCCAGATACAGCCCCTTCCCTTTTTCCCTTGCTGCCGCCCACGGTCTGTCCTTTGCCCCCAGCa                                                       >  1:853793/1‑150 (MQ=255)
                                                                                            cGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGAAAAAACCCGTCCGGGGTAAAAACCTCGCCCCGGAAAACGCCCTTCTCTTTTTCCGTTGCCGCGCCCCGCGTTCTGTTCtt                                                                  >  2:422692/1‑138 (MQ=255)
                                                                                             gAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGTAATACACCGCGACCAGATACAACCCTTCCGCCTTTGCCGGTGCTGCCCCCCGCGGTCTGTCCTTTTCCGCCCGCaa                                                      >  2:202638/1‑149 (MQ=255)
                                                                                                cATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATAATACCGGGCAGGGGAATCCACCCCGACCAGATACAGCCCTTCCCCCTTTTCCGTTGCTGCCGCCAGCGGTCTGTCCCTTTCCCCCCGCAActc                                                   >  2:1162727/1‑149 (MQ=255)
                                                                                                 attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAAATTATACCGGCCCGGGGAAACCCCCACGACCCGATAAAGCCCTTTCGCCTTTGCCGGTGCGGCCGCCCGCGGTCTGTTCTTTtccgcccccc                                                       >  1:148091/1‑140 (MQ=255)
                                                                                                 attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCCGTCAGGGGAATCCACCCCGACCCGATAAAGCCCTTCCGCTTTTGCCGTTGCTGCCCCCCGCGCTCTGTCCTTTGCCGCCCGCAActc                                                   >  2:458926/1‑148 (MQ=255)
                                                                                                    attCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTa                                              <  1:247289/150‑1 (MQ=255)
                                                                                                    attCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTAt                                             <  2:1080252/151‑1 (MQ=255)
                                                                                                       ctcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATc                                            <  2:736405/149‑1 (MQ=255)
                                                                                                         cGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCa                                          <  1:729367/149‑1 (MQ=255)
                                                                                                         cGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGc                                        <  2:1126469/151‑1 (MQ=255)
                                                                                                         cGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCACCAGCAACTCTGCTATCCAGc                                        <  2:1038612/151‑1 (MQ=255)
                                                                                                           ttAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGct                                       <  1:1078893/150‑1 (MQ=255)
                                                                                                                 cGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGc                                <  2:416534/151‑1 (MQ=255)
                                                                                                                   ccAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACATACCCGGCAGGGGAATTCACCCCGGCCAGAAACAACCCCTCCCCTTTTTCCGTTGCTGCCCCCCGCGGTCTTTTCCTTTCCCCCCGCAACTCCGCTATTCCGCTCTcccg                                 >  2:153422/1‑146 (MQ=255)
                                                                                                                    cAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATTATACCGGGCAGGGGAAACCACCCCGACCAGATACAACCCCTCCCCTTTTTCCGGTTCTGCCCCCCGCGGTCTGTCCTTTGCCGCCCGCCACCCTGCTATCCAGCTCCCCGGCTgg                             >  1:1198106/1‑151 (MQ=255)
                                                                                                                      gAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGtt                           <  1:878758/151‑1 (MQ=255)
                                                                                                                      gAAATAGCGGGCCCATTGGCGGCTTTTGGAGAACAAACCGGTCAGGGTAATCCACCCCCGCCACAAACAGCCCTTCCGCTTTTGCCGTTTCTTCCCCCCGCGGTCTGGCCTTTTCCGCCAGCAACTCTTCTTTCCAGCTCTcccggtg                              >  2:414267/1‑143 (MQ=255)
                                                                                                                         aTAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTtgtg                        <  2:198597/151‑1 (MQ=255)
                                                                                                                         aTAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTtgtg                        <  1:422691/151‑1 (MQ=255)
                                                                                                                           aGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTgg                       <  2:1078672/150‑1 (MQ=255)
                                                                                                                             cGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCa                    <  2:432106/151‑1 (MQ=255)
                                                                                                                             cGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCa                    <  1:941451/151‑1 (MQ=255)
                                                                                                                               ggCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCAcc                  <  2:579170/151‑1 (MQ=255)
                                                                                                                                   cATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACt              <  1:234690/151‑1 (MQ=255)
                                                                                                                                        gCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCAt         <  2:241023/151‑1 (MQ=255)
                                                                                                                                           gTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAg      <  2:674459/151‑1 (MQ=255)
                                                                                                                                              tttGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAgg     <  1:1095492/150‑1 (MQ=255)
                                                                                                                                               ttGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGc    >  1:307971/1‑150 (MQ=255)
                                                                                                                                                tGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCt   <  1:338256/150‑1 (MQ=255)
                                                                                                                                                tGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTg  >  2:328116/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                                       
CGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTG  >  NC_000913/2434717‑2435010

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