Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 781,061 C→T 8.7% intergenic (+143/‑7) lysW → / → valZ tRNA‑Lys/tRNA‑Val

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913781,0610CT 8.7% 376.4 / 15.4 114intergenic (+143/‑7)lysW/valZtRNA‑Lys/tRNA‑Val
Reads supporting (aligned to +/- strand):  ref base C (56/48);  new base T (5/5);  total (61/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 121 total aligned reads displayed.

CCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTG  >  NC_000913/780912‑781210
                                                                                                                                                     |                                                                                                                                                     
cccACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTGGTATCCAgcgc                                                                                                                                                       >  2:2421444/1‑150 (MQ=255)
cccACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGGTTGTAGTATCCAgcgc                                                                                                                                                       >  1:721253/1‑150 (MQ=255)
 ccaccaGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCa                                                                                                                                                      <  2:2447914/150‑1 (MQ=255)
        ttttAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCgg                                                                                                                                               >  1:2037319/1‑150 (MQ=255)
        ttttAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCgg                                                                                                                                               <  1:2073460/150‑1 (MQ=255)
        ttttAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCgg                                                                                                                                               <  2:1018109/150‑1 (MQ=255)
        ttttAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCgg                                                                                                                                               <  1:2671588/150‑1 (MQ=255)
        ttttAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCgg                                                                                                                                               >  1:81978/1‑150 (MQ=255)
        ttttAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCgg                                                                                                                                               >  2:263700/1‑150 (MQ=255)
        ttttAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCgg                                                                                                                                               <  2:2587379/150‑1 (MQ=255)
         tttAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCggg                                                                                                                                              <  1:393930/150‑1 (MQ=255)
         tttAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCggg                                                                                                                                              <  1:619541/150‑1 (MQ=255)
            aaCATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGa                                                                                                                                           <  1:2379258/150‑1 (MQ=255)
                tCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAg                                                                                                                                       >  1:2322492/1‑150 (MQ=255)
                tCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAg                                                                                                                                       >  1:2144462/1‑150 (MQ=255)
                tCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAg                                                                                                                                       >  2:106761/1‑150 (MQ=255)
                  gAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCt                                                                                                                                     <  1:2787177/150‑1 (MQ=255)
                  gAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGGTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTCGCt                                                                                                                                     >  1:2034377/1‑150 (MQ=255)
                     gACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCTCGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAg                                                                                                                                  <  2:1068311/150‑1 (MQ=255)
                     gACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAg                                                                                                                                  <  2:1068209/150‑1 (MQ=255)
                          ttGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTTGGGCGAGCGACGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTggg                                                                                                                             <  2:2364414/149‑1 (MQ=255)
                                aGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCa                                                                                                                       <  2:2519870/150‑1 (MQ=255)
                                         gATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTa                                                                                                              <  2:924709/150‑1 (MQ=255)
                                         gATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCGTa                                                                                                              >  2:1036846/1‑150 (MQ=255)
                                             aCGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAag                                                                                                          >  2:2664384/1‑150 (MQ=255)
                                               gTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAagga                                                                                                        <  1:724496/150‑1 (MQ=255)
                                               gTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAagga                                                                                                        <  2:1152216/150‑1 (MQ=255)
                                               gTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAagga                                                                                                        <  2:1459780/150‑1 (MQ=255)
                                                ttGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAaggag                                                                                                       <  2:516929/150‑1 (MQ=255)
                                                  gCGTCATCAACGGCCCGTAGGGAGAGCGAAGCGAGTCATCCTGCACGACCCACCACTTATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAggg                                                                                                     <  2:322697/150‑1 (MQ=255)
                                                  gCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAggg                                                                                                     <  2:2763003/150‑1 (MQ=255)
                                                  gCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAggg                                                                                                     <  2:313292/150‑1 (MQ=255)
                                                   cGTCAGCATCGGCGAGTTGGGAGAGCTACTCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAgggg                                                                                                    <  1:1414865/150‑1 (MQ=25)
                                                    gTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAggggg                                                                                                   >  2:2039011/1‑150 (MQ=255)
                                                    gTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGAGTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAggggg                                                                                                   >  1:1796800/1‑150 (MQ=255)
                                                       aGCAACGGCCCGTAGGGCCAGCGAAGCGAGTCATCCTGCACGACCCACCACTCATGACGTTGGGTTCGGTGGAAGTAGTTTGGAGTATACAGCGCAGTATCGGGTGATGAGCTCAGCTGGGAGAGCACCTCCCTTAACAGGAGGGGGTcg                                                                                                >  2:1987546/1‑150 (MQ=37)
                                                           aCGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGAAGGTCGGAACGGAGGAAGGAGTTTGTAGTATCCAGCTCAGTATCGGGTGATAAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTcggcgg                                                                                            >  2:974198/1‑150 (MQ=25)
                                                              gCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTTGGTTCGGTGTAAGTATATTGAAGTAACCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCACTTACAAGGACGGGGTCGGCGGTTc                                                                                         >  1:1648549/1‑150 (MQ=25)
                                                              gCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTc                                                                                         >  1:2624695/1‑150 (MQ=255)
                                                                 cGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGAt                                                                                      <  1:942281/150‑1 (MQ=255)
                                                                     gggCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCCGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCg                                                                                  <  1:314042/150‑1 (MQ=255)
                                                                         gagcgaAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGtcat                                                                              <  2:2322492/150‑1 (MQ=255)
                                                                         gagcgaAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGtcat                                                                              <  2:2144462/150‑1 (MQ=255)
                                                                            cgaagcgaGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATcac                                                                           >  1:2582175/1‑150 (MQ=255)
                                                                            cgaagcgaGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATcac                                                                           >  1:1854149/1‑150 (MQ=255)
                                                                                    gTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACt                                                                   >  1:1684233/1‑150 (MQ=255)
                                                                                     tCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTc                                                                  >  1:2745606/1‑150 (MQ=255)
                                                                                     tCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGAAGGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGATCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTc                                                                  >  1:1639124/1‑150 (MQ=255)
                                                                                          cTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTc                                                             >  2:2586942/1‑150 (MQ=255)
                                                                                          cTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTc                                                             >  2:2586815/1‑150 (MQ=255)
                                                                                          cTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTc                                                             >  1:1935032/1‑150 (MQ=255)
                                                                                            gCACGAGCCTCCAGTATTGGCGGGGGGTTCTGTGGTTGTAGATTGTCCTCTCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGt                                                           <  1:2768622/150‑1 (MQ=14)
                                                                                               cGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAg                                                        <  2:1103468/150‑1 (MQ=255)
                                                                                                   ccaccaCTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCa                                                    >  1:1197698/1‑150 (MQ=255)
                                                                                                        acTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTgg                                               >  1:2018948/1‑150 (MQ=255)
                                                                                                        acTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTgg                                               >  2:2416230/1‑150 (MQ=255)
                                                                                                        acTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTgg                                               >  2:1116788/1‑150 (MQ=255)
                                                                                                          tcaTGCCGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTa                                             <  1:1190530/148‑1 (MQ=255)
                                                                                                             tGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACATGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTagag                                          <  1:1139758/150‑1 (MQ=255)
                                                                                                             tGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTagag                                          <  1:504277/150‑1 (MQ=255)
                                                                                                               aCGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCa                                        >  1:1703571/1‑150 (MQ=255)
                                                                                                               aCGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCa                                        >  1:1703459/1‑150 (MQ=255)
                                                                                                                  gTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGtt                                     >  1:1636280/1‑150 (MQ=255)
                                                                                                                   tGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTg                                    >  2:881328/1‑150 (MQ=255)
                                                                                                                   tGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTAGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTg                                    >  2:1456287/1‑150 (MQ=255)
                                                                                                                     ggTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGAc                                  >  2:1666559/1‑150 (MQ=255)
                                                                                                                          ggTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTa                             <  1:2574488/150‑1 (MQ=255)
                                                                                                                            tGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGTTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAAt                           <  1:1736046/150‑1 (MQ=255)
                                                                                                                             ggAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATc                          <  2:265452/150‑1 (MQ=255)
                                                                                                                              gAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCa                         >  2:2366983/1‑150 (MQ=255)
                                                                                                                              gAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCa                         >  2:2366970/1‑150 (MQ=255)
                                                                                                                               aaGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCaa                        <  2:1916313/150‑1 (MQ=255)
                                                                                                                                 gtagtTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAAtt                      <  2:2037319/150‑1 (MQ=40)
                                                                                                                                  tttgttGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTg                     <  2:2515432/146‑1 (MQ=14)
                                                                                                                                    gtTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGt                   <  1:1455611/150‑1 (MQ=38)
                                                                                                                                    gtTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGt                   >  2:2611439/1‑150 (MQ=38)
                                                                                                                                    gtTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGt                   <  1:1455605/150‑1 (MQ=38)
                                                                                                                                      ttGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGt                                              <  2:403501/121‑1 (MQ=17)
                                                                                                                                      ttGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGt                                              >  1:403501/1‑121 (MQ=17)
                                                                                                                                      ttGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGgaga                                                                                                                          <  1:1986218/45‑1 (MQ=38)
                                                                                                                                      ttGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCg                 >  1:1404816/1‑150 (MQ=37)
                                                                                                                                         tagtaTCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTTACTTTTAATCAATTGGTCGctg              >  2:321673/1‑148 (MQ=12)
                                                                                                                                         tagtaTCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAg              >  1:874837/1‑150 (MQ=37)
                                                                                                                                          agtaTCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATGAATTGGTCGCAgg             >  1:96755/1‑150 (MQ=37)
                                                                                                                                            taTCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGtt           >  1:2931220/1‑150 (MQ=37)
                                                                                                                                             aTCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTc          <  2:2582175/150‑1 (MQ=37)
                                                                                                                                             aTCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTc          <  2:835882/150‑1 (MQ=37)
                                                                                                                                               ccAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGa        >  1:1631949/1‑150 (MQ=37)
                                                                                                                                               ccAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGa        <  2:1684233/150‑1 (MQ=37)
                                                                                                                                               ccAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGa        <  2:1684219/150‑1 (MQ=37)
                                                                                                                                                cAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGt                                                           <  1:1657395/98‑1 (MQ=18)
                                                                                                                                                cAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGt                                                           >  2:1657395/1‑98 (MQ=18)
                                                                                                                                                cAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAg                                       >  1:1593813/1‑118 (MQ=17)
                                                                                                                                                cAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAg                                       <  2:1593813/118‑1 (MQ=17)
                                                                                                                                                cAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAg                                       <  1:802059/118‑1 (MQ=17)
                                                                                                                                                 aGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAAt      >  2:627062/1‑150 (MQ=37)
                                                                                                                                                  gcgcAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATc     <  1:1022793/150‑1 (MQ=37)
                                                                                                                                                  gcgcAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATc     <  1:106761/150‑1 (MQ=37)
                                                                                                                                                   cgcAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATcc    <  1:2039011/150‑1 (MQ=37)
                                                                                                                                                    gcAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCt   <  1:2664384/150‑1 (MQ=35)
                                                                                                                                                     |                                                                                                                                                     
CCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTG  >  NC_000913/780912‑781210

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