Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,434,8820TG42.9% 157.4 / 25.7 91Y252S (TAC→TCC) truAtRNA pseudouridine(38‑40) synthase
Reads supporting (aligned to +/- strand):  ref base T (9/43);  new base G (39/0);  total (48/43)
Fisher's exact test for biased strand distribution p-value = 2.10e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAC  >  NC_000913/2434738‑2435027
                                                                                                                                                |                                                                                                                                                  
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                            <  1:1458761/150‑1 (MQ=255)
                                                                                                                  ggCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTcc                            <  2:489019/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)
                                                                                                                            gAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGcc                 <  2:1198577/151‑1 (MQ=255)
                                                                                                                            gAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGcc                 <  1:1198927/151‑1 (MQ=255)
                                                                                                                             aaGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCg                >  1:999520/1‑151 (MQ=255)
                                                                                                                               gATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGAc              <  1:165586/151‑1 (MQ=255)
                                                                                                                                aTCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACa             <  2:514788/151‑1 (MQ=255)
                                                                                                                                  cATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACaa            <  2:742929/150‑1 (MQ=255)
                                                                                                                                  cATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAat           <  2:894384/151‑1 (MQ=255)
                                                                                                                                    tACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAatat         >  2:714869/1‑151 (MQ=255)
                                                                                                                                     aCCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATAtt        <  2:467858/151‑1 (MQ=255)
                                                                                                                                      ccGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATAtt        <  2:456968/150‑1 (MQ=255)
                                                                                                                                      ccGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTc       >  2:449367/1‑151 (MQ=255)
                                                                                                                                       cGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTcc      <  2:670777/151‑1 (MQ=255)
                                                                                                                                       cGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTcc      >  2:1302434/1‑151 (MQ=255)
                                                                                                                                         gTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTg    <  2:276232/151‑1 (MQ=255)
                                                                                                                                           cAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAc  <  2:576565/151‑1 (MQ=255)
                                                                                                                                            aGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAc  <  1:1301592/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAC  >  NC_000913/2434738‑2435027

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: