Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 557,635 A→T L36Q (CTG→CAG)  folD ← bifunctional methylenetetrahydrofolate dehydrogenase/methenyltetrahydrofolate cyclohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913557,6350AT97.9% 159.4 / NA 47L36Q (CTG→CAG) folDbifunctional methylenetetrahydrofolate dehydrogenase/methenyltetrahydrofolate cyclohydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (20/26);  minor base G (0/1);  total (20/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCAGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCT  >  NC_000913/557391‑557876
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
gccggaaTGAATACGTCCCAGCTCTTTGACTTTACAAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTg                                                                                                                                                                                                                                             <  1:358737/245‑1 (MQ=255)
        gAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGc                                                                                                                                                                                                                                     >  1:376881/1‑251 (MQ=255)
          aTACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGACTGACATCGATGGTGTTGTAGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGc                                                                                                                                                                                                                                     >  1:178965/1‑249 (MQ=255)
              gTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCg                                                                                                                                                                                                                               >  1:64147/1‑251 (MQ=255)
                                    aacatccgccgGTAAGGGCAGGTGAACCAGAATGCCATCGAAGGTGTTGCGGGCGTTCAGCGTTACGATAATCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCTGGAGGTCTTAATAGCGGGCGACTAACCCGAATTCTTCACAAGCCTAGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCGCAACGGTCGGTCCTGGTGCCCGCAGGCCGGGTGCAATACGCGCCTTAACtt                                                                                                                                                                                                         <  2:398989/246‑1 (MQ=255)
                                           gccgGTAACGGCAGTTGAACCAGAATGCCATCGATGGTCTTGTCAGCATTCAGCGTATTGATAAGCTCCAGCAGCTCCGCTTCGCTGGAGGTTTCCGGGAGGTCATAAGATCGGGAGAAGTACCCTACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCAATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGTCCGAACTTTTTGAGc                                                                                                                                                                                                  >  2:317621/1‑251 (MQ=255)
                                            ccgGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGc                                                                                                                                                                                                  >  1:284136/1‑250 (MQ=255)
                                               gTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATGTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACt                                                                                                                                                                                              <  2:392745/251‑1 (MQ=255)
                                                    ggCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCaga                                                                                                                                                                                         <  1:385162/251‑1 (MQ=255)
                                                            aaCCAGAAGGCCATCGATGGTGTTGACGGTATTGATAGTAGAGATAATCTCCAGAAGGTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCTGGCGACGAACCCGACTTCTGCACAAGCCTGGCGCTTGCTTTCGACATAAATTTGCGATGCATGGTTACTACCCACCAGTTCAACGGCCTGTCCTGCTACCCGCAGTCCGTCTGCAATACCCGCCTGATCTTTTTGAGCAACTTCAGAGCGCCCct                                                                                                                                                                                 <  2:390164/251‑1 (MQ=255)
                                                               cAGATTGCCATCGATGGTGTTGTCGGCATTGAGGGTATCAATAAGCTCCAGGAGCTCCGCTTCGCTGGTGGTTGCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCGTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACACACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCATAGCGCACctgct                                                                                                                                                                              <  2:206702/251‑1 (MQ=255)
                                                                      ccaatgATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTGCGCTGGTGTTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGATTTCTTCACAAGCCTTGCGTTTGCTTGCGAGATAAATTTGCAATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTTCTGCGCAAt                                                                                                                                                                       <  2:325907/246‑1 (MQ=255)
                                                                            aTGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCAGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCCGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGtttt                                                                                                                                                                 >  1:260832/1‑251 (MQ=255)
                                                                                    gtCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGtttt                                                                                                                                                                 >  2:86023/1‑243 (MQ=255)
                                                                                    gtCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGtttt                                                                                                                                                                 <  1:86023/243‑1 (MQ=255)
                                                                                             cAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCttt                                                                                                                                                 <  1:143145/250‑1 (MQ=255)
                                                                                              aGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTGTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCATGTCCTGGTGCCCGCAGTACGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTtgc                                                                                                                                               >  2:259228/1‑251 (MQ=255)
                                                                                                 gTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGGGGTTGACGTGAGGTCATAAGAGCGTGAGACGAACCCTACTTCGTCACAAGCCTTGCCTTTGCTGGCGACATAAATGTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCTCGCAGTCCGGCTGCAATACGCGTCTGAACTGTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAGTCGTTTTTCCTTCAAGAATCTTtgctgc                                                                                                                                            <  2:77449/251‑1 (MQ=255)
                                                                                                            ccccAGCAGCTCCGCTTCGATGGTGGTTTCCGGAAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTATCTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGACCGCACCTGCTGCACAATCGTTTTACCTTCAATAATCTTTGCTGCCATCagagaga                                                                                                                                 <  1:310868/249‑1 (MQ=255)
                                                                                                                    gcTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTcc                                                                                                                          >  1:81273/1‑250 (MQ=255)
                                                                                                                             cGCTGGTGGTTTCCGGGAGGTCATAAGAGGGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACg                                                                                                                <  2:284136/251‑1 (MQ=255)
                                                                                                                                  gtggtTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTCCTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCCCCTGCTGCGCAATCGTTTTCCCGTCAATAATCTTTGCTGCCATCAGATAGAGGATTCCATCTGTTACGTAGAt                                                                                                           <  2:61383/251‑1 (MQ=255)
                                                                                                                                    ggtTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCg                                                                                                         <  1:145497/251‑1 (MQ=255)
                                                                                                                                    ggtTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCAGCTGTTACGTAGATc                                                                                                          <  1:86815/250‑1 (MQ=255)
                                                                                                                                       ttCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCCCAAGACTTGCGTTTGCTTGCGACAAAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGAAATCGTTTGACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGa                                                                                                        >  1:102045/1‑249 (MQ=255)
                                                                                                                                         ccGGGAGGTCAGAAGAGCGGGTGAAGAACCCGACTTCTTCCCAACCCGTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCATCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGTGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAgg                                                                                                     <  2:324619/250‑1 (MQ=255)
                                                                                                                                             gAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATg                                                                                                >  2:37575/1‑251 (MQ=255)
                                                                                                                                                    aaaGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTAt                                                                                         <  2:377374/250‑1 (MQ=255)
                                                                                                                                                     aaGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTAtt                                                                                        >  1:352174/1‑251 (MQ=255)
                                                                                                                                                          cGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTc                                                                                   <  1:234256/251‑1 (MQ=255)
                                                                                                                                                                  gAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCgg                                                                           <  1:236416/251‑1 (MQ=255)
                                                                                                                                                                            tcttcACACGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGt                                                                 >  2:237569/1‑251 (MQ=255)
                                                                                                                                                                             cttcACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTc                                                                >  2:204762/1‑251 (MQ=255)
                                                                                                                                                                               tcACAAGGCTTGCGTTTGCTTGCTACATAAATTTGCGATGCAGGGTTACTACCCCCCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGt                                                                 <  1:369989/248‑1 (MQ=255)
                                                                                                                                                                                 acaAGCCTTGCGTTTGCTTGCGACATAACTTTGCGCTGCACGGTTACTACCCACCACCGCAACGGCCTGTCCTGGTGCCCGCACTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGATCGCACCTGCTGCGCACTCATTTTACCGTCAATAATCTTTGCTGCACTCAGAGAGAGGATTCCATCTGTTACGGAGATCGATGGGGATGCGCCTATTTTTTCCGAAACTGGGCGCGCTGTCAGgt                                                            >  2:182704/1‑251 (MQ=255)
                                                                                                                                                                                   aaGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGgttt                                                          >  1:299312/1‑251 (MQ=255)
                                                                                                                                                                                    aGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCCCCAGCACAACTGCCTGTCCTGGTGCCCGCAGTCCTGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTTCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGgtttc                                                         <  1:259228/251‑1 (MQ=255)
                                                                                                                                                                                          gCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGgtttcgtttc                                                    <  1:17875/250‑1 (MQ=255)
                                                                                                                                                                                                   tgcGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTAt                                            <  1:92919/249‑1 (MQ=255)
                                                                                                                                                                                                                  cGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACctct                           <  1:3344/251‑1 (MQ=255)
                                                                                                                                                                                                                         gggTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAg                    <  2:299312/251‑1 (MQ=255)
                                                                                                                                                                                                                            ttACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTg                 >  2:11143/1‑251 (MQ=255)
                                                                                                                                                                                                                            ttACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTg                 >  1:179748/1‑251 (MQ=255)
                                                                                                                                                                                                                               ctacCCACCAGCACAACGGCCTGTCCTGTTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATGATCTTTGCTGCCATCAGAGCGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAgg              >  2:380481/1‑251 (MQ=255)
                                                                                                                                                                                                                                     accaGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCAcc        >  2:146919/1‑251 (MQ=255)
                                                                                                                                                                                                                                           acaACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTGCCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTTCGTAGATCGAATGGGATGCGCCTATTTTGTCAGATTCGGGGCGCGCTGTCTGGTTTCGTTTCAGATTTATTGCGTCAAGCGACCTCCGGCGAAGGTGATGCGTTCCGTCGCt  <  2:178965/251‑1 (MQ=255)
                                                                                                                                                                                                                                           acaACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCt  <  1:204762/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
GTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCAGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCT  >  NC_000913/557391‑557876

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