Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 908,169 C→T R38R (AGG→AGA ybjT ← putative NAD(P)‑binding protein YbjT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913908,1690CT97.1% 233.6 / ‑5.3 71R38R (AGG→AGAybjTputative NAD(P)‑binding protein YbjT
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (39/29);  total (40/31)
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

GCTCATGTGGCGGGGCCTGCAACGAACTAAGAAAGATTAATTGCTTAACTGGTACTTCACGTAGCGCATCGCGGACGTTGAGAGCCACCTGGCGCTCCTGAGCGATAAAATCGCCGCCTTCGCCCATGCTGTGCACCAGAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGCCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTT  >  NC_000913/907928‑908411
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   
gCTCATGTGGCGGGGCCTGCAACGAACTAAGAAAGATTAATTGCTTAACTGGTACTTCACGTAGCGCATCGCGGACGTTGAGAGCCACCTGGCGCTCCTGAGCGATAAAATCGCCGCCTTCGCCCATGCTGTGCACCAGAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGAc                                                                                                                                                                                                                                             <  2:154124/250‑1 (MQ=255)
         gCGGGGCCTGCAACGAACTAAGAAAGATTAATTGCGTAACTGGTACTTCACGTAGCGCATCGCGGACGTTGAGAGCCACCTGGCGCTCCTGAGCGATAAAATCGCCGCCTTCGCCCATGCTGTGCACCAGAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTAGCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTgcc                                                                                                                                                                                                                                   <  2:81125/251‑1 (MQ=255)
                cTGCAACGAACTAAGAAAGATTAATTGCTTAACTGGTACTTCACGTAGCGCATCGCGGACGTTGAGAGCCACCTGGCGCTCCTGAGCGATAAAATCGCCGCCTTCGCCCATGCTGTGCACCAGAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTgccgccgccc                                                                                                                                                                                                                            >  2:523355/1‑250 (MQ=255)
                             agaaagaTTAATTGCTTAACTGGTACTTCACGTAGCGCATCGCGGACGTTGAGAGCCACCTGGCGCTCCTGAGCGATAAAATCGCCGCCTTCGCCCATGCTGTGCACCAGAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGccc                                                                                                                                                                                                              >  2:397791/1‑251 (MQ=255)
                              gaaagaTTAATTGCTTAACGGGTACTTCACGTAGCGCATCGCGGACGTTGAGAGCCACCTGGCGCTCCTGAGCGATAAAATCGCCGCCTTCGCCCATGCTGTGCACCAGAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCt                                                                                                                                                                                                             <  2:388656/251‑1 (MQ=255)
                               aaagaTTAATTGCTTAACTGGTACTTCACGTAGCGCATCGCGGACGTTGAGAGCCACCTGGCGCTCCTGAGCGATAAAATCGCCGCCTTCGCCCATGCTGTGCACCAGAAAATACACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCtt                                                                                                                                                                                                            >  1:173778/1‑251 (MQ=255)
                               aaagaTTAATTGCTTAACTGGTACTTCACGTAGCGCATCGCGGACGTTGAGAGCAACCTGGCGCTCCTGAGCGATAAAATCGCCGCCTTCGCCCATGCTGTGCACCAGAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGCTCTGATGCCCtt                                                                                                                                                                                                            >  1:69349/1‑251 (MQ=255)
                                 agaTTAATTGCTTAACTGGTACTTCACGTAGCGCATCGCGGACGGTGAGAGCCACCTGGCGCTCCTGAGCGATAAAATCGCCGCCTTCGCCCATGCTGTGCACCAGAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTgc                                                                                                                                                                                                          <  2:78601/251‑1 (MQ=255)
                                      aaTTGCTTAACTGGTACTTCACGTAGCGCATCGCGGACGTTGAGAGCCACCTGGCGCGCCTGATCGATAAAATCGCCGCCTTCGCCCATGCTGTGCACCAGAAAATAGACCGTATCGATATCCTGCAACAGTGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTgctggc                                                                                                                                                                                                      <  2:433061/250‑1 (MQ=255)
                                            ttAACTGGTACTTCACGTAGCGCATCGCGGACGTTGAGAGCCACCTGGCGCTCCTGAGCGATAAAATCGCCGCCTTCGCCCATGCTGTGCACCAGAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAgtgg                                                                                                                                                                                               >  2:92245/1‑250 (MQ=255)
                                               aCTGGTACTTCACGTAGCGCATCGCGGACGTTGAGAGCCACCTGGCGCTCCTGAGCGATAAAATCGCCGCCTTCGCCCATGCTGTGCACCAGAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGTTGAGAGTgcg                                                                                                                                                                                            >  2:378224/1‑251 (MQ=255)
                                                     aCTTCACGTAGCGCATGGCGGACGGTGAGAGCCACCTTGCGCTCCTGAGCGATAAAATCGCCGCCTGCGCCCATGCTGTGCACCAGAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTATGATATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCTCACCac                                                                                                                                                                                      <  2:299250/251‑2 (MQ=255)
                                                      cTTCACGTAGCGCATCGCGGACGTTGAGAGCCACCTGGCGCTCCTGAGCGATAAAATCGCCGCCTTCGCCCATGCTGTGCACCAGAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTgcgcccc                                                                                                                                                                                        >  2:276123/1‑245 (MQ=255)
                                                          aCGTAGCGCATCGCGGACGATGAGAGCCACCTGGCGCTCCTGAGCGATAAAATCGCCGCCTTCGCCCATGCTGTGCACCAGAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCt                                                                                                                                                                                 >  1:134243/1‑251 (MQ=255)
                                                           cGTAGCGCATCGCGGACGTTGAGAGCCACCTGGCGCTCCTGAGCGATAAAATCGCCGCCTTCGCCCATGCTGTGCACCAGAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTg                                                                                                                                                                                >  2:523794/1‑251 (MQ=255)
                                                                 gcATCGCGGACGTTGAGAGCCACCTGGCGCTCCTGAGCGATAAAATCGCCGCCTTCGCCCATGCTGTGCACCAGAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAAt                                                                                                                                                                          <  1:549905/251‑1 (MQ=255)
                                                                                agagCCACCTGGCGCTCCTGAGCGATAAAATCGCCGCCTTCGCCCATGCTGTGCACCAGAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCGGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACGGGCAcc                                                                                                                                                           <  1:386580/251‑1 (MQ=255)
                                                                                     cACCTGGCGCTCCTGAGCGATAAAATCGCCGCCTTCGCCCATGCTGTGCACCAGAAAATAGACCGTATCGATATCCTGCAACAGAGCCGGTAGGTTATCCGGCCAGCTGAGATCGCCTTTCTGGCAACTGACGTTTGCCAGATGCAGCTTTGAAAGTCTGTCGACATGACGAGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGccacgaga                                                                                                                                                       >  1:313921/1‑243 (MQ=255)
                                                                                       ccTGGCGCCCCTGGGCGATAAAAGCGCCGCCTTCGCCCATGCTGTGCACCAGAAAATAGACCGTATCGAGATCCTGCAACAGGGCCGGTAGGTTATCCGGCCGGCTGAGATCGACGTTATGGCAACTGACGTTTGCCAGTTGCAGTTTTGCAAGGCTGGCGACATGACGTGCCGGCGCCAGGAGCTGATGCCCATGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATTTAGCCACTGGCACCGAGAACt                                                                                                                                                    <  2:47177/251‑1 (MQ=255)
                                                                                            cgcTCCTGAACGATAAAAACGCCGCCTTCGCCCATGCTGTGCACCAGAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGGAGGTTATCCGGCCTGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACGGGCACCGAGAACTAAAAt                                                                                                                                               <  1:419648/251‑1 (MQ=255)
                                                                                               tCCTGAGCGATAAAATCGCCGCCTTCGCCCATGCTGTGCACCAGAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCg                                                                                                                                            >  1:54241/1‑251 (MQ=255)
                                                                                                 cTGAGCGCTAACATCGCCGCCTTCGCCCATGCTGTGCGCCAGAAAACCGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCAAGCTGAGATCGCCTTTATGCCAACTGACGTTTGCCAGTTGCACCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTACCCACTGGCACCGAGAACTATAATGCGtt                                                                                                                                          >  2:83462/1‑251 (MQ=255)
                                                                                                        agaaaaTCGCCGCCTTCGCCCATGCTGTGAACCAGAAAATAGACCTTATCGAGATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGACCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGGCGACATTAAGTACCGCCGCCAGGATCTGATGCCCGTGCTGGCTGAGTTTGTGCACCATATGCTGACCAATGGAGCCACTGGCACCGAGAACCAAAATTCTTGGCGGCAc                                                                                                                                   <  2:360635/249‑1 (MQ=255)
                                                                                                          aaaaTCGCCGCCTTCGCCCATGCTGTGCACCAGAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGt                                                                                                                                 <  1:313071/251‑1 (MQ=255)
                                                                                                              gcgccgcCTTCGCCCATGCTGTGCACCAGAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGtctct                                                                                                                             <  2:470600/250‑1 (MQ=255)
                                                                                                                  cgccTTCGCCCATGCTGTGCACCAGAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTcc                                                                                                                         <  1:460676/251‑1 (MQ=255)
                                                                                                                    ccTTCGCCCATGCTGTGCACCAGAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGGCGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCtt                                                                                                                       <  1:38646/251‑1 (MQ=255)
                                                                                                                        cGCCCATGCTGTGCACCAGAAAATCGACCGTATCGATATCCTGCAACAGGGTCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAGCTAAAATGCGTTGCGGCACGTCTCTCTCCTTAAcg                                                                                                                   <  2:13788/251‑1 (MQ=255)
                                                                                                                        cGCCCATGCTGTGCACCAGAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGAGCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATTCGTTGCGGCACTTCTCTCTTCTTAAcg                                                                                                                   <  2:366373/251‑1 (MQ=255)
                                                                                                                         gCCCATGCTGTGCACCAGAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGCCTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAAcgc                                                                                                                  >  1:465237/1‑251 (MQ=255)
                                                                                                                               gCTGTGCACCAGAAAATAGCCCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAgg                                                                                                            >  1:113192/1‑251 (MQ=255)
                                                                                                                                cTGTGCACCAGAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGCGCACTAACATGCGTTGCGGCACGTCTCTCTCCTTTACGCGCCacg                                                                                                            >  2:126695/1‑248 (MQ=255)
                                                                                                                                     cACCAGAACATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTATGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTACCGCGCCAGGAATGc                                                                                                       >  1:183793/1‑250 (MQ=255)
                                                                                                                                     cACCAGAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGc                                                                                                       >  1:510591/1‑250 (MQ=255)
                                                                                                                                     cACCAGAAAATAGACCGTATCGATATCCGGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGAGCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAAGGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCa                                                                                                      <  2:165395/251‑1 (MQ=255)
                                                                                                                                       ccAGAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACg                                                                                                    <  1:158503/251‑1 (MQ=255)
                                                                                                                                          gAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCCCCAGGAATGCACGccc                                                                                                 >  2:102698/1‑250 (MQ=255)
                                                                                                                                               tAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTggg                                                                                            >  1:79106/1‑251 (MQ=255)
                                                                                                                                                   ccGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGggcggc                                                                                        <  1:359245/251‑1 (MQ=255)
                                                                                                                                                      tatagATATACTGCACCAGGGACGGGAGGTTACCCGGCCAGCTAAAATCGACTTTATGGCAACTGACGTTTGCCAGTTGAAGCTTTGCACGCATGTCGACATGCCGTGCCGCAGCACGGATCTGACGCCCTTGATGGCTGAGTGTGAGCCCCAGAAGCTGCCCACTGTAGCCACTGGCACCGAGAACTAAAATGCGGTGCGGCACGTCTCTCTCCTTCACGCGCCCGGAATCCACGCCAGTGGGGGGCGAc                                                                                     >  2:534242/5‑251 (MQ=255)
                                                                                                                                                       aTCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACt                                                                                    <  1:523794/251‑1 (MQ=255)
                                                                                                                                                         cGATATCCTGCAACAGGGCCGTTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCGTGCTGGCTGAGTGTGCGCACCAGATGCGGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTGGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTTGGCGGCGACTTc                                                                                  <  2:542007/251‑1 (MQ=255)
                                                                                                                                                          gATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTcc                                                                                 >  2:135611/1‑251 (MQ=255)
                                                                                                                                                          gATATCCCGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTcc                                                                                 <  2:510591/251‑1 (MQ=255)
                                                                                                                                                            tatCCTGCAACAGGGCCGGTAGGATATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCGGCTTTGCAAGTCTGTCGACATGACGTGCAGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCACCAGGAATGCACGCCAGTGGGCGGCGCCTTCCGc                                                                               >  1:347606/1‑251 (MQ=255)
                                                                                                                                                               ccTGCAACAGGGCCGGCAGGTTATCCGGCCAGCTGAGATCGACCTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGcc                                                                              >  2:268739/1‑249 (MQ=255)
                                                                                                                                                                cTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAgt                                                                           >  2:259301/1‑251 (MQ=255)
                                                                                                                                                                  gCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGCCATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAgttg                                                                         >  1:335997/1‑251 (MQ=255)
                                                                                                                                                                  gCAACAGAGCCGGAACGTTATCCGGCAAGCTCAGATCGACTTTATGGCAACCGCCGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGACGACAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGACGCTGACCAATGTAGCCACTGGCACCGAGACCTAAAATGCGTTGCGGCCCGTCTCTCTCCTTAACGCGCCCGGAATGCACGCCAGTGGGCGGCGCCTTccccccgctg                                                                         >  2:226631/1‑243 (MQ=255)
                                                                                                                                                                                       ccGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATg                                                    >  1:279323/1‑251 (MQ=255)
                                                                                                                                                                                        cGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACTTCAAGATGc                                                   <  1:157051/251‑1 (MQ=255)
                                                                                                                                                                                         ggCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCg                                                  >  2:306377/1‑251 (MQ=255)
                                                                                                                                                                                           ccAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTc                                                <  1:196787/251‑1 (MQ=255)
                                                                                                                                                                                            cAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCa                                               <  1:476684/251‑1 (MQ=255)
                                                                                                                                                                                            cAGCTGAGATCGACTTTATGGCAACTGACGATTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGATGAGTGTGCGCCCCAGCTGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGAAGTCAAGATGCGTc                                                >  2:344363/1‑250 (MQ=255)
                                                                                                                                                                                                 gagaTCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAg                                           <  1:135611/250‑1 (MQ=255)
                                                                                                                                                                                                 gagaCCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAgg                                          <  2:230680/251‑1 (MQ=255)
                                                                                                                                                                                                    aTCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCgg                                       >  2:208638/1‑251 (MQ=255)
                                                                                                                                                                                                       gACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGAc                                     >  1:95744/1‑250 (MQ=255)
                                                                                                                                                                                                              tGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCCCGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCg                             >  1:295658/1‑251 (MQ=255)
                                                                                                                                                                                                                   aCTGACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCg                        >  1:60382/1‑251 (MQ=255)
                                                                                                                                                                                                                      gcctgttGCCATTTGCAGCTTTGCAGTCCTGTCGGCATGACTTGCCGCCGCCGGAGTCGGATGCCCTTGCTGGCTGAGTGAGCGCACCAAATGCCGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGTCACGGATCTCTCCTTAAGGCGCCAGGAATGCGCGCCAGTGTGCGGCCCCTTCCGCCAGTTGTTCTCGCGAGACGTCAAGAGGCGGCTCCAGGCGGATAATCGGCGAGGAGtt                      <  2:182428/245‑1 (MQ=255)
                                                                                                                                                                                                                      gACGTTTGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGtt                      >  1:472203/1‑250 (MQ=255)
                                                                                                                                                                                                                            tGCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTAAAGAGGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGc               >  1:531035/1‑251 (MQ=255)
                                                                                                                                                                                                                             gCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGc               <  2:54241/250‑1 (MQ=255)
                                                                                                                                                                                                                             gCCAGTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGc               >  1:347191/1‑250 (MQ=255)
                                                                                                                                                                                                                                 gTTGCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGtt          >  1:238243/1‑251 (MQ=255)
                                                                                                                                                                                                                                    gCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTtct       >  1:70829/1‑251 (MQ=255)
                                                                                                                                                                                                                                    gCAGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTtct       >  2:206789/1‑251 (MQ=255)
                                                                                                                                                                                                                                      aGCTTTGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCg     >  1:55718/1‑251 (MQ=255)
                                                                                                                                                                                                                                          ttGCAAGTCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCtt  <  2:154446/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   
GCTCATGTGGCGGGGCCTGCAACGAACTAAGAAAGATTAATTGCTTAACTGGTACTTCACGTAGCGCATCGCGGACGTTGAGAGCCACCTGGCGCTCCTGAGCGATAAAATCGCCGCCTTCGCCCATGCTGTGCACCAGAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGCCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTT  >  NC_000913/907928‑908411

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