Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,046,172 A→T 100% K256N (AAA→AAT hemN → coproporphyrinogen III oxidase, SAM and NAD(P)H dependent, oxygen‑independent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,046,1720AT95.2% 430.0 / 8.8 135K256N (AAA→AAThemNcoproporphyrinogen III oxidase, SAM and NAD(P)H dependent, oxygen‑independent
Reads supporting (aligned to +/- strand):  ref base A (3/4);  new base T (66/62);  total (69/66)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 137 total aligned reads displayed.

CGACCTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGG  >  CP009273/4046028‑4046315
                                                                                                                                                |                                                                                                                                               
cGACCTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGAt                                                                                                                                              >  2:1277442/1‑148 (MQ=255)
cGACCTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGAt                                                                                                                                              >  2:1690360/1‑148 (MQ=255)
cGACCTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCACAAAATCAATGAt                                                                                                                                              >  1:3113649/1‑148 (MQ=255)
 gACCTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAatgatg                                                                                                                                             <  1:3269244/148‑1 (MQ=255)
 gACCTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAatgatg                                                                                                                                             <  2:3108730/148‑1 (MQ=255)
  aCCTGATTTACGGCCTGCCGATACAGACGCCGGAGTGTTTCGCCTTTACCCTGAAACGTGTGGCGGATCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGc                                                                                                                                            <  2:3025268/148‑1 (MQ=255)
  aCCTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGc                                                                                                                                            <  2:1843961/148‑1 (MQ=255)
   ccTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCt                                                                                                                                           >  1:3119587/1‑148 (MQ=255)
   ccTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCt                                                                                                                                           >  1:157566/1‑148 (MQ=255)
         ttACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTg                                                                                                                                     <  2:2662539/148‑1 (MQ=255)
          tACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGc                                                                                                                                    <  1:14274/148‑1 (MQ=255)
          tACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGc                                                                                                                                    <  1:14260/148‑1 (MQ=255)
              gCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAg                                                                                                                                >  1:2704144/1‑148 (MQ=255)
                   ccGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCgc                                                                                                                           <  2:2404949/148‑1 (MQ=255)
                       aaCAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCgcagca                                                                                                                       >  1:425823/1‑148 (MQ=255)
                          aGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCaaaa                                                                                                                    <  2:2699018/148‑1 (MQ=255)
                            aCGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAAc                                                                                                                  >  2:1982726/1‑148 (MQ=255)
                               ccGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCg                                                                                                               <  2:531333/148‑1 (MQ=255)
                               ccGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCg                                                                                                               <  1:1960430/148‑1 (MQ=255)
                               ccGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCg                                                                                                               <  1:291987/148‑1 (MQ=255)
                                 ggAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGat                                                                                                             <  1:3043155/148‑1 (MQ=255)
                                 ggAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGat                                                                                                             <  2:73515/148‑1 (MQ=255)
                                 ggAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGat                                                                                                             <  1:1817715/148‑1 (MQ=255)
                                 ggAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGat                                                                                                             <  1:1567248/148‑1 (MQ=255)
                                 ggAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGat                                                                                                             <  1:146142/148‑1 (MQ=255)
                                  gagaGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGata                                                                                                            >  2:778346/1‑148 (MQ=255)
                                  gagaGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGata                                                                                                            >  2:778351/1‑148 (MQ=255)
                                  gagaGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGata                                                                                                            <  1:737905/148‑1 (MQ=255)
                                     agTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATcc                                                                                                         >  1:1328339/1‑148 (MQ=255)
                                       tttCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTg                                                                                                       <  2:58473/148‑1 (MQ=255)
                                       tttCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTg                                                                                                       >  2:70381/1‑148 (MQ=255)
                                       tttCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGGCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTg                                                                                                       >  1:1199326/1‑148 (MQ=255)
                                            ccTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGa                                                                                                  >  1:2333096/1‑148 (MQ=255)
                                             cTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGaa                                                                                                 <  1:74329/148‑1 (MQ=255)
                                               ttACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAAc                                                                                               >  1:3143511/1‑148 (MQ=255)
                                               ttACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAAc                                                                                               >  2:1956672/1‑148 (MQ=255)
                                               ttACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAAc                                                                                               >  2:1956675/1‑148 (MQ=255)
                                                 aCCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAAcct                                                                                             <  2:2743071/148‑2 (MQ=255)
                                                 aCCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCa                                                                                             <  2:416258/148‑1 (MQ=255)
                                                 aCCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCa                                                                                             <  2:416266/148‑1 (MQ=255)
                                                 aCCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCa                                                                                             <  2:1916079/148‑1 (MQ=255)
                                                  cccTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCAt                                                                                            >  2:358145/1‑148 (MQ=255)
                                                   ccTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATc                                                                                           >  2:3140480/1‑148 (MQ=255)
                                                     tGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGc                                                                                         <  2:2643491/148‑1 (MQ=255)
                                                       aaaCGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCt                                                                                       >  2:2564808/1‑148 (MQ=255)
                                                         aCGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTc                                                                                     >  1:1801840/1‑148 (MQ=255)
                                                                 cGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCa                                                                             >  2:579202/1‑148 (MQ=255)
                                                                  ggAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCaa                                                                            <  1:1798608/148‑1 (MQ=255)
                                                                  ggAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCaa                                                                            <  1:602928/148‑1 (MQ=255)
                                                                     gCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCg                                                                         <  1:2499827/148‑1 (MQ=255)
                                                                      cTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCgg                                                                        >  1:990839/1‑148 (MQ=255)
                                                                      cTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCgg                                                                        >  1:990851/1‑148 (MQ=255)
                                                                       tGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCggg                                                                       <  1:2061551/148‑1 (MQ=255)
                                                                         aaCCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCt                                                                     <  1:255607/148‑1 (MQ=255)
                                                                         aaCCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCt                                                                     <  2:2937604/148‑1 (MQ=255)
                                                                           ccccGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTAt                                                                   <  1:590928/148‑1 (MQ=255)
                                                                            cccGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATc                                                                  <  2:2435321/148‑1 (MQ=255)
                                                                             ccGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCa                                                                 >  1:1139115/1‑148 (MQ=255)
                                                                              cGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTCTCAg                                                                >  1:3208728/1‑148 (MQ=255)
                                                                              cGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAg                                                                <  2:2886020/148‑1 (MQ=255)
                                                                               gATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGt                                                               >  1:337441/1‑148 (MQ=255)
                                                                                      tGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCgg                                                        <  1:214155/148‑1 (MQ=255)
                                                                                      tGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCgg                                                        <  1:1883122/148‑1 (MQ=255)
                                                                                       gAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGt                                                       <  2:1978727/148‑1 (MQ=255)
                                                                                          tgtCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATg                                                    >  2:2103524/1‑148 (MQ=255)
                                                                                           gtCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATgg                                                   >  2:1208469/1‑148 (MQ=255)
                                                                                            tCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGa                                                  <  2:157566/148‑1 (MQ=255)
                                                                                               ttAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCa                                               >  1:1564574/1‑148 (MQ=255)
                                                                                                   cTACGCGCATCTGCCGACCATTTTTGCTGCTCCGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACttt                                           >  2:1234991/1‑148 (MQ=255)
                                                                                                   cTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACttt                                           >  1:1821044/1‑148 (MQ=255)
                                                                                                    tACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCACTCGGGCTATCAGTTTATCGGTATGGATCACTTTg                                          >  2:2412501/1‑148 (MQ=255)
                                                                                                     aCGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGc                                         <  1:761069/148‑1 (MQ=255)
                                                                                                            tCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCg                                  >  2:2536070/1‑148 (MQ=255)
                                                                                                            tCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCg                                  >  1:3099275/1‑148 (MQ=255)
                                                                                                            tCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCg                                  >  2:2536283/1‑148 (MQ=255)
                                                                                                             cTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCAGAGTCCGCAGAAAAAACTCGATATCCTGCAGGAAACAATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCTGTATGGATCACTTTGCCAGTCAgg                                 >  2:961663/1‑148 (MQ=255)
                                                                                                                ccGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAAACATCTCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATg                              <  2:2751708/148‑1 (MQ=255)
                                                                                                                   aCCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACg                           >  2:2202833/1‑148 (MQ=255)
                                                                                                                    ccATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGa                          >  1:487350/1‑148 (MQ=255)
                                                                                                                    ccATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGa                          <  1:2498589/148‑1 (MQ=255)
                                                                                                                    ccATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGa                          >  2:549351/1‑148 (MQ=255)
                                                                                                                    ccATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGa                          >  2:548329/1‑148 (MQ=255)
                                                                                                                     cATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAg                         >  1:2581483/1‑148 (MQ=255)
                                                                                                                     cATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAg                         >  1:925478/1‑148 (MQ=255)
                                                                                                                      aTTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGc                        >  2:2040705/1‑148 (MQ=255)
                                                                                                                       tttttGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCt                       >  1:3169661/1‑148 (MQ=255)
                                                                                                                        ttttGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTg                      <  2:3088195/148‑1 (MQ=255)
                                                                                                                        ttttGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTg                      <  1:1539879/148‑1 (MQ=255)
                                                                                                                         tttGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTgg                     <  2:425823/148‑1 (MQ=255)
                                                                                                                         tttGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTgg                     <  1:2404700/148‑1 (MQ=255)
                                                                                                                           tgctgcTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCg                   >  2:1843805/1‑148 (MQ=255)
                                                                                                                               gctCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTgg               >  1:2355822/1‑148 (MQ=255)
                                                                                                                               gctCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTgg               >  1:380350/1‑148 (MQ=255)
                                                                                                                                ctCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCTATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGc              >  2:1401116/1‑148 (MQ=255)
                                                                                                                                    gcgcAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAg          <  2:990839/148‑1 (MQ=255)
                                                                                                                                          aaTCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGaa    >  1:1656158/1‑148 (MQ=255)
                                                                                                                                            tCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAgg  >  2:740295/1‑148 (MQ=255)
                                                                                                                                            tCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAgg  >  2:739431/1‑148 (MQ=255)
                                                                                                                                            tCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAgg  <  1:2565357/148‑1 (MQ=255)
                                                                                                                                            tCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAgg  <  2:337500/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
CGACCTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGG  >  CP009273/4046028‑4046315

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