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

AAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCT  >  NC_000913/2434712‑2435009
                                                                                                                                                     |                                                                                                                                                     
aaaaTCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTg                                                                                                                                                       <  1:759940/150‑1 (MQ=255)
  aaTCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGaa                                                                                                                                                    <  1:1170161/151‑1 (MQ=255)
        tGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCa                                                                                                                                               >  1:236597/1‑150 (MQ=255)
           ggAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAcc                                                                                                                                           <  1:749618/151‑1 (MQ=255)
            gAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATAcc                                                                                                                                           >  1:1110546/1‑150 (MQ=255)
            gAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCg                                                                                                                                          <  1:1041731/151‑1 (MQ=255)
             aaataaatCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACATAcccg                                                                                                                                         >  2:1301557/1‑149 (MQ=255)
              aataaatCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCAt                                                                                                                                              >  2:885400/1‑145 (MQ=255)
               ataaatCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCaaaaccg                                                                                                                                         >  1:300592/1‑143 (MQ=255)
                taaatCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATAcccg                                                                                                                                         >  1:1123268/1‑146 (MQ=255)
                   atCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAggg                                                                                                                                   <  2:336924/151‑1 (MQ=255)
                    tCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGTCAgggg                                                                                                                                  >  2:343622/1‑150 (MQ=255)
                      aGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATAACGGTCAggg                                                                                                                                   >  1:2377/1‑148 (MQ=255)
                      aGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGGCAggg                                                                                                                                   >  2:1455293/1‑148 (MQ=255)
                          ccATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAATcc                                                                                                                             >  1:489026/1‑150 (MQ=255)
                          ccATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACATACCGGGCAGGGGAAttc                                                                                                                             >  1:17146/1‑148 (MQ=255)
                          ccATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGCTCATAACGGTCAGGGGaa                                                                                                                                >  2:1029815/1‑147 (MQ=255)
                          ccATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATcc                                                                                                                             <  1:138402/150‑1 (MQ=255)
                          ccATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGACCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCCGGGGAATcc                                                                                                                             >  1:456975/1‑150 (MQ=255)
                           cATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCAc                                                                                                                           <  1:786975/151‑1 (MQ=255)
                            atatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAAATCATACCCGGCCGGGGAATCCAcc                                                                                                                          >  2:864628/1‑151 (MQ=255)
                             tatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAATCCAcc                                                                                                                          >  1:254646/1‑150 (MQ=255)
                             tatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATCCAcc                                                                                                                          >  1:276235/1‑150 (MQ=255)
                              atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACcgc                                                                                                                        <  1:1343643/151‑1 (MQ=255)
                                  gCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGAcc                                                                                                                    <  2:814517/151‑1 (MQ=255)
                                   ccTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGGCAGGGGAATCCACcccg                                                                                                                       >  1:1098071/1‑147 (MQ=255)
                                    cTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACATAACCGTCAGGGGAATCCACCCCGAcc                                                                                                                    >  2:1391652/1‑149 (MQ=255)
                                      aTGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAACGGTCAGGGGGAACCACCCCGGcc                                                                                                                    >  2:1458717/1‑147 (MQ=255)
                                      aTGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCCGGCCGGGGAATCCACCGCGAcc                                                                                                                    >  1:1286156/1‑147 (MQ=255)
                                         ccTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCAAAACCGGCAGGGGAATCCACCGCGACCAGATa                                                                                                               >  1:802062/1‑149 (MQ=255)
                                           ttGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACa                                                                                                             >  1:374521/1‑149 (MQ=255)
                                            tGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAg                                                                                                           <  2:582257/149‑1 (MQ=255)
                                             gTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGccc                                                                                                        <  2:37236/151‑1 (MQ=255)
                                                cTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCCGGCCGGGGAATCCACCCCGACCCGATACAGcccct                                                                                                      >  2:380186/1‑148 (MQ=255)
                                                  cGTCATCTCTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTc                                                                                                     <  1:1074809/149‑1 (MQ=255)
                                                     cATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCt                                                                                                 <  2:1089795/150‑1 (MQ=255)
                                                        cTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCAGGGGAATCCACCCCGGCCAGATACAACCCTTCCGCtttt                                                                                              >  1:467865/1‑150 (MQ=255)
                                                          ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATAATACCCGGCAGGGGAATCCACCCCGACCACATACAGCCCCTCCCCcttt                                                                                              >  2:986997/1‑148 (MQ=255)
                                                          ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCtttt                                                                                              <  1:141297/148‑1 (MQ=255)
                                                          ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTg                                                                                             <  2:1384610/149‑1 (MQ=255)
                                                                  gTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTTCAGATGATACCGCCGGGGGGAATTCACCCCGAACAGAAAAATGCCCTTCTCCtttgtc                                                                                            >  2:1255681/1‑140 (MQ=255)
                                                                    aaCTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCCGGGGAAACCCCCCCGGCCCGATACAGCCCTTCCGCTTTTGCCGTTGCTgcc                                                                                  >  2:401649/1‑150 (MQ=255)
                                                                    aaCTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTgcc                                                                                  <  1:1152390/150‑1 (MQ=255)
                                                                      cTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAACCGGCCGGGGAATCCACCCCGACCAGATACAACCCTTCCGCTTTTTCCGTTGCTgccgcc                                                                               >  2:1021058/1‑151 (MQ=255)
                                                                        aCGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCCGTCAGGGGAAACCCCCGCGACCAGATACAGCCCTTCCCCTTTTTCCGTTGCTgccccc                                                                               >  1:159714/1‑149 (MQ=255)
                                                                         cGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGc                                                                            <  2:1000759/151‑1 (MQ=255)
                                                                          gTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATTATAACCGTCAGGGGAATCCACCCCGGCCAGATACAGCCCTTTCGCCTTTTCCGTTGCTGCCGCCa                                                                              >  2:977134/1‑148 (MQ=255)
                                                                             ggCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGtt                                                                         <  2:115511/150‑1 (MQ=255)
                                                                              gCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCCTCCGCTTTTGCCGCTGCTGCCGCCcgcggt                                                                         >  1:1229929/1‑147 (MQ=255)
                                                                               cTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATTATACCGGGCAGGGGAAACCACCGCGGCCAGATACAACCCCTCCCCTTTTTCCGGTGCTGGCGCCAGCggt                                                                         >  2:799798/1‑146 (MQ=255)
                                                                                tACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGTAATCCACCCCGACCAGATACAGCCCCTCCGCTTTTGTCGTTGCTGCCGCCCGCGTTCt                                                                       >  1:1202576/1‑149 (MQ=255)
                                                                                tACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATAACGGGCAGGGGAAACCACCCCGACCAGAAAAAACCCTTCCGCCTTTTCCCTTGCTGCCGCCcgcggt                                                                         >  1:599014/1‑145 (MQ=255)
                                                                                   aGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACATACCGGGCAGGGGAAATCACCGCGACCAGATAAAACCCTTCCCCCTTTTTCGTTGCTGCCCCCAGCGGTCTGTcc                                                                   >  1:1249284/1‑150 (MQ=255)
                                                                                     ccTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAGCCCCTCCGCTTTTTCCGTTGCTGCCGCCAGCGGTCTGTCCttt                                                                >  1:466637/1‑151 (MQ=255)
                                                                                         aaaTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAAATCATAACGGTCAGGGGAAACCACCGCGACCAGATACAGACCTTCCGCTTTTGCCGTTGCTGCCCCCAGCGGTCTTTCCTTTgcc                                                             >  1:1107546/1‑150 (MQ=255)
                                                                                            tGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCAGTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTGCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTgcctcc                                                          <  2:223871/150‑1 (MQ=255)
                                                                                             gTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAg                                                        <  1:275354/151‑1 (MQ=255)
                                                                                                 cGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATTATATCGGTCAGGGGAATCCACCGCGACCAGATAAAGCCCCTCCGCTTTTGTCGTTGCTGCCGCCAGCGGTCTGTCCTTTTccgcccgc                                                       >  2:1107856/1‑148 (MQ=255)
                                                                                                    gCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAActct                                                 <  2:1081744/151‑1 (MQ=255)
                                                                                                          ttCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATCCACCCCGACCCGATACAGCCCTTCCGCTTTTTCCGGTGCTGCCCCCCGCGTTCTGTCCTTTGCCGCCAGCCACCCTGCTAt                                            >  2:376017/1‑150 (MQ=255)
                                                                                                            ctcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCCTACCCGTCAGGGGAAACCACCGCGACCAGATACAGCCCCTCCGCTTTTGCCGTTGCTGCCGCCCGCGGTCTGTCCTTTTCCGCCCGCAAATCTGCTATcc                                          >  1:1118835/1‑150 (MQ=255)
                                                                                                            ctcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGGAATCCACCCCGACCAGATACAGCCCTTCCCCTTTTTCCGTTGCTGCCGCCCGCGGTTTTTCCTTTTCCCCCCGCAAATCTGCTATcc                                          >  1:236053/1‑150 (MQ=255)
                                                                                                              cGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATcc                                          <  1:1021083/148‑1 (MQ=255)
                                                                                                              cGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAg                                        <  2:254643/150‑1 (MQ=255)
                                                                                                               gTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGct                                      <  2:1123238/151‑1 (MQ=255)
                                                                                                               gTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATAAAACCCCTCCGCTTTTTCCGTTGCTGCCCCCAGCGGTCTGTCCTTTTTCGCCAGCAACTCTGCTAttc                                          >  1:949746/1‑145 (MQ=255)
                                                                                                                  aGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAATATACCGGTCAGGGGAATTCACCGCGAACAGAAACAACCCTTCCGCCTTTGTCGGTTCTGCCCCCCGCCTTTTGTCCTTTGCCCCCAGCAAATCTGCCTTCCACctc                                     >  1:977960/1‑149 (MQ=255)
                                                                                                                  aGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCtctc                                   <  2:1406142/151‑1 (MQ=255)
                                                                                                                   gTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTcc                                  <  1:254542/151‑1 (MQ=255)
                                                                                                                     ccGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGAACATACCCGGCCGGGGAATCCACCCCCGCCCGATACAGCCCCTCCCCTTTTTCCGTTGCCGCCCCCCCCGGTTTTTCCTTTTTCGCCCGCAACTCTTCTATCCAGctcctc                                  >  2:32867/1‑146 (MQ=255)
                                                                                                                         cAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTgg                            <  1:21240/151‑1 (MQ=255)
                                                                                                                          aGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGt                           <  2:1249252/151‑1 (MQ=255)
                                                                                                                           gAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCAGGGGAATACACCCCGACCAGATACAACCCTTCCGCTTTTTCCGTTTCTGCCCCCCGCGTTCTGTCCTTTTCCGCCCGCCACCCTGCTATTCAGCCCTCCGGCTgg                            >  1:305672/1‑149 (MQ=255)
                                                                                                                           gAAATAGCGGGCCCATTGGCGGTTTTTGAAGATAATAACGGGCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTTCCGTTGCTGCCGCCAGCGGTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTgg                            >  1:1036250/1‑149 (MQ=255)
                                                                                                                           gAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGtt                          <  1:278210/151‑1 (MQ=255)
                                                                                                                              aTAGCGGGCCCATTGGCGGTTTTTGGAGATCAAACCCGTCAGGGGAATCCACCCCGACCAGATACAACCCTTCCGCTTTTTCCGGTGCTGCCGCCCGCGGTCTTTCCTTTTCCGCCCGCCACCCTGCTATCCACCTCTCCCGCTGGgtgt                        >  2:261793/1‑150 (MQ=255)
                                                                                                                               tAGCGGGCCCATTGGCGGTTTTTGGAAATCATACCGGGCAGGGGAATCCACCCCGACCCGATACAACCCTTCCGCCTTTTCCGGTGCTTCCCCCCGCGGTCTTTTCTTTTGCCCCAGCCACCCTGCTTTTCACCTCTCCGGCTGGgtgt                        >  2:1151434/1‑149 (MQ=255)
                                                                                                                               tAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCCGTCAGGGGAATTCACCGCGACCAGATACAGCCCCTCCGCCTTTTCCGGTGCTGCCGCCCGCGGTCTTTCCTTTGCCGCCCGCAACTCTGCTATCCAGCCCTCCCGCTGGGTGTgg                      >  2:325749/1‑151 (MQ=255)
                                                                                                                                aGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTgg                      <  1:864650/150‑1 (MQ=255)
                                                                                                                                 gCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTgg                      <  1:54006/149‑1 (MQ=255)
                                                                                                                                 gCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGc                    <  1:343627/151‑1 (MQ=255)
                                                                                                                                         aTTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACtt            <  1:32867/151‑1 (MQ=255)
                                                                                                                                           tGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTcc          >  1:949088/1‑151 (MQ=255)
                                                                                                                                            ggCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTcc          <  2:489019/150‑1 (MQ=255)
                                                                                                                                            ggCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTcc          <  1:1458761/150‑1 (MQ=255)
                                                                                                                                            ggCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTcc          >  2:897227/1‑150 (MQ=255)
                                                                                                                                                gTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATc       >  1:818189/1‑149 (MQ=255)
                                                                                                                                                   ttGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGc   <  2:1172079/150‑1 (MQ=255)
                                                                                                                                                   ttGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCt  >  2:272880/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCT  >  NC_000913/2434712‑2435009

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