Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,434,8820TG42.6% 94.4 / 19.3 62Y252S (TAC→TCC) truAtRNA pseudouridine(38‑40) synthase
Reads supporting (aligned to +/- strand):  ref base T (10/25);  new base G (26/0);  total (37/25)
Fisher's exact test for biased strand distribution p-value = 3.22e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCAT  >  NC_000913/2434733‑2435030
                                                                                                                                                     |                                                                                                                                                    
cAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTa                                                                                                                                                     <  1:506452/151‑1 (MQ=255)
 aGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTaa                                                                                                                                                    <  2:417711/151‑1 (MQ=255)
    tCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCAGGGGAATcc                                                                                                                                                 >  1:288175/1‑151 (MQ=255)
     ccATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAACCGTCAGGGGAATcc                                                                                                                                                 >  2:560870/1‑150 (MQ=255)
     ccATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCa                                                                                                                                                <  2:637334/151‑1 (MQ=255)
         atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCCGGCAGGGGAAaccacccc                                                                                                                                            >  2:650886/1‑149 (MQ=255)
         atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACcg                                                                                                                                             <  2:323177/150‑1 (MQ=255)
         atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACATTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGGAATTCAcc                                                                                                                                              >  1:360379/1‑149 (MQ=255)
           ttGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATAACCGGCAGGGGAATCCAcccc                                                                                                                                            >  1:412704/1‑147 (MQ=255)
            tGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGa                                                                                                                                          <  2:62334/150‑1 (MQ=255)
              ccTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGGTCATACCGGGCAGGGGAAAACACcgcggc                                                                                                                                         >  1:137186/1‑147 (MQ=255)
              ccTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGATTTTGCAGATCATACCGGGCAAGGTAAACCACCGCGccc                                                                                                                                        >  1:441099/1‑150 (MQ=255)
                tATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGAACATACCGGGCAGGGGAAACCACCGCGGcc                                                                                                                                        >  2:352727/1‑148 (MQ=255)
                       tGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCCGTCAGGGGAATCCACCGCGACCAGATACAGcc                                                                                                                              >  1:182054/1‑151 (MQ=255)
                       tGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATAACGGTCAGGGGAATCCACCGCGACCAGATACAGcc                                                                                                                              >  1:659369/1‑151 (MQ=255)
                        gTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGccc                                                                                                                             <  2:996048/151‑1 (MQ=255)
                          aCTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACTGCCTGAATTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCt                                                                                                                            <  1:149586/150‑1 (MQ=255)
                                 aTCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAATCCACCCCGGCCAGATACAACCCTTCCGCttt                                                                                                                    >  1:37134/1‑151 (MQ=255)
                                       cAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCACTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGt                                                                                                              <  2:814721/151‑1 (MQ=255)
                                          gCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTtgc                                                                                                           <  1:585532/151‑1 (MQ=255)
                                           cTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAACGGTTCGGGGGATACACCCCGACCCGATTAAACCCTTTCCCCTTTTCCggt                                                                                                             >  2:695734/1‑146 (MQ=255)
                                           cTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGAACATACCGGTCAGGGGAATCCACCGCGGCCCGATAAAACCCCTCCCCTTTTTCCggt                                                                                                             >  1:182240/1‑146 (MQ=255)
                                             gTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTtgctgc                                                                                                        <  2:994044/151‑1 (MQ=255)
                                             gTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTtgctg                                                                                                         <  1:1096068/150‑1 (MQ=255)
                                                  tACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAATACACCCCGACCAGATACAGCCCCTCCGCCTTTTCCGTTGCTgccccc                                                                                                    >  2:1021635/1‑150 (MQ=255)
                                                   aCGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCa                                                                                                   <  1:593697/150‑1 (MQ=255)
                                                    ctttGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAACGGGCAGGGGAAACCACCGCGGCCCGATACAACCCTTCCCCTTTTTCCGTTGCTgccgcc                                                                                                    >  1:594963/3‑148 (MQ=255)
                                                    cGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCCGGGGAATCCACAGCGACCAAATACACCCCTTCCGCCTTTTCCGTTGCTgccgc                                                                                                     >  1:630683/1‑147 (MQ=255)
                                                        ggCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGAACATAACGGGCAGGGGAATTCCCCGCCACCAGATACAGACCTTCCCCTTTTTCCGTTGCTGCCGCCAGCggt                                                                                              >  1:890077/1‑148 (MQ=255)
                                                           tACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAACCGGCAGGGGAATCCACCGCGACCAGATACAGCCCCTCCCCTTTTTCCGTTGCTTCCGCCAGCGTTCTGt                                                                                          >  1:967138/1‑151 (MQ=255)
                                                             cAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATAATACCCGGCAGGGGAAACCACCGCGACCAGATACACCCCCTCCGCCTTTTCCCTTGCTGCCCCCAGCGGTCt                                                                                            >  2:529871/1‑147 (MQ=255)
                                                                   gAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGGCAGGGGAAACCACCCCGACCCGATAAAGCCCCTCCCCTTTTTCCGTTTCTGCCCCCCGCGTTCTGTCCTTTtcc                                                                                  >  1:1024794/1‑151 (MQ=255)
                                                                        gTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTTAAGATCATACCCGTCAGGGGAATCCACCCCGCCCCGATACAGCCCCTCCCCCTTTTCCGTTTCTGCCGCCCGCGGTCTTTCCTTTtccccc                                                                               >  2:423910/1‑149 (MQ=255)
                                                                        gTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCa                                                                              <  2:242584/150‑1 (MQ=255)
                                                                          tCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGTAAATCACCCCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCCCCCGCCTTCTGTCCTTTGCCGCCAGc                                                                            >  1:717560/1‑150 (MQ=255)
                                                                           ccGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCa                                                                           <  1:736431/150‑1 (MQ=255)
                                                                                cATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATAATACCGGGCAGGGGAATCCACCGCGACCAGATACAACCCCTCCCCTTTTTCCGGTGCTTCCCCCCGCGTTCTTTCCTTTTCCCCCCGcca                                                                          >  2:959168/1‑144 (MQ=255)
                                                                                 attattCTCGTTAGTCCGCCAGAAATAGCGGTCCCATTGGCGGTTTTTGAAGAACATAACAGTGAGGGCAAAACAGCGGGCACAGAAACAGCCTGGTCGCCTTTCCCGGCCCCCCCGCCCGCAGTCTGTCCTTTTccgcccgctccccttc                                                                    >  1:481321/1‑143 (MQ=255)
                                                                                  ttattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGGCAGGGGAATACCCCCCGACCCGATAAAGCCCCTCCGCCTTTTCCGTTGCTGCCGCCAGCGTTCTTTCCCTTTCCGCCCGCAACTCTg                                                                     >  2:736052/1‑149 (MQ=255)
                                                                                   tattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAACCCTTCCGCTTTTTCCGGTGCTGCCGCCAGCGGTTTTTCCTTTTCCGCCAGCAACTCTg                                                                     >  1:78459/1‑148 (MQ=255)
                                                                                         cGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCCGCAACTCTGCTATCCAGc                                                            >  1:203811/1‑151 (MQ=255)
                                                                                                  gCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCt                                                   <  2:479621/151‑1 (MQ=255)
                                                                                                   ccAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTg                                                  <  2:711395/151‑1 (MQ=255)
                                                                                                   ccAGAAATAGCGGGCCCATTGGCGGCTTTTGGAGATCATACCCGGCAGGGGAAACCACCCCGACCCCATACAGCCCCTCCCCCTTTTCCGGTGCTGCCCCCCGCGCTTTGTCCTTTTCCGCCCGCAACCCTGCTTTCCCGCTCTcccg                                                     >  2:463446/1‑146 (MQ=255)
                                                                                                    cAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTTCCGTTGCTGCCGCCAGCGGTCTTTCCTTTTCCGCCAGCAACTCTGCTATCCAACTCTCCGGCTgg                                                 >  1:31224/1‑151 (MQ=255)
                                                                                                                       ggCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTcc                               <  1:579932/150‑1 (MQ=255)
                                                                                                                           gTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATc                            >  2:1063342/1‑149 (MQ=255)
                                                                                                                             tttGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAgg                         <  1:259078/150‑1 (MQ=255)
                                                                                                                               tGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCt                       <  2:609673/150‑1 (MQ=255)
                                                                                                                                ggAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGc                     >  2:422304/1‑151 (MQ=255)
                                                                                                                                     aTCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACa                <  2:100591/151‑1 (MQ=255)
                                                                                                                                      tCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACa                >  1:745541/1‑150 (MQ=255)
                                                                                                                                        aTACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAat              >  1:598509/1‑150 (MQ=255)
                                                                                                                                          aCCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATAtt           >  2:830265/1‑151 (MQ=255)
                                                                                                                                            cGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTcc         <  2:801168/151‑1 (MQ=255)
                                                                                                                                                cAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGa      <  2:883151/150‑1 (MQ=255)
                                                                                                                                                  gggTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAcc    >  1:1046812/1‑150 (MQ=255)
                                                                                                                                                  gggTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCa   <  2:261543/151‑1 (MQ=255)
                                                                                                                                                   ggTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCat  <  2:412689/151‑1 (MQ=255)
                                                                                                                                                   ggTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTCGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCat  <  2:334262/151‑1 (MQ=255)
                                                                                                                                                    gTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCat  >  2:425889/1‑150 (MQ=255)
                                                                                                                                                    gTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCat  >  2:973641/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCAT  >  NC_000913/2434733‑2435030

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