Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 863,376 A→T 6.3% A170A (GCA→GCT hemH → Ferrochelatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451863,3760AT 6.3% 218.5 / 3.4 64A170A (GCA→GCThemHFerrochelatase
Reads supporting (aligned to +/- strand):  ref base A (35/24);  new base T (2/2);  total (38/26)
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.66e-01

GCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGAAGTAAGTTCGGTTTGCTATCGCGGGCAATGCCTGGCTACGGCCAAAGCCTTTGCGCAGAGCATG  >  AE015451/863152‑863589
                                                                                                                                                                                                                                |                                                                                                                                                                                                                     
gCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCg                                                                                                                                                                                                            >  1:1063660/1‑236 (MQ=255)
  cATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGc                                                                                                                                                                                                          <  1:782553/236‑1 (MQ=255)
   aTGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCg                                                                                                                                                                                                            >  2:973055/1‑233 (MQ=255)
           cGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGcagccagc                                                                                                                                                                                                                 <  2:476092/220‑1 (MQ=255)
                         tGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCGGGAGCAGGGGTACGac                                                                                                                                                                                   >  1:306469/1‑236 (MQ=255)
                               aTGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGccc                                                                                                                                                                                                             <  2:632399/204‑1 (MQ=255)
                                  tgtTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAgcgc                                                                                                                                                                                                               >  2:795893/1‑199 (MQ=255)
                                        cgcgGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACt                                                                                                                                                                        >  2:1068701/1‑232 (MQ=255)
                                          cgGCTGGCCGCGCAAGGTGTGCGCAAGGTTACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCGCTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGcag                                                                                                                                                                                                                      <  2:142506/184‑1 (MQ=255)
                                             cTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTa                                                                                                                                                                                       <  2:1100390/212‑1 (MQ=255)
                                             cTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTa                                                                                                                                                                                       >  1:1100391/1‑212 (MQ=255)
                                                           tgtgCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATGGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGcctcgatgta                                                                                                                                                                                                      <  1:611737/183‑9 (MQ=255)
                                                                                               gTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAg                                                                                                              >  1:676869/1‑235 (MQ=255)
                                                                                                ttCGCCGACAGCACGGTGCCCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTgaag                                                                                                                                      <  1:120766/210‑1 (MQ=255)
                                                                                                             cGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTg                                                                                                                                                        >  1:314927/1‑179 (MQ=255)
                                                                                                             cGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTcc                                                                                                                           >  2:129287/1‑208 (MQ=255)
                                                                                                             cGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTgcgcgc                                                                                                >  1:88702/1‑235 (MQ=255)
                                                                                                             cGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTgcgc                                                                                                  >  2:693072/1‑233 (MQ=255)
                                                                                                                tGACCACTGTGATCGAGCAGGCAAAACGGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGcctcgatg                                                                                                                                                                                                        <  2:932785/128‑7 (MQ=255)
                                                                                                                tGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGGGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGcc                                                                                            <  1:449904/236‑1 (MQ=255)
                                                                                                                 gACCCCTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCCCCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTa                                                                                                                                                                                       <  1:108503/144‑1 (MQ=255)
                                                                                                                  aCCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGa                                                                                          <  1:992390/236‑1 (MQ=255)
                                                                                                                    cACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGAt                                                                                         <  1:475067/235‑1 (MQ=255)
                                                                                                                       tgtgATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATtgc                                                                                      <  2:676868/235‑1 (MQ=255)
                                                                                                                        gtgATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACg                                                                                                                                                >  1:277993/1‑176 (MQ=255)
                                                                                                                         tgATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCTGCCAGCGCCCgg                                                                                                                                                                                                           >  2:216092/1‑116 (MQ=255)
                                                                                                                          gATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATtgct                                                                                     >  2:489432/1‑233 (MQ=255)
                                                                                                                           aTCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGa                                                                                                                         <  2:605472/196‑1 (MQ=255)
                                                                                                                           aTCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCAc                                                                                                           >  2:948151/1‑210 (MQ=255)
                                                                                                                           aTCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATACGCGTGCTGCAGCCGCTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGtt                                                                                                                                                          >  1:829206/1‑163 (MQ=255)
                                                                                                                             cGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTgg                                                                                                                                                                                                 <  2:531379/122‑1 (MQ=255)
                                                                                                                             cGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTgg                                                                                                                                                                                                 >  1:531380/1‑122 (MQ=255)
                                                                                                                                 cAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGc                                                                           >  1:972375/1‑236 (MQ=255)
                                                                                                                                    gCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCTGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAAc                                                                        >  2:146123/1‑236 (MQ=255)
                                                                                                                                    gCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTgaa                                                                                                                                       >  2:324536/1‑173 (MQ=255)
                                                                                                                                                 tCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTgcgcgc                                                                                                >  2:958534/1‑199 (MQ=255)
                                                                                                                                                              gTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGtt                                                                                                                              >  2:581850/1‑156 (MQ=255)
                                                                                                                                                               gtACCGTTGCAGATGCGCGTGCTACAGCCGTACTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCt                                                                                                      >  2:1042795/2‑179 (MQ=255)
                                                                                                                                                                 aCCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGAAGTAAGTTCGGTTTGCTATcgcg                                           <  1:972946/236‑1 (MQ=255)
                                                                                                                                                                         aGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCCGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGACCGGCACCTGAAGAAACTGTTTc                                                                                                                            >  2:1033482/1‑147 (MQ=255)
                                                                                                                                                                         aGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTAACTGGAGCAGGGCTACGACCAGCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTgaa                                                                                                                                       >  1:175439/1‑136 (MQ=255)
                                                                                                                                                                                        aGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGAAGTAAGTTCGGTTTGCTATCGCGGGCAATGcc                                  >  1:196326/1‑222 (MQ=255)
                                                                                                                                                                                          ccGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCTGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATtgc                                                                                      <  2:494368/168‑1 (MQ=255)
                                                                                                                                                                                          ccGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGAAGTAAGTTCGGTTTGCTATCGCGGGCAATGCCTGGCTACGGCCa                      >  2:580058/1‑232 (MQ=255)
                                                                                                                                                                                           cGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTa                                                                                                                                                                          <  1:485690/83‑1 (MQ=255)
                                                                                                                                                                                               cTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGAAGTAAGTTCGGTTTGCTATCGCGGGCAATGCCTGGCTACGGCCAAAGCCttt              <  1:928896/235‑1 (MQ=255)
                                                                                                                                                                                                tATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTGTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGACGTAAGTTCgg                                                       >  2:354998/1‑193 (MQ=255)
                                                                                                                                                                                                  tGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGAAGTAAGTTCGGTTTg                                                   >  2:663279/1‑195 (MQ=255)
                                                                                                                                                                                                  tGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGAAGTAAGTTCGGTTTGCTATCGCGGGCAATGCCTGGCTACGGCCAAAGCCTTTgc            <  1:61870/234‑1 (MQ=255)
                                                                                                                                                                                                   gAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGAAGTAAGTTCGGTTTGCTATCGCGGGCAATGCCTGGCTACGGCCAAAGCCTTTg             <  1:496226/232‑1 (MQ=255)
                                                                                                                                                                                                          cGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGAAGTAAGTTCGGTTTGCTATCGCGGGCAATGCCTGGCTACGGCCAAAGCCTTTGCGCAGAGc     >  2:974621/1‑233 (MQ=255)
                                                                                                                                                                                                          cGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGAAGTAAGTTCGGTTTGCTATCGCGGGCAATGCCTGGCTACGGCCAAAGCCTTTGCGCAGAGCATg  >  1:131325/1‑236 (MQ=255)
                                                                                                                                                                                                            gCGTCCATCGAGGCGCTGGCAGCCAGCGCCCGGCCGt                                                                                                                                                                                                       <  2:297839/37‑1 (MQ=255)
                                                                                                                                                                                                            gCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGt                                                                                                                                                                                                       >  1:297840/1‑37 (MQ=255)
                                                                                                                                                                                                            gCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGACGAGACTGt                                                                                                                               >  2:720500/1‑109 (MQ=255)
                                                                                                                                                                                                             cGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGAAGTAAGTTCGGTTTGCTATCGCGGGCaa                                      <  1:955725/197‑1 (MQ=255)
                                                                                                                                                                                                             cGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGAAGTAAGTTCGGTTTGCTATCGCGGGCaa                                      >  2:955724/1‑197 (MQ=255)
                                                                                                                                                                                                                  agcGAGGCGCGGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTgc                                                                                                    <  2:118533/128‑1 (MQ=255)
                                                                                                                                                                                                                     cgggCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCg                                                                                           >  1:363526/3‑136 (MQ=255)
                                                                                                                                                                                                                      aGGCGCTGGCTGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGAAGTAAGTTCGGTTTGCTATCGCgg                                          <  1:146123/184‑1 (MQ=255)
                                                                                                                                                                                                                         ggcTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTa                                                                                                                                                                                       <  2:358449/39‑1 (MQ=255)
                                                                                                                                                                                                                          gcTGGCAGCCAGCGCCCGGCCGTACCTGGGGCAGGGtt                                                                                                                                                                                        >  2:204439/1‑38 (MQ=255)
                                                                                                                                                                                                                          gcTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGac                                                                 >  2:88637/1‑156 (MQ=255)
                                                                                                                                                                                                                          gcTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGAAGTAAGTTCGGTTTGCTATCGCGGGCAATGCCTGGCTACGGCCAAAGc                  >  2:922796/1‑204 (MQ=255)
                                                                                                                                                                                                                             ggCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATg                                                                                                                                                             <  2:277992/62‑1 (MQ=255)
                                                                                                                                                                                                                              gCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGAAGTAAGTTCGGTTTGCTATCGCGGGCAATGCCTGGCTACGGCCAAAGCCTTTgc            <  1:929083/206‑1 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                     
GCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGAAGTAAGTTCGGTTTGCTATCGCGGGCAATGCCTGGCTACGGCCAAAGCCTTTGCGCAGAGCATG  >  AE015451/863152‑863589

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