Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,989,202 A→T 5.6% C209* (TGT→TGA hemD ← uroporphyrinogen III synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,989,2020AT 5.6% 319.4 / 5.1 103C209* (TGT→TGAhemDuroporphyrinogen III synthase
Reads supporting (aligned to +/- strand):  ref base A (41/56);  new base T (3/3);  total (44/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 104 total aligned reads displayed.

CTTGTTCCGTCATTATGGCTTCCTGTTATGAGAGTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGCTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCAT  >  NC_000913/3989054‑3989349
                                                                                                                                                    |                                                                                                                                                   
cTTGTTCCGTCATTATGGCTTCCTGTTATGAGAGTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGa                                                                                                                                                     <  2:1309230/149‑1 (MQ=255)
 ttGTTCCGTCATTATGGCTTCCTGTTATGAGAGTTATTGTAATGCCCGTAAAAGCGCAGCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGAc                                                                                                                                                    <  1:1551473/149‑1 (MQ=255)
  tGTTCCGTCATTATGGCTTCCTGTTATGAGTGTTTTTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCTGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCTCCAAACGCTCACTGTCGACTATTGGTCGTCa                                                                                                                                                   <  1:733222/149‑1 (MQ=255)
   gTTCCGTCATTATGGCTTCCTGTTATGAGAGTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAg                                                                                                                                                  >  1:1215548/1‑149 (MQ=255)
    ttCCGTCATTATGGCTTCCTGTTATGAGAGTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGTCGACTATTGGTCGTCAgt                                                                                                                                                 <  1:1534578/149‑1 (MQ=255)
       cGTCATTATGGCTTCCTGTTATGAGAGTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTgta                                                                                                                                              >  2:1214909/1‑149 (MQ=255)
       cGTCATTATGGCTTCCTGTTATGAGAGTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTgta                                                                                                                                              <  2:1259558/149‑1 (MQ=255)
       cGTCATTATGGCTTCCTGTTATGAGAGTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAAATGTCGACAgtgtt                                                                                                                                              >  2:1486257/1‑148 (MQ=255)
        gTCATTATGGCTTCCTGTTATGAGAGTTATTGTATTGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTTCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGTCGTCTATTGGTCGTCAgtgagc                                                                                                                                             <  2:607151/149‑4 (MQ=255)
            ttATGGCTTCCTGTTATGTGAGTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAAc                                                                                                                                         <  2:1227188/149‑1 (MQ=255)
            ttATGGCTTCCTGTTATGAGAGTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAAc                                                                                                                                         <  1:606852/149‑1 (MQ=255)
            ttATGGCTTCCTGTTATGAGAGTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAAc                                                                                                                                         >  1:720384/1‑149 (MQ=255)
            ttATGGCTTCCTGTTATGAGAGTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAAc                                                                                                                                         <  2:1541594/149‑1 (MQ=255)
            ttATGGCTTCCTGTTATGAGAGTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAAc                                                                                                                                         <  2:1527941/149‑1 (MQ=255)
               tGGCTTCCTGTTATGAGAGTTATTGTAATGCCCGTAAAAGCGCAGCGGGGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAg                                                                                                                                      <  1:1195393/149‑1 (MQ=255)
                   ttCCTGTTATGAGAGTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTAGACAGTGTAGTAACCAGTGca                                                                                                                                  >  2:791453/1‑148 (MQ=255)
                     ccTGTTATGAGAGTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCa                                                                                                                                <  1:197283/149‑1 (MQ=255)
                      cTGTTATGAGAGTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCAc                                                                                                                               <  1:1421036/149‑1 (MQ=255)
                       tGTTATGAGAGTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGTGTTTCGCCAAACGCTCACTGACGACCAATAGTCCACAGTGTAGTACCCAGTGCtgccg                                                                                                                              >  2:74273/1‑145 (MQ=255)
                        gTTATGAGAGTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGa                                                                                                                             <  1:1226151/149‑1 (MQ=255)
                        gTTATGAGAGTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGa                                                                                                                             <  1:934108/149‑1 (MQ=255)
                         ttATGAGAGTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAgtg                                                                                                                                                >  1:1481364/1‑129 (MQ=255)
                         ttATGAGAGTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAgtg                                                                                                                                                <  2:1481364/129‑1 (MQ=255)
                            tGAGAGTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCATGCCAGCCAAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGAACAATAGTCGACAGTGTAGTAACCAGTGATCACGATagc                                                                                                                         >  2:1411216/1‑147 (MQ=255)
                             gagaGTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCa                                                                                                                        <  1:1522331/149‑1 (MQ=255)
                             gagaGTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCa                                                                                                                        >  2:320879/1‑149 (MQ=255)
                              agagTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCAt                                                                                                                       <  2:457791/149‑1 (MQ=255)
                                agTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCTGCGTTATCGGCGACGTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCAAtg                                                                                                                     >  1:188127/1‑149 (MQ=255)
                                agTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATtg                                                                                                                     >  1:787948/1‑149 (MQ=255)
                                agTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATtg                                                                                                                     >  2:971420/1‑149 (MQ=255)
                                 gTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGAGTAGTAACCAGTGCTCACGATACCATtgt                                                                                                                    >  2:782751/1‑149 (MQ=255)
                                  ttATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATtgtg                                                                                                                   <  2:1323243/149‑1 (MQ=255)
                                  ttATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATtg                                                                                                                     >  1:265090/1‑147 (MQ=255)
                                  ttATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATtg                                                                                                                     <  2:265090/147‑1 (MQ=255)
                                  ttATTGTAATGCCCGTAAAAGCGAATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATtgtg                                                                                                                   <  1:46973/149‑1 (MQ=255)
                                    aTTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTggg                                                                                                                 >  1:34665/1‑149 (MQ=255)
                                        tAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCAtt                                                                                                                      >  2:898448/1‑140 (MQ=255)
                                            gCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGa                                                                                                         <  1:742412/149‑1 (MQ=255)
                                             cccGTAAATGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTTGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGAc                                                                                                        <  1:162527/149‑1 (MQ=255)
                                             cccGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCAAGATACCATTGTGGGATCAGCGAc                                                                                                        <  2:393354/149‑1 (MQ=255)
                                             cccGTAAAAGCGCATCGTTGTCAGAGTTATCGGCGACATTAATGTCTTGCCAGCCAAGTTCACGTGCGAGTTTCTCAAAACGCTCAATGACTACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGAc                                                                                                        >  2:33405/1‑149 (MQ=255)
                                               cGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCa                                                                                                      >  1:984912/1‑149 (MQ=255)
                                                    aaGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGtt                                                                                                 <  1:1131680/149‑1 (MQ=255)
                                                    aaGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGtt                                                                                                 <  2:787948/149‑1 (MQ=255)
                                                      gcgcATCGTTGTCAGCTTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTTCTCACGATACCATTGTGGGATCAGCGACCAGAGTtgc                                                                                               <  1:330496/149‑1 (MQ=255)
                                                       cgcATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTtgct                                                                                              >  1:782340/1‑149 (MQ=255)
                                                       cgcATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTtgct                                                                                              <  2:20859/149‑1 (MQ=255)
                                                       cgcATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTtgct                                                                                              <  2:34665/149‑1 (MQ=255)
                                                        gcATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTtgctg                                                                                             <  1:782751/149‑1 (MQ=255)
                                                           tCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCaa                                                                                          >  1:1188195/1‑149 (MQ=255)
                                                                ggcAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACTACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACAttt                                                                                     <  2:461423/147‑1 (MQ=255)
                                                                gTCAGCGTTTTCTGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTGGCTGCAACAttt                                                                                     <  2:1098019/149‑1 (MQ=255)
                                                                gTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACAttt                                                                                     <  2:1529791/149‑1 (MQ=255)
                                                                 tCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTtgct                                                                                              <  2:1149891/139‑1 (MQ=255)
                                                                 tCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTtgct                                                                                              >  1:1149891/1‑139 (MQ=255)
                                                                 tCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTc                                                                                    >  2:720974/1‑149 (MQ=255)
                                                                 tCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTc                                                                                    >  2:54325/1‑149 (MQ=255)
                                                                        tATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTg                                                                             <  1:1486257/149‑1 (MQ=255)
                                                                        tATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTg                                                                             <  2:1391760/149‑1 (MQ=255)
                                                                         aTCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCACTAGTCGACAGTGTAGTAACCAGTGCTCACGATACCCTTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTgg                                                                            <  1:1177241/149‑1 (MQ=255)
                                                                             gCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTaac                                                                        >  1:585555/1‑149 (MQ=255)
                                                                               gACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTaacaa                                                                      <  1:1050627/149‑1 (MQ=255)
                                                                               gACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTaacaa                                                                      <  1:720974/149‑1 (MQ=255)
                                                                               gACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTaacaa                                                                      >  2:1359607/1‑149 (MQ=255)
                                                                                  cTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGa                                                                   <  1:320879/149‑1 (MQ=255)
                                                                                  cTTAATGGCTTGCCAGCCCAGTTCCCGGGCGAGGTTCGCCAAACGCTCACAGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGa                                                                   <  1:66361/149‑1 (MQ=255)
                                                                                   ttAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGAc                                                                  <  2:720384/149‑1 (MQ=255)
                                                                                   ttAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGAc                                                                  >  1:1388331/1‑149 (MQ=255)
                                                                                       tGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATc                                                              >  1:802160/1‑149 (MQ=255)
                                                                                         tCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGt                                                            >  2:552370/1‑149 (MQ=255)
                                                                                         tCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGt                                                            <  2:781790/149‑1 (MQ=255)
                                                                                                     aGTTCCCGCGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTGACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTTACAACGACAATCGTCACCTCGCGGGc                                                >  2:264161/1‑149 (MQ=255)
                                                                                                       ttCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGCtt                                              <  1:54325/149‑1 (MQ=255)
                                                                                                           cGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACAGATGGTAACAACGACCATCGTCCCCTCGCGGTCTTGCCa                                          >  2:1254297/1‑149 (MQ=255)
                                                                                                            gggCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCACTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGCATGCCAg                                         >  2:1599048/1‑149 (MQ=255)
                                                                                                             ggCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGCTTGCCAgc                                        <  2:1147479/149‑1 (MQ=255)
                                                                                                                gAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCTGGCAAGCCcgcgag                                     <  1:733586/149‑3 (MQ=255)
                                                                                                                gAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGCTTGCCAGCGCa                                     >  1:628257/1‑149 (MQ=255)
                                                                                                                 aGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGCTTGCCAGCGCAt                                    <  1:971420/149‑1 (MQ=255)
                                                                                                                  gTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGCTTGCCAGCGCATc                                   >  1:529486/1‑149 (MQ=255)
                                                                                                                    ttCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGCTTGCCAGCGCATCGc                                 <  2:782340/149‑1 (MQ=255)
                                                                                                                       gCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGCTTGCCAGCGCATCGcttc                              <  1:1271348/149‑1 (MQ=255)
                                                                                                                        ccAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGCTTGCCAGCGCATCGcttct                             >  2:489033/1‑149 (MQ=255)
                                                                                                                          gttactacACTGTCGACTATTGGTCGTCAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACCACGACCATCGTCACCTCGCGGGCTTGCCAGCGCATCGCttcttc                           >  1:419968/8‑149 (MQ=255)
                                                                                                                            aCGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGCTTGCCAGCGCATCGCTTCTTCTg                         <  2:585555/149‑1 (MQ=255)
                                                                                                                              gCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGCTTGCCAGCGCATCGCTTCTTCTGCa                       <  1:1254297/149‑1 (MQ=255)
                                                                                                                              gCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGATGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGCTTGCCAGCGCATCGCTTCTTCTGCa                       <  1:33405/149‑1 (MQ=255)
                                                                                                                               tacACTGTCGACTATTGGTCGTCAGAGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGAACATCGTCACCTCGCGGGCTTGCCAGCGCATCGCTTCTTCTGCAc                      >  1:1013320/3‑149 (MQ=255)
                                                                                                                               cTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGCTTGCCAGCGCATCGCTTCTTCTGCAc                      >  2:1529940/1‑149 (MQ=255)
                                                                                                                                 cACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGCTTGCCAGCGCATCGCTTCTTCTGCACCa                    <  1:1214909/149‑1 (MQ=255)
                                                                                                                                 cACTGACGACAAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGCTTGCCAGCGCATCGCTTCTTCTGCACCa                    <  1:161543/149‑1 (MQ=255)
                                                                                                                                  aCTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCa                                                          >  2:250074/1‑110 (MQ=255)
                                                                                                                                  aCTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCa                                                          <  1:250074/110‑1 (MQ=255)
                                                                                                                                   cTGTCGACTATTGGTCGTCAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACAGCTGGTAACAACGACCATCGTCACCTCGCGGGCTTGCCAGCGCATCGCTTCTTCTGCACCATc                  >  1:1579099/1‑149 (MQ=255)
                                                                                                                                    tGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGCTTGCCAGCGCATCGCTTCTTCTGCACCATCg                 >  1:613459/1‑149 (MQ=255)
                                                                                                                                      acgacCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGCTTGCCAGCGCATCGCTTCTTCTGCACCATCGTa               >  2:1182186/1‑149 (MQ=255)
                                                                                                                                      acgacCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGCTTGCCAGCGCATAGCTTCTTCTGCACCATCGTa               >  2:1423103/1‑149 (MQ=255)
                                                                                                                                       cgacCATTAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCTCTGGTAACAACGACCATCGTCACCTCGCGGGCTTGCCAGCGCATCGCTTCTTCTGCACCATCGTaa              <  1:1411216/149‑1 (MQ=255)
                                                                                                                                       cgacCAATAGTCGACTGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGCTTGCCAGCGCATCGCTTCTTCTGCACCATCGTaa              <  1:791453/149‑1 (MQ=255)
                                                                                                                                                   gACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGCTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCAt  <  2:1565133/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTTGTTCCGTCATTATGGCTTCCTGTTATGAGAGTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCGACCTTAATGTCTTGCCAGCCCAGTTCCCGGGCGAGTTTCGCCAAACGCTCACTGACGACCAATAGTCGACAGTGTAGTAACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGCTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCAT  >  NC_000913/3989054‑3989349

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