Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,617,616 T→C A98A (GCT→GCC deoC → deoxyribose‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,617,6160TC96.4% 166.7 / ‑2.3 56A98A (GCT→GCCdeoCdeoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (25/29);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.17e-01

CTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCC  >  NC_000913/4617371‑4617858
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        
cTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCcgatg                                                                                                                                                                                                                                                      >  1:593478/1‑250 (MQ=255)
 catctccccTGAATGACGACGACACCGACGAAAAATAAGTTGCTCTTGGTCATCAGCCCAAACCCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGaa                                                                                                                                                                                                                                                    <  1:39305/245‑1 (MQ=255)
      accCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGaa                                                                                                                                                                                                                                                    >  2:39305/1‑246 (MQ=255)
            aaTGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGtggg                                                                                                                                                                                                                                         >  1:219543/1‑248 (MQ=255)
                              tagaAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGTAAATCCGTATTGCTACGGTAACCAACTTCCAACACGGTAACGATGACCTCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGAAGTGGAGTTCCCGTAACGCGCGCt                                                                                                                                                                                                                       >  2:509213/2‑251 (MQ=255)
                                gaAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTg                                                                                                                                                                                                                      >  1:244064/1‑250 (MQ=255)
                                    agtaattgctctgagtCATCAGGCCAAAACCCCGGTCGGCAATACCACCGCTATCTGTCTCTATCCTCGCTTTATCCCGATAGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGCAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGGTGACATCGAAATCGCGCGGGCAGCAACCCGTGCGGCAATCTCCTACGGTGACGATGAAGTATACGTGGTGGTCCCGAACCGCGCGCTGacgg                                                                                                                                                                                                                  >  2:219610/15‑248 (MQ=255)
                                          tgctcggtgtCATCAGCCCAAACCCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGt                                                                                                                                                                                                             <  1:281345/244‑1 (MQ=255)
                                           gctcTGTGTCACCAGGCCAAAACCCCGGTCGTCAATTCCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTa                                                                                                                                                                                                            <  1:255707/247‑1 (MQ=255)
                                              cTGTGACATCAGCCGAAACCCCCGGTAGACAATACCGCCGCTAGCTGTATCTATCCTCGCCTTCTCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCGACATACCACACGGTTACGATGCAATCGAAATCGCGCTGACAGAAACCCGTACGGCAATCGCATTCGGGGCCGATGACGTTGATGTGGTGTtaccgtaccgaccgctaccggcggggtacg                                                                                                                                                                                                         >  2:124547/1‑230 (MQ=255)
                                                  gtCATCAGGCCAAAACCCCGGTCGGCAAATCCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAg                                                                                                                                                                                                    <  1:630271/251‑1 (MQ=255)
                                                  gtCAACAGGCAAAACCCCCGTTCGGCATAACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAg                                                                                                                                                                                                    <  1:224296/251‑1 (MQ=255)
                                                     atcaGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGa                                                                                                                                                                                                                     >  1:380394/1‑231 (MQ=255)
                                                       cggtccacacccccGGTCGGAAATCCCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAggttgg                                                                                                                                                                                               <  2:146221/239‑1 (MQ=255)
                                                               acccccgcggggaaacaCCGCCGCGAGCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGAcc                                                                                                                                                                                      <  2:341916/235‑1 (MQ=255)
                                                                         ggCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGc                                                                                                                                                                            >  1:103645/1‑251 (MQ=255)
                                                                          gaaaTACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGcc                                                                                                                                                                           <  1:211407/249‑1 (MQ=255)
                                                                                ccgccgCTAACTGTATCTAACCTCGCCTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCGCCCCGTAAATCCGTACTGCTACGTGAACCAACTTCCCACGCGGTAACCATGACATCGAAATCGCGCTGGCAGAAACCCGCGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGAGGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTaaa                                                                                                                                                                     <  2:275221/251‑1 (MQ=255)
                                                                                 cgccgcTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGtca                                                                                                                                                                      >  1:557682/1‑247 (MQ=255)
                                                                                   ccgcTATCGGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAagag                                                                                                                                                                  <  2:219543/251‑1 (MQ=255)
                                                                                            gTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGTAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGAGCCGATGCAGTTGACGTGGTGTTCACGTCCAGCGCGCTGATGTCGGGAAACGCGCAGGGGGGGTTTTACCTGGTGAAAGACTGTCAAGAGGCTTGcgcg                                                                                                                                                         >  2:483050/1‑251 (MQ=255)
                                                                                            gTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGcgcg                                                                                                                                                         <  1:395715/251‑1 (MQ=255)
                                                                                              atctatcCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGGGCAGGTTGGTTTTTGCCCGGGTAAAGCCTGAAAAGAgggtttgcgggg                                                                                                                                                       >  2:232742/1‑241 (MQ=255)
                                                                                                ctatcCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGGGAAAGCCTGTAAAGAGGCTTTCGCGGc                                                                                                                                                      >  1:483978/1‑249 (MQ=255)
                                                                                                        gCTTTTTCCCGATTGCTCGCAAACCACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTAAGGTCACCAAATTCCCACACGGTCACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGAAATCGCCTACGGAGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACCAGCCGGCTGGTTTTGACCTGGTGCAAGCCAGTACAGAGGCTTgcgggggcgcggatgt                                                                                                                                              >  2:612311/1‑237 (MQ=255)
                                                                                                             aTCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGGAACGAGCCGGTTGGGTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTActgct                                                                                                                                       >  1:427582/1‑251 (MQ=255)
                                                                                                                         cGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGatca                                                                                                                            <  1:36444/250‑1 (MQ=255)
                                                                                                                           cAAAACACTGCAAGCGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAAATTCCCACACGGTAACGATGACATCGTAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGCTGGCGGGAAACGAGCAGGTAGGGTTTGACCTGGTGAAAGCCTGTGAAGAGGATTGGGCGGCAGGGGATGGACTGCTGCAGATGATCATCg                                                                                                                         >  1:370434/1‑251 (MQ=255)
                                                                                                                           cAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGGTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTg                                                                                                                                                             >  1:196722/1‑216 (MQ=255)
                                                                                                                              aacactGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGaaa                                                                                                                      <  2:548546/247‑1 (MQ=255)
                                                                                                                              aacactGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGGTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGAGAGCATCGaa                                                                                                                       >  1:622815/5‑250 (MQ=255)
                                                                                                                                    gAAAGAGCAGTGCACCCCGGAAATCCGTAATGCTACGGTAACCAACTCCCCACACGGTAACGATTACATCGAAATCTGGCTGGCGGAACCCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAAGGTACTGCTTAAAGTGATCATCGAAACCGGCg                                                                                                                <  2:597761/251‑1 (MQ=255)
                                                                                                                                    gAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGGTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCgg                                                                                                                  <  1:612311/249‑1 (MQ=255)
                                                                                                                                      aaGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGGATGTACTGCTGAAAGGGATCATCGAAAACGGCGa                                                                                                               >  1:416419/1‑250 (MQ=255)
                                                                                                                                      aaGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTGGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGaa                                                                                                              <  2:483978/251‑1 (MQ=255)
                                                                                                                                       agcgCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAAAGAGCAGGGTGGTTTTTACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGGATGGACTGCTGAAAGTGATCATCAAAACCGGCGAAc                                                                                                             >  2:294121/4‑251 (MQ=255)
                                                                                                                                                        aaaTCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTg                                                                                           <  2:244064/251‑1 (MQ=255)
                                                                                                                                                                       cGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGGGAAAGCCTGTAAAGAGGCTTGCGCGGCAGGGGATGTAATGCTGAAAGTGATCATCGAAACCGGCGCACCGAAAAAGGAAGCGCTGagcccgaaagaggat                                                                            >  2:312757/1‑237 (MQ=255)
                                                                                                                                                                        ggTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTAACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGCGCACGTTGTGTTTGACCTGGt                                                                                                                                                                                  >  1:232235/1‑150 (MQ=255)
                                                                                                                                                                               aaCTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCCAAACCGGCGAACTTGAAGACGAAGCGCGGATCCGTAAAGCGTCTGaaaattc                                                                   >  1:296277/1‑247 (MQ=255)
                                                                                                                                                                                        cacGGTAATGATGACATCGAAATCGCGCTGGCAGAAACCAGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGc                                                         <  1:500037/251‑1 (MQ=255)
                                                                                                                                                                                         aaggTAACGATGACATCGAAATCGCGTTGGCAGAATCCCGTGCGGAAATCGCCTACGGTGCCGATGAAGGTGACGTGGTGTGCCCGTACCGCGCGCTGAGGGAGGGTAACGAGCAGGTAGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGAGCAACGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGGCTGAAATCTCCAGCAAAgag                                                        <  2:43902/249‑3 (MQ=255)
                                                                                                                                                                                              atcgaggactgcgaaatcgCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGAAGTGGTTTTCCCGTACCGCGCGCTGAGGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGACTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGCCGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTg                                                    <  2:92871/235‑1 (MQ=255)
                                                                                                                                                                                              aaCGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGc                                                   <  1:116409/251‑1 (MQ=255)
                                                                                                                                                                                                 gatgacatcgaaatcGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGCTGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAAAGACTTGAGCGGACGGGGAATTACTTCTGAAAGGGATGATCGAAACCGGCGAACTGACAGGCGCAGAGCCGAGCCGGAAAACGGCTGAAATCTCCATCCAAGACGGTGCggg                                                >  2:578196/4‑250 (MQ=255)
                                                                                                                                                                                                  atgacatcgaaatcGCGCTGGCAGAAGCCCGGGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGAc                                               <  2:367427/249‑1 (MQ=255)
                                                                                                                                                                                                   tgacatcgaaatcGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTAATGCTGAAAGGGGTCATCGAAACCGGCGAACTGAAAGACGCAGCGCTGGTCCGTAAAGCGTCTGAAATCTCCATCCAAACGGGGGCgggct                                              >  1:597892/2‑248 (MQ=255)
                                                                                                                                                                                                                           aaCCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCCTCAAAACCTCTACCGGTaaa                      >  1:248797/1‑250 (MQ=255)
                                                                                                                                                                                                                              ccGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTgg                  <  1:623554/251‑1 (MQ=255)
                                                                                                                                                                                                                                 tGCGGCAATCGCCTACGGGGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGATAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATTTACTTCTGAAAGAGATCATCCACACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAACGTCTTAAATCTCCATCCAAGCGTGTGCGGACTTCCTCAAAACCTCTAcccgtcaagatgctt               >  1:193063/1‑238 (MQ=255)
                                                                                                                                                                                                                                  gCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCtgt              <  1:635825/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gCGGCAATCGCCTACAGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGATGAGGGCGGGTAACGAGCAGGGTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCtgt              <  2:315042/251‑1 (MQ=255)
                                                                                                                                                                                                                                   cGGCAATCGCCTACGGTGCCGATGAAGATGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGGTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCCGGGCCACGAATTTACTGCTGGAAGGGGTCATTGAAATCGGCGAAATCAAAAACGAAACGCGTGTCCAGAAAGCGTCGTGAAACTCCATCAAAAGCGGGGCGGACCTCCAACAAAACCTTTACGGTTAAGTgggtggg             >  2:523746/1‑246 (MQ=255)
                                                                                                                                                                                                                                      cAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCCCTACCGGTAAAGTGGCTGTGAAc          <  2:144480/251‑1 (MQ=255)
                                                                                                                                                                                                                                      cAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAAcg         <  1:142028/251‑1 (MQ=255)
                                                                                                                                                                                                                                               aCGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGcc  <  1:604532/250‑1 (MQ=255)
                                                                                                                                                                                                                                                cGGTGCCGATGAAGGTTAAGTGGTGTTCCCGTACCGGGCTCTTATGGCGGGGAACGAGGAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCCGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGTGTCTGAAATCTCCATCAAAGCGGGGGCGGACGTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGcc  <  2:416419/249‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        
CTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCC  >  NC_000913/4617371‑4617858

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