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,1720AT97.4% 484.0 / 2.6 153K256N (AAA→AAThemNcoproporphyrinogen III oxidase, SAM and NAD(P)H dependent, oxygen‑independent
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base T (75/73);  total (78/75)
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.58e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 153 total aligned reads displayed.

TCGACCTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGC  >  CP009273/4046027‑4046316
                                                                                                                                                 |                                                                                                                                                
tCGACCTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGa                                                                                                                                                <  1:621910/148‑1 (MQ=255)
tCGACCTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGa                                                                                                                                                >  1:2744787/1‑148 (MQ=255)
    ccTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCt                                                                                                                                            >  1:2036693/1‑148 (MQ=255)
    ccTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCt                                                                                                                                            >  1:753586/1‑148 (MQ=255)
    ccTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCt                                                                                                                                            >  1:753579/1‑148 (MQ=255)
      tGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGa                                                                                                                                          <  1:2286775/148‑1 (MQ=255)
      tGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGa                                                                                                                                          >  1:1330142/1‑148 (MQ=255)
        aTTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGAcc                                                                                                                                        >  2:3243768/1‑148 (MQ=255)
          ttACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTg                                                                                                                                      >  2:151141/1‑148 (MQ=255)
           tACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGc                                                                                                                                     <  1:1502995/148‑1 (MQ=255)
           tACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGc                                                                                                                                     >  1:1621229/1‑148 (MQ=255)
             cGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCg                                                                                                                                   <  2:3522465/148‑1 (MQ=255)
              ggCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGa                                                                                                                                  >  2:2314801/1‑148 (MQ=255)
               gCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAg                                                                                                                                 >  1:2252920/1‑148 (MQ=255)
                ccTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGt                                                                                                                                <  1:2883704/148‑1 (MQ=255)
                ccTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGt                                                                                                                                <  1:3015605/148‑1 (MQ=255)
                ccTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGt                                                                                                                                <  1:3015615/148‑1 (MQ=255)
                 cTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTc                                                                                                                               <  1:122483/148‑1 (MQ=255)
                            gACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCaaaaa                                                                                                                    <  2:580825/148‑1 (MQ=255)
                            gACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCaaaaa                                                                                                                    <  2:580869/148‑1 (MQ=255)
                             aCGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAAc                                                                                                                   >  1:602723/1‑148 (MQ=255)
                             aCGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAAc                                                                                                                   >  1:344370/1‑148 (MQ=255)
                              cGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACt                                                                                                                  <  1:873699/148‑1 (MQ=255)
                                ccGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCg                                                                                                                <  1:1692170/148‑1 (MQ=255)
                                ccGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCg                                                                                                                >  1:1571059/1‑148 (MQ=255)
                                   gagaGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCGTCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGata                                                                                                             >  1:2094640/1‑148 (MQ=255)
                                   gagaGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGata                                                                                                             <  2:1680256/148‑1 (MQ=255)
                                    agagTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGatat                                                                                                            >  2:1044665/1‑148 (MQ=255)
                                    agagTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGatat                                                                                                            >  1:1201115/1‑148 (MQ=255)
                                    agagTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGatat                                                                                                            >  1:1200441/1‑148 (MQ=255)
                                     gagTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATc                                                                                                           <  1:1584728/148‑1 (MQ=255)
                                     gagTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATc                                                                                                           <  2:133267/148‑1 (MQ=255)
                                      agTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATcc                                                                                                          >  2:1149304/1‑148 (MQ=255)
                                        tttCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCCGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTg                                                                                                        >  2:2624103/1‑148 (MQ=255)
                                        tttCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTg                                                                                                        >  2:2624097/1‑148 (MQ=255)
                                               tttACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGaaa                                                                                                 <  1:3515221/148‑1 (MQ=255)
                                                ttACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAAc                                                                                                >  1:2411668/1‑148 (MQ=255)
                                                 tACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAGATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAAcc                                                                                               >  2:2268129/1‑148 (MQ=255)
                                                   cccTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCAt                                                                                             <  2:2608195/148‑1 (MQ=255)
                                                           cGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTcc                                                                                     <  2:256230/148‑1 (MQ=255)
                                                           cGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTcc                                                                                     >  2:2381267/1‑148 (MQ=255)
                                                              gtgGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCGTGa                                                                                  <  2:112294/148‑1 (MQ=255)
                                                              gtgGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCGTGa                                                                                  <  2:1035175/148‑1 (MQ=255)
                                                                  cGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCa                                                                              <  2:2911625/148‑1 (MQ=255)
                                                                   ggAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCaa                                                                             <  2:1631825/148‑1 (MQ=255)
                                                                    gAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAAt                                                                            >  1:2895416/1‑148 (MQ=255)
                                                                    gAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAAt                                                                            >  2:2666393/1‑148 (MQ=255)
                                                                       cTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCgg                                                                         >  1:3516796/1‑148 (MQ=255)
                                                                        tGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCggg                                                                        <  2:3190647/148‑1 (MQ=255)
                                                                        tGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCggg                                                                        <  2:605415/148‑1 (MQ=255)
                                                                         gAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGc                                                                       <  2:316013/148‑1 (MQ=255)
                                                                         gAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGc                                                                       <  2:287320/148‑1 (MQ=255)
                                                                         gAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGc                                                                       <  2:692803/148‑1 (MQ=255)
                                                                          aaCCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCt                                                                      <  2:3149176/148‑1 (MQ=255)
                                                                          aaCCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCt                                                                      >  2:682274/1‑148 (MQ=255)
                                                                           aCCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTa                                                                     >  1:1046148/1‑148 (MQ=255)
                                                                            ccccGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTAt                                                                    >  1:3335454/1‑148 (MQ=255)
                                                                                  tCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGttt                                                              >  2:2085512/1‑148 (MQ=255)
                                                                                   cGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTa                                                             <  2:1200441/148‑1 (MQ=255)
                                                                                   cGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTa                                                             <  2:1201115/148‑1 (MQ=255)
                                                                                   cGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTa                                                             >  2:3319003/1‑148 (MQ=255)
                                                                                    gTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTAt                                                            <  2:893531/148‑1 (MQ=255)
                                                                                     tCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATc                                                           >  1:3349692/1‑148 (MQ=255)
                                                                                     tCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATc                                                           >  1:1713396/1‑148 (MQ=255)
                                                                                        gAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGt                                                        <  2:993914/148‑1 (MQ=255)
                                                                                           tgtCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATg                                                     <  1:1179737/148‑1 (MQ=255)
                                                                                           tgtCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATg                                                     <  2:1262876/148‑1 (MQ=255)
                                                                                               tttAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATc                                                 <  1:2283683/148‑1 (MQ=255)
                                                                                               tttAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATc                                                 >  1:635337/1‑148 (MQ=255)
                                                                                                ttAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCa                                                <  2:2959723/148‑1 (MQ=255)
                                                                                                ttAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCa                                                <  2:2960567/148‑1 (MQ=255)
                                                                                                    cTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACttt                                            >  1:1293492/1‑148 (MQ=255)
                                                                                                      aCGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGTCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGc                                          >  2:2117078/1‑148 (MQ=255)
                                                                                                      aCGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGc                                          >  2:1345828/1‑148 (MQ=255)
                                                                                                      aCGCGCATCTGCCGACCATTTTGGCGGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGc                                          <  1:2139512/148‑1 (MQ=255)
                                                                                                        gcgcATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGccc                                        <  2:3533717/148‑1 (MQ=255)
                                                                                                            aTCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTcc                                    <  1:1303095/148‑1 (MQ=255)
                                                                                                            aTCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTcc                                    >  1:769987/1‑148 (MQ=255)
                                                                                                            aTCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTcc                                    <  1:3072355/148‑1 (MQ=255)
                                                                                                            aTCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTcc                                    <  1:3073651/148‑1 (MQ=255)
                                                                                                              cTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCgg                                  >  2:2047251/1‑148 (MQ=255)
                                                                                                                gCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGAt                                <  2:1330142/148‑1 (MQ=255)
                                                                                                                  cGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGa                              >  1:3418900/1‑148 (MQ=255)
                                                                                                                    aCCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACg                            <  2:2387203/148‑1 (MQ=255)
                                                                                                                     ccATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGa                           >  2:72516/1‑148 (MQ=255)
                                                                                                                         ttttGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTg                       >  2:1053430/1‑148 (MQ=255)
                                                                                                                          tttGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTgg                      >  1:321699/1‑148 (MQ=255)
                                                                                                                               tgctCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAACAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTg                 <  2:1825070/148‑1 (MQ=255)
                                                                                                                                gctCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTgg                >  1:569120/1‑148 (MQ=255)
                                                                                                                                    aGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCa            >  1:2717528/1‑148 (MQ=255)
                                                                                                                                        cAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGt        <  2:2744787/148‑1 (MQ=255)
                                                                                                                                        cAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGt        >  1:1854486/1‑148 (MQ=255)
                                                                                                                                        cAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGt        >  1:1653538/1‑148 (MQ=255)
                                                                                                                                          aaaTCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGa      >  2:461119/1‑148 (MQ=255)
                                                                                                                                          aaaTCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGa      <  2:302343/148‑1 (MQ=255)
                                                                                                                                          aaaTCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGa      <  2:302470/148‑1 (MQ=255)
                                                                                                                                           aaTCAACGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGaa     >  1:1492992/1‑148 (MQ=255)
                                                                                                                                            aTCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAg    <  2:807241/148‑1 (MQ=255)
                                                                                                                                              cAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGc  >  2:416578/1‑148 (MQ=255)
                                                                                                                                              cAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGc  >  1:2621782/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TCGACCTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGC  >  CP009273/4046027‑4046316

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: