Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,715,921 T→A 7.3% intergenic (‑32/+27) pdxY ← / ← tyrS pyridoxamine kinase/tyrosyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,715,9210TA 7.3% 219.6 / 4.8 77intergenic (‑32/+27)pdxY/tyrSpyridoxamine kinase/tyrosyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (35/36);  new base A (3/3);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

GAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGT  >  NC_000913/1715772‑1716070
                                                                                                                                                     |                                                                                                                                                     
gAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGAttt                                                                                                                                                       >  1:1542680/1‑150 (MQ=255)
   cGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTtctc                                                                                                                                                    <  1:53241/150‑1 (MQ=255)
   cGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTtctc                                                                                                                                                    <  1:852796/150‑1 (MQ=255)
   cGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTtctc                                                                                                                                                    <  2:1371127/150‑1 (MQ=255)
   cGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTtctc                                                                                                                                                    >  2:704477/1‑150 (MQ=255)
     gtgtTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAACCAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGAGTtctcgc                                                                                                                                                  >  1:1033107/1‑148 (MQ=255)
     gtgtTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCAc                                                                                                                                                  <  2:1415786/150‑1 (MQ=255)
     gtgtTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGAGATTTCTCAc                                                                                                                                                  <  2:1258690/150‑1 (MQ=255)
          cAGCGGCCAGACGTTCGCGCCCTGGCGGCGCATCGGAATCTCTGCCGCACTGTTACCCGCATAACCATAAACAAAGTGAGGCGGGATAACGAGAATATTTTTCATCATGTACGCTCTATATAAACACGGGGGAGAGATTAATCACTcccc                                                                                                                                             >  1:626605/1‑150 (MQ=255)
          cAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAAAATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGAGAAATCACTcccc                                                                                                                                             <  2:3586/150‑1 (MQ=255)
           aGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAATCTCTGCCGCACTGTTACCCGCATGACCATATACAACGTGAGACTGGATAGCGAGAATATTTTTCATGATGTTCCCTGTATAAAAACCAGGGGGGGTGAGTTCTCACTccccg                                                                                                                                            >  1:29918/1‑149 (MQ=255)
              ggCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGGGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCttt                                                                                                                                         >  2:371700/1‑150 (MQ=255)
              ggCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGAGAAATCACTCCCCCtgg                                                                                                                                         <  2:787137/150‑3 (MQ=255)
               gCCAGACGTTCGCGCCCAGTCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACGATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGTGGAGTGATTTCTCACTCCCCCTTTc                                                                                                                                        <  1:477183/150‑1 (MQ=255)
                ccAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGAGAAATCACTCCCCCtggtt                                                                                                                                       <  1:1390244/150‑5 (MQ=255)
                  aGACGGTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCAc                                                                                                                                     <  2:612451/150‑1 (MQ=255)
                     cGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTa                                                                                                                                  >  1:440898/1‑150 (MQ=255)
                             cccAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCatta                                                                                                                          >  2:758394/1‑150 (MQ=255)
                                    ggcgCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGTGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAg                                                                                                                   <  1:382082/150‑1 (MQ=255)
                                    ggcgCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAg                                                                                                                   >  2:961068/1‑150 (MQ=255)
                                     gcgcATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTACCCCTTTCCACTTAATGCATTATTTCCAGc                                                                                                                  >  2:830907/1‑150 (MQ=255)
                                          tCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATc                                                                                                             >  1:659489/1‑150 (MQ=255)
                                           cGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCa                                                                                                            <  2:678731/150‑1 (MQ=255)
                                             gAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGa                                                                                                          <  1:961068/150‑1 (MQ=255)
                                                   ctGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTaa                                                                                                    >  1:461994/1‑150 (MQ=255)
                                                       cGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCt                                                                                                <  2:450444/150‑1 (MQ=255)
                                                        gCACTGTTACCCGCATGACCATAAACAACGTGAGAATGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCtt                                                                                               >  1:754437/1‑150 (MQ=255)
                                                          aCTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAAATCtttt                                                                                             >  1:381334/1‑150 (MQ=255)
                                                           cTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCttttt                                                                                            >  1:560668/1‑150 (MQ=255)
                                                           cTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCttttt                                                                                            <  2:851397/150‑1 (MQ=255)
                                                             gTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTAc                                                                                          >  2:1593716/1‑150 (MQ=255)
                                                             gTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTAc                                                                                          <  1:879605/150‑1 (MQ=255)
                                                               tACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACcg                                                                                        >  1:274816/1‑150 (MQ=255)
                                                                aCCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTAcctt                                                                                       >  2:167301/1‑148 (MQ=255)
                                                                   cGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGAc                                                                                    >  1:150677/1‑150 (MQ=255)
                                                                      aTGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTATTTTTACAGCGAAGCa                                                                                 >  1:321562/1‑150 (MQ=255)
                                                                           cATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAg                                                                            <  1:577839/150‑1 (MQ=255)
                                                                            aTAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAgg                                                                           <  1:110393/150‑1 (MQ=255)
                                                                             tAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGt                                                                          <  2:1099107/150‑1 (MQ=255)
                                                                              aaaCAACGAGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTa                                                                         <  1:1447828/150‑1 (MQ=255)
                                                                                     gTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGAc                                                                  <  1:845545/150‑1 (MQ=255)
                                                                                     gTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGAc                                                                  <  2:440898/150‑1 (MQ=255)
                                                                                      tGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGAcc                                                                 <  2:1648469/150‑1 (MQ=255)
                                                                                      tGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGAcc                                                                 >  2:592841/1‑150 (MQ=255)
                                                                                      tGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGAcc                                                                 >  2:778700/1‑150 (MQ=255)
                                                                                      tGAGAATGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGAcc                                                                 <  2:754437/150‑1 (MQ=255)
                                                                                        tgaCTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTATCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCaa                                                               <  2:1409336/149‑1 (MQ=255)
                                                                                              gATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGAc                                                         <  2:659489/150‑1 (MQ=255)
                                                                                                  gCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtc                                                     >  2:558640/1‑150 (MQ=255)
                                                                                                  gCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtc                                                     >  1:554423/1‑150 (MQ=255)
                                                                                                   cGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtct                                                    >  2:1335032/1‑150 (MQ=255)
                                                                                                   cGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATAATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtct                                                    >  2:710533/1‑150 (MQ=255)
                                                                                                    gagaATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtctt                                                   >  2:630668/1‑150 (MQ=255)
                                                                                                     agaATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtcttc                                                  >  1:1527317/1‑150 (MQ=255)
                                                                                                         tatTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtcttcttct                                              >  2:760157/1‑150 (MQ=255)
                                                                                                           tttttCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCttt                                            <  1:1278706/150‑1 (MQ=255)
                                                                                                              ttCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTaaa                                         <  1:524425/150‑1 (MQ=255)
                                                                                                               tcatcaTGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAaag                                        <  2:461994/150‑1 (MQ=255)
                                                                                                                catcatGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAaaga                                       >  1:716383/1‑150 (MQ=255)
                                                                                                                 atcatGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTACTTCTTTAaagaa                                      >  2:494905/1‑150 (MQ=255)
                                                                                                                  tcatGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAaagaag                                     <  1:778700/150‑1 (MQ=255)
                                                                                                                      gTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTAtt                                 >  1:376400/1‑150 (MQ=255)
                                                                                                                        tCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCa                               <  1:760157/150‑1 (MQ=255)
                                                                                                                         cccTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAg                              >  2:751647/1‑150 (MQ=255)
                                                                                                                           cTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtc                                                     <  2:418784/125‑1 (MQ=255)
                                                                                                                           cTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtc                                                     >  1:418784/1‑125 (MQ=255)
                                                                                                                             gTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATc                          <  2:716383/150‑1 (MQ=255)
                                                                                                                               ataAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCgg                        <  2:150677/150‑1 (MQ=255)
                                                                                                                                 aaaaaCCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGAc                      >  2:833947/1‑150 (MQ=255)
                                                                                                                                 aaaaaCCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAAATCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGAc                      <  1:925360/150‑1 (MQ=255)
                                                                                                                                   aaaCCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTg                    <  2:1527478/150‑1 (MQ=255)
                                                                                                                                   aaaCCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTg                    <  1:1335032/150‑1 (MQ=255)
                                                                                                                                    aaCCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGt                   >  2:1614085/1‑150 (MQ=255)
                                                                                                                                        aGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGTACTGttttt               >  2:1551649/1‑150 (MQ=255)
                                                                                                                                        aGGGGGAGTGAGAAATCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGttttt               >  2:583113/1‑150 (MQ=255)
                                                                                                                                          gtggAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATATCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCa             <  1:751647/148‑1 (MQ=255)
                                                                                                                                          ggggAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCa             <  1:215451/150‑1 (MQ=255)
                                                                                                                                          ggggAGTGAGAAATCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCa             >  1:717863/1‑150 (MQ=255)
                                                                                                                                           gggAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCAc            >  1:900934/1‑150 (MQ=255)
                                                                                                                                           gggAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCAc            <  2:994836/150‑1 (MQ=255)
                                                                                                                                                tGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTa       <  1:682367/150‑1 (MQ=255)
                                                                                                                                                tGATTTATCACTCCCCCTTTCCACTTAATTCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAATGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTa       <  2:1178768/150‑1 (MQ=255)
                                                                                                                                                  aTTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAAt     >  2:1478752/1‑150 (MQ=255)
                                                                                                                                                   tttCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATg    >  2:1328326/1‑150 (MQ=255)
                                                                                                                                                    ttCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATgg   >  2:825804/1‑150 (MQ=255)
                                                                                                                                                     tctcACTCACCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATACTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGCATCGGACTGTTTTTCACCGTTAATGGt  >  1:1357733/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGT  >  NC_000913/1715772‑1716070

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: