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

GGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATA  >  NC_000913/2434735‑2435031
                                                                                                                                                   |                                                                                                                                                      
ggTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAAt                                                                                                                                                     <  2:1219013/151‑1 (MQ=255)
 gTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGTCAGGGGaa                                                                                                                                                      >  1:789775/1‑149 (MQ=255)
    cATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATacac                                                                                                                                                 >  1:101283/1‑151 (MQ=255)
     atatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGAAAATAACGGGCAGGGGAATcccac                                                                                                                                                >  1:977985/1‑148 (MQ=255)
       atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACcgc                                                                                                                                              <  2:852227/151‑1 (MQ=255)
       atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACcg                                                                                                                                               <  1:840112/150‑1 (MQ=255)
           gCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGAGCATACCAGGCAGGGGAATTCAcc                                                                                                                                                >  1:27114/1‑145 (MQ=255)
           gCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTTGGAGATCATAACGGTCAGGGGAAACCACCCCCAcc                                                                                                                                          >  1:173656/1‑151 (MQ=255)
              tATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCAAACCGGGCAGGGGAAACCACCGCGAcc                                                                                                                                          >  1:590138/1‑148 (MQ=255)
              tATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGGCCGGGGAAACCACCCCCAcc                                                                                                                                          >  1:175576/1‑148 (MQ=255)
                 gCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAATCCACCCCGAcc                                                                                                                                          >  1:525237/1‑145 (MQ=255)
                      gTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGcc                                                                                                                               <  2:1115071/150‑1 (MQ=255)
                        aCTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCtt                                                                                                                            <  2:2874/151‑1 (MQ=255)
                        aCTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTGGGAGATCATAACGGGAAGGGGAATCCACCCCGACCCGATACAAccc                                                                                                                              >  2:688236/1‑149 (MQ=255)
                         cTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAACGGTCAGGGGAATCCACCCCGACCAGATACAGCCCtt                                                                                                                            >  1:125694/1‑150 (MQ=255)
                         cTCGTCATCTTTCAGGCCGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCCGTCAGGGGAATCCACCGCGACCAGATACAGCCCtt                                                                                                                            >  2:839180/1‑150 (MQ=255)
                           cGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTc                                                                                                                           <  2:934145/149‑1 (MQ=255)
                              cATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCt                                                                                                                       <  2:1299591/150‑1 (MQ=255)
                              cATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCt                                                                                                                       <  2:255277/150‑1 (MQ=255)
                                   ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCttt                                                                                                                     <  2:585378/147‑1 (MQ=255)
                                   ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGc                                                                                                                  <  1:1135655/150‑1 (MQ=255)
                                   ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAACCCTTCCGCTTTTGcc                                                                                                                 >  2:341332/1‑151 (MQ=255)
                                    tCAGGCTGTAACTACGTTGGCTACAGCTTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGcc                                                                                                                 <  2:1081278/150‑1 (MQ=255)
                                    tCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGcc                                                                                                                 <  1:1286450/150‑1 (MQ=255)
                                     cAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGt                                                                                                               <  1:624648/151‑1 (MQ=255)
                                        gCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTtgc                                                                                                            <  2:663409/151‑1 (MQ=255)
                                        gCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTTCCGtt                                                                                                              >  2:951196/1‑149 (MQ=255)
                                         cTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATTATACCGGGCAGGGGAATCCACCCCCCACAGATACAACCCCTTCCCCtttttcggt                                                                                                               >  1:1249689/1‑143 (MQ=255)
                                         cTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGGAAACCACCCCGACCAGATAAAACCCTTCCGCTTTTTCCGtt                                                                                                              >  2:41065/1‑148 (MQ=255)
                                          tGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAAACCACCGCGACCAGATACAACCCTTCCCCCTTTTCCggt                                                                                                              >  1:1058539/1‑145 (MQ=255)
                                          tGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATAATAACGGTCAGGGGAATCCAACGCGACCAGAAACAGCCCTTCCGCTTTTGGCGGtgct                                                                                                           >  2:959630/1‑150 (MQ=255)
                                           gTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTtgct                                                                                                           <  2:893620/149‑1 (MQ=255)
                                                tACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGGAAACCACCGCGACCAGATACAACCCTTCCGCTTTTGCCGTTTCTgccgcc                                                                                                     >  2:1238242/1‑150 (MQ=255)
                                                 aCGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAACCCTTCCGCTTTTGCCGTTGCTgccccc                                                                                                     >  2:429780/1‑149 (MQ=255)
                                                  cGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATTCTACCTGGCCGGGCGATCAAACGCGGCAAGATAGAGCCCCTCCCCTTTTGCCCAtgctgtcccc                                                                                                     >  2:731371/1‑143 (MQ=255)
                                                   gTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCCCTTTTTCCGTTTCTGCCGCCcgcg                                                                                                 >  2:769056/1‑151 (MQ=255)
                                                    ttGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGAACATTCCCGGCCGGGGAAAACAACCCGACCCAAAAAAGCCCTTCCCCCTTTTCCGCTTGTGCCCCCcgcg                                                                                                 >  1:567905/1‑150 (MQ=255)
                                                     tGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCCGTCAGGGTAATCCACCGCGACCCGATACAGCCCCTCCGCTTTTGCCGTTGGTGCCCCCAGCggt                                                                                               >  1:160380/1‑149 (MQ=255)
                                                      ggCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGtt                                                                                               <  1:438413/150‑1 (MQ=255)
                                                       gCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCt                                                                                             <  2:1334144/151‑1 (MQ=255)
                                                       gCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCt                                                                                             <  1:1205211/151‑1 (MQ=255)
                                                             gCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCtt                                                                                       <  1:48388/151‑1 (MQ=255)
                                                                   aaTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGGCAGGGGAATCCACCGCGGCCAGATACAGCCCTTCCGCTTTTTCCGTTGCTGCCGCCCGCGTTCTGTCCTTTgcc                                                                                   >  2:1085362/1‑149 (MQ=255)
                                                                     tGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAACGGTCAGGGTAATCCACCGCGACCAGATAAAGCCCCTCCGCTTTTTCCGTTGCTGCCCCCCGCGTTCTGTCCTTTtccgcc                                                                                >  2:225741/1‑150 (MQ=255)
                                                                           gAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGGCCGGGGAATCCACCGCGGCCAGATACAGCCCCTCCGCCTTTGCCGTTGCTGCCCCCCGCGGTCTGTCCCTTGCCGCCCGCAAc                                                                          >  1:977931/1‑150 (MQ=255)
                                                                           gAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGAACATACCGGGCAGGGGAATCCACCGCGACCAGATAAAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCCGCGGTCTGTCCTTTTCCGCCAGCaa                                                                           >  2:1150757/1‑149 (MQ=255)
                                                                           gAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTCGAAAAACATTTCCGGCAGGGGAAACCACCCCGAACCAATAAAGCCCCTCCCCCTTTTCCGTTTCTGCCGCCAGCGTTCTGGCCTTTgccccc                                                                                >  2:974381/1‑144 (MQ=255)
                                                                              cATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGCAGATCATACCGGGCAGGGGAATCCACCGCGCCCAGATACAGCCCCTCCCCTTTTGCCGTTGCTGCCGCCAGCGCTCTGTCCTTTtccgcc                                                                                >  1:1278616/1‑141 (MQ=255)
                                                                                   ttCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCAGGGGAATCCACCCCGACCCGATACAGCCCTTCCGCCTTTTCCGTTGCTGCCGCCCGCGGTCTGTCCTTTTCCGCCAGCAACTCTTCTAt                                                                  >  2:1264270/1‑150 (MQ=255)
                                                                                    tctcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCGGGGGAATCCACCGCCACCCGAAACAGCCCCTCCCCCTTTTCCGGTGCTTCCGCCCGCGTTCTGTTCTTTGCCGCCCGCCActct                                                                       >  1:836990/1‑144 (MQ=255)
                                                                                     ctcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCCCTTTTTCCGTTGCTGCCGCCCGCGTTCTGTCCTTTTCCGCCCGCAACTCTGCTATcc                                                                >  2:711960/1‑150 (MQ=255)
                                                                                        gTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGGCAGGGGAATACCCCCCGCCCAGATACAGCCCCTCCCCTTTTTCCGGTGCTGCCGCCAGCGGTCTGTTCTTTGCCCCCAGCCACTCTGCTATCCAGc                                                             >  1:699625/1‑150 (MQ=255)
                                                                                          tAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGGAAACCACCGCGACCAAATACACCCCTTCCGCTTTTGCCCTTGCTGCCGCCAGCGTTCTGTCCTTTTCCGCCAGCAACATTGCCATCCCGctct                                                          >  2:1139617/1‑151 (MQ=255)
                                                                                          tAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGctc                                                           <  1:1267642/150‑1 (MQ=255)
                                                                                           aGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGAACATACCGGGTAGGGGAATCCACCGCGACCAGATACAGCCCTTCCCCCTTTGCCGTTTCTGTCCCCAGCGGTCTGTCCTTTGCCGCCCGCAACtcttctcttta                                                               >  1:528528/1‑140 (MQ=255)
                                                                                            gTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGAACATACCGGGCCGGGGAAACCAACGCGAACAGATACAACCCTTCCGCTTTTTCCGTTGCTGCCGCCAGCGTTTTGTCCTTTGCCGCCAGCAACACTGCCATCCCGCCCTcc                                                        >  2:1233100/1‑151 (MQ=255)
                                                                                                     aaaTAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTTCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTTCTATCCAGCTCTCCGGGTgggt                                                >  2:287637/1‑148 (MQ=255)
                                                                                                     aaaGAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGtt                                                <  1:1089288/150‑1 (MQ=255)
                                                                                                      aaTAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATAAAGCCCTTCCGCTTTTTCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGtt                                                >  2:273299/1‑149 (MQ=255)
                                                                                                          gCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTggg                                           <  1:1047553/150‑1 (MQ=255)
                                                                                                            gggCCCATTGGCGGTTTTGGAAGATCATACCGGGCAGGGTAATCCACCCCGACCCGATACAGCCCCTCCCCTTTTGCCGGTGCTGCCGCCCGCGCTCTGTCCTTTTCCGCCCGCAACCCTGCTCTTCAGCTCTCCGGCTGGTTTTggg                                           >  1:784407/1‑148 (MQ=255)
                                                                                                              gCCCATTGGCGGTTTTGGAAGGTCATAACCGGCAGGGGAATCCACCGCGACCCGATACAGCCCTTCCCCTTTTGCCGTTTCTGCCCCCAGCGCTCTGTCCCTTTCCGCCCGCAACTCTTCTATTCAGCTCTCCGGCTGGGTTTGGGGAc                                        >  1:643667/1‑149 (MQ=255)
                                                                                                                     ggCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTcc                                >  2:344085/1‑150 (MQ=255)
                                                                                                                        ggTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCa                            <  2:102084/151‑1 (MQ=255)
                                                                                                                         gTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATc                             >  1:1200246/1‑149 (MQ=255)
                                                                                                                         gTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATc                             >  2:219536/1‑149 (MQ=255)
                                                                                                                             tGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCt                        <  2:926458/150‑1 (MQ=255)
                                                                                                                              ggAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTg                       <  2:318718/150‑1 (MQ=255)
                                                                                                                              ggAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGc                      >  2:102232/1‑151 (MQ=255)
                                                                                                                               gAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGcc                     >  1:297630/1‑151 (MQ=255)
                                                                                                                                aaGATCATACCGGTCAGGGTAATCCACCGCGACTAGATACAGCCCTTCCGCTTTGGCTGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACCTCCATCAGGCTGc                      >  1:946263/1‑149 (MQ=255)
                                                                                                                                 aGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCg                    >  2:530087/1‑150 (MQ=255)
                                                                                                                                         ccGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATAtt            <  1:712920/150‑1 (MQ=255)
                                                                                                                                         ccGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTc           >  2:124335/1‑151 (MQ=255)
                                                                                                                                           ggTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCt         <  1:50529/151‑1 (MQ=255)
                                                                                                                                           ggTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCt         >  1:646010/1‑151 (MQ=255)
                                                                                                                                            gTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTg        <  1:55648/151‑1 (MQ=255)
                                                                                                                                            gTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTg        <  2:558669/151‑1 (MQ=255)
                                                                                                                                            gTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCATg        <  2:593720/151‑1 (MQ=255)
                                                                                                                                                  gTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCata  <  1:287864/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATA  >  NC_000913/2434735‑2435031

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: