Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 557,144 C→A V200L (GTA→TTA)  folD ← bifunctional methylenetetrahydrofolate dehydrogenase/methenyltetrahydrofolate cyclohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913557,1440CA96.1% 166.8 / NA 51V200L (GTA→TTA) folDbifunctional methylenetetrahydrofolate dehydrogenase/methenyltetrahydrofolate cyclohydrolase
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (25/24);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACGCACGCCTGTAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTT  >  NC_000913/556901‑557390
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
cAACGCACGCCTGTAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTAAGtgatg                                                                                                                                                                                                                                                  >  1:17411/1‑250 (MQ=255)
cAACGCACGCCTGTAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACAGCAACGCTCAATAGATCGGCATTTTCTAAGtgatg                                                                                                                                                                                                                                                  >  1:451467/1‑250 (MQ=255)
             tAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGAttttt                                                                                                                                                                                                                                      >  1:227284/1‑249 (MQ=255)
                 gtgtTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGaa                                                                                                                                                                                                                                <  1:146681/251‑1 (MQ=255)
                 gtgtTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGa                                                                                                                                                                                                                                 >  1:381266/1‑250 (MQ=255)
                  tgtTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGa                                                                                                                                                                                                                                 >  1:84616/1‑249 (MQ=255)
                          aTCAGCGTGTCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTGCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTGGCCAACGGCAACGATCAATAGATCGGCTTTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGt                                                                                                                                                                                                                       <  2:431858/251‑1 (MQ=255)
                            cAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTACTGTCTGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCCTTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGCTCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCACCTATCATTCGATCGGCCTTTTCTACGTGATGACGCAGATTTTTAGTGAATCGGTGattcc                                                                                                                                                                                                                     >  2:65862/1‑247 (MQ=255)
                            cAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTc                                                                                                                                                                                                                      >  1:355253/1‑250 (MQ=255)
                                  ggCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAAGCACAATTGCGCCTGCTTTGAGCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAg                                                                                                                                                                                                               <  1:438745/251‑1 (MQ=255)
                                    cAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAgtg                                                                                                                                                                                                             <  2:265439/251‑1 (MQ=255)
                                     aaCCGTCATCGGGACAACGACGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAACCACGACGAAGCCCGCAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGACCAATAGATCGGCATTTTCTACGACATGACGCCGATTTTTAGTTAAGCGGAGAGTCACTGAAgtgg                                                                                                                                                                                                            >  2:346494/1‑251 (MQ=255)
                                              cGGGCCAACGCCGCCTGGAACAGGCGTAATGTATGATGCGCGTTTAGCCGCTTCTTCAATCACTACGTCGCCCACAACTTTGCCATTTGCCAGACGGGTGAGGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTTGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCTCTGTAGTGGTGCAACCTg                                                                                                                                                                                                   <  2:328877/251‑1 (MQ=255)
                                                 tccAACGCCGCCGGGAACCGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGCCGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGACTAAAGCCTGGCTTGCCAACGGCAACAATCAATAGATCGGAATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCCCTGTAGTGGTGCAACCTGccc                                                                                                                                                                                                >  2:83959/2‑250 (MQ=255)
                                                 gCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCca                                                                                                                                                                                                <  1:177386/251‑1 (MQ=255)
                                                           ccgGGAACAGGCGTAATGGATGAGGCGTGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTc                                                                                                                                                                                      <  2:429609/251‑1 (MQ=255)
                                                               gAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATg                                                                                                                                                                                  <  1:357210/251‑1 (MQ=255)
                                                                              aTGAGGCGCGTTTAGCCGCGTATTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGcc                                                                                                                                                                    >  2:23595/1‑250 (MQ=255)
                                                                                       gTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTc                                                                                                                                                          >  1:485983/1‑251 (MQ=255)
                                                                                           aGCCGCGTCTTCAAACACGACGTCGCCCCCCACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAACGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTACATGCTCATCGGGCGGCCAACGATATTCGATg                                                                                                                                                      >  2:191936/1‑251 (MQ=255)
                                                                                                        aaCACGACGTCGCCCACAACTTTACCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATcaccac                                                                                                                                         <  1:188832/251‑1 (MQ=255)
                                                                                                         acacGACGTCGCCCACAACTTTGCCATTTTCCAGCCGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGCTGCGCCAATCACCACg                                                                                                                                        >  2:130017/1‑251 (MQ=255)
                                                                                                                  cGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAgg                                                                                                                               <  1:428842/251‑1 (MQ=255)
                                                                                                                    cccACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATACAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGcc                                                                                                                             >  1:305260/1‑251 (MQ=255)
                                                                                                                     ccACAACTTTGCCATTTTCCAGAAGGTTGATGCCGACATCAATCACAAGTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATGTTCTAAGGGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAGCGATATTCGATGCGCCAATCGCCACGGTGTTGAGGCCg                                                                                                                            <  2:213983/251‑1 (MQ=255)
                                                                                                                              tGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATc                                                                                                                    >  1:344829/1‑250 (MQ=255)
                                                                                                                               gCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCaa                                                                                                                  <  1:254321/251‑1 (MQ=255)
                                                                                                                                       ccAGACGGTGGATGCGGACATCAATCACTATTGCGATGTCTTTGATCCAGTTACCGGGAATAAAGCCTGGATTGACAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGAAGATTTTTAGTGAAGCGGTGAGTTACTGTAGTGGTGCAACCTGCCAGCAGCAGTTTCATGCTCATAGGGCGGACAACGATATTCGATGCGACAATCACCACGGCGGTGAGGCCGAAGGTATCAAtgttg                                                                                                             <  1:60195/248‑1 (MQ=255)
                                                                                                                                              gTTGATGCCGGTATCAATCACAATTGCGCCTGCATTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGTCAACGATCAATAGATCGGCATTTTCGAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCTAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCaa                                                                                                   <  2:278841/251‑1 (MQ=255)
                                                                                                                                                   tGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCCTGCTCATCGGGCGGCCAACGCTATTCGCTGCGCCAATCACCACGGCGTTGAGGCCGAGGGTATCAATGTTGTAACGCTCAaccagc                                                                                              >  2:15945/1‑251 (MQ=255)
                                                                                                                                                     ccGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGCTGCGCCAATCACCACGGCGTTGAGGCCGAAGGAATCAATGTTGTAACGCTCAAGCAGCGt                                                                                            >  2:309365/1‑251 (MQ=255)
                                                                                                                                                        aCATCAATCACAATAGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCAAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTACCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGACGTTGAGGCCGAAGGTATCAATGTTTTAACGCTCAAGCAGCGTGAc                                                                                         >  2:312500/1‑251 (MQ=255)
                                                                                                                                                                   aaTTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACcgcgcg                                                                              <  1:254810/251‑1 (MQ=255)
                                                                                                                                                                          ccTTCTTTGATCCAGTCACCGGGAATAAAGCCTTGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGGCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCa                                                                       <  1:65862/251‑1 (MQ=255)
                                                                                                                                                                             tcttTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAg                                                                      <  2:160178/249‑1 (MQ=255)
                                                                                                                                                                                 tGATCCAGTCACCGGGAATAAAGCCTGGATTGCCAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGc                                                                >  2:408187/1‑251 (MQ=255)
                                                                                                                                                                                     ccAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGAc                                                            <  1:457366/251‑1 (MQ=255)
                                                                                                                                                                                             cGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCAGTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGgcgcg                                                    <  1:181842/251‑1 (MQ=255)
                                                                                                                                                                                               ggAATAAAGCCTGGCTTGCCCACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGAAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCCACACGGCGTTGAGGCCGAAGGTCTCAATGATGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGACGCAGGGACGCAGACGCGgggcgcg                                                  >  1:2731/1‑251 (MQ=255)
                                                                                                                                                                                                     aaGCCTGGCTTGCCAACGGCAACGATCCATAGCTCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCAGGCGGCCACCGATATTCGATGCGCAAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGGGCAGGGACGCCGACGCGGCGCGCGCTGGc                                             >  2:320755/1‑250 (MQ=255)
                                                                                                                                                                                                      aGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCATCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGcac                                           <  1:176458/251‑1 (MQ=255)
                                                                                                                                                                                                      aGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGcat                                           <  1:176459/251‑2 (MQ=255)
                                                                                                                                                                                                             cTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACAACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAgacgac                                    >  1:168385/1‑251 (MQ=255)
                                                                                                                                                                                                                gCCATCGGCAACGATCAATAGACCGGCCTTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGGGCAACCTGCCAGCAGCAGTTCCATGCGCATCGGGCGGCCAACGATATGCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATTCCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCTCAgacgacaga                                 <  2:65623/251‑4 (MQ=255)
                                                                                                                                                                                                                gCCAACGGCAACGATAAATCGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGAAGCAGTTCCATGCTCATCGGGCGGCCAACGCTATTCGATGCGCCAATCACCACGCCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCCGCGTGACGATAACGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGAcgaccg                                  >  2:495058/1‑250 (MQ=255)
                                                                                                                                                                                                                    aCGGCAACGATCAATAGATCGGCATTTTCTAAGTGATGACGCAGCTTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGtt                             >  2:136016/1‑251 (MQ=255)
                                                                                                                                                                                                                                 aTAGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTCTCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATggca                 >  2:171294/1‑248 (MQ=255)
                                                                                                                                                                                                                                   aGATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGCTATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAAcc              >  1:275804/1‑251 (MQ=255)
                                                                                                                                                                                                                                    gATCGGCATTTTCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTAAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCGGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCg             <  1:23595/251‑1 (MQ=255)
                                                                                                                                                                                                                                               tCTAAGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACTGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGTATGGAAACCGTCCACGTCtt   <  2:412415/250‑1 (MQ=255)
                                                                                                                                                                                                                                               tCTAAGTGATGACGCAGAATTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATAGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGACGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCCCAGACGACCGACGTTGAAAGGATGGAAACCGTCCACGTCttt  >  2:326859/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
CAACGCACGCCTGTAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTT  >  NC_000913/556901‑557390

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