Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,208,581 C→A 14.7% E84D (GAG→GAT)  mug ← G/U mismatch‑specific DNA glycosylase; xanthine DNA glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,208,5810CA14.7% 872.1 / 85.4 259E84D (GAG→GAT) mugG/U mismatch‑specific DNA glycosylase; xanthine DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base C (127/94);  new base A (20/18);  total (147/112)
Fisher's exact test for biased strand distribution p-value = 5.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 261 total aligned reads displayed.

TGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGATATATCACCTTCCAGAAGCG  >  CP009273/3208432‑3208727
                                                                                                                                                     |                                                                                                                                                  
  cGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATc                                                                                                                                                  <  2:1387217/150‑1 (MQ=255)
    tCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCt                                                                                                                                                >  2:187578/1‑150 (MQ=255)
     tgAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTg                                                                                                                                               >  1:1576170/2‑150 (MQ=255)
     cGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTg                                                                                                                                               >  1:4216760/1‑150 (MQ=255)
       agccAATGGTGAGCGTTTGTTTCCCCCACTGTGCGCCGCGCTGGCTGAATCCCTGTTCATAGGCTTGTTTGCCCAGAATCGCCCACGCCTGCGGCTGAGAAGCTTCAATTTTTTCTATCAGCTTTCGCCCGCCTGCGTGTAGCTCCTGCt                                                                                                                                             >  1:4246873/3‑150 (MQ=255)
        aCCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCtt                                                                                                                                            >  2:3502469/1‑150 (MQ=255)
        aCCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCtt                                                                                                                                            >  1:3190540/1‑150 (MQ=255)
         ccAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCttt                                                                                                                                           <  2:3006090/150‑1 (MQ=255)
           aaTGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGa                                                                                                                                         >  1:4099117/1‑150 (MQ=255)
            aTGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGaa                                                                                                                                        >  1:1211564/1‑150 (MQ=255)
            aTGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGaa                                                                                                                                        <  2:3853801/150‑1 (MQ=255)
             tGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGaaa                                                                                                                                       <  1:3870921/150‑1 (MQ=255)
               gTGAGCGTTTGTTTCACCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACt                                                                                                                                     <  1:1960886/150‑1 (MQ=255)
                 gAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTc                                                                                                                                   >  1:1764117/1‑150 (MQ=255)
                  aGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCa                                                                                                                                  <  1:2739206/150‑1 (MQ=255)
                    cGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCAtt                                                                                                                                >  1:3485270/1‑150 (MQ=255)
                    cGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCAtt                                                                                                                                <  2:1666345/150‑1 (MQ=255)
                    cGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCAtt                                                                                                                                <  2:3182719/150‑1 (MQ=255)
                     gtttgtttCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTg                                                                                                                               >  1:3348717/1‑150 (MQ=255)
                     gattgtttCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTg                                                                                                                               <  1:2480567/148‑1 (MQ=255)
                       ttgtttCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGc                                                                                                                             <  2:1253701/150‑1 (MQ=255)
                        tgtttCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGGTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCt                                                                                                                            >  1:4370153/1‑150 (MQ=255)
                          tttCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTg                                                                                                                          >  2:376878/1‑150 (MQ=255)
                           ttCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGc                                                                                                                         >  1:2367400/1‑150 (MQ=255)
                               cccACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCg                                                                                                                     <  2:280112/150‑1 (MQ=255)
                                 cACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGtt                                                                                                                   <  1:4027545/150‑1 (MQ=255)
                                  aCTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTg                                                                                                                  >  1:1107884/1‑150 (MQ=255)
                                   cTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTgg                                                                                                                 <  2:2962312/150‑1 (MQ=255)
                                   cTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTgg                                                                                                                 >  2:1143689/1‑150 (MQ=255)
                                   cTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTgg                                                                                                                 >  2:486617/1‑150 (MQ=255)
                                    tgtgCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGa                                                                                                                >  1:1614466/1‑150 (MQ=255)
                                      tgCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACg                                                                                                              >  2:1773723/1‑150 (MQ=255)
                                       gCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACgg                                                                                                             >  2:247388/1‑150 (MQ=255)
                                         aCCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTc                                                                                                           >  1:3954435/1‑150 (MQ=255)
                                           cgcgCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTa                                                                                                         <  2:3080607/150‑1 (MQ=255)
                                               ctggctgAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAg                                                                                                     >  2:965774/1‑150 (MQ=255)
                                                tggctgAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGt                                                                                                    >  1:938787/1‑150 (MQ=255)
                                                    tgAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTgg                                                                                                >  1:1025501/1‑150 (MQ=255)
                                                       aTCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGa                                                                                             <  2:4062367/150‑1 (MQ=255)
                                                       aTCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGa                                                                                             <  2:4114962/150‑1 (MQ=255)
                                                        tCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGAc                                                                                            <  1:4382928/150‑1 (MQ=255)
                                                           cTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGcc                                                                                         <  1:2625522/150‑1 (MQ=255)
                                                            tGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCca                                                                                        >  1:2881308/1‑150 (MQ=255)
                                                              ttCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCcaca                                                                                      >  2:1037405/1‑150 (MQ=255)
                                                              ttCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCcaca                                                                                      <  2:1644555/150‑1 (MQ=255)
                                                              ttCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCcaca                                                                                      <  2:1851341/150‑1 (MQ=255)
                                                                cATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAAc                                                                                    >  2:3983688/1‑150 (MQ=255)
                                                                 atatGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACg                                                                                   >  1:447548/1‑150 (MQ=255)
                                                                 atatGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACg                                                                                   >  1:2119427/1‑150 (MQ=255)
                                                                  tatGCTTGATTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTATCAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGa                                                                                  >  2:2617132/1‑150 (MQ=255)
                                                                   atGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGAt                                                                                 >  1:3212754/1‑150 (MQ=255)
                                                                   atGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGAt                                                                                 >  1:1329214/1‑150 (MQ=255)
                                                                         gTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCca                                                                           >  1:4079092/1‑150 (MQ=255)
                                                                           ttGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCcagc                                                                         >  1:2165561/1‑150 (MQ=255)
                                                                              cccAGACTCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGa                                                                      <  1:3531308/150‑1 (MQ=255)
                                                                                    aTCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTTGCTCCTGATTTGAAACTTTTTTGGCTTGCACAGTTGGACGGTATTCCAGTTTGGTGACGCCTCAACGATAATCCAGCAGATGCtgt                                                                <  2:600171/150‑1 (MQ=255)
                                                                                     tCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCtgtg                                                               >  1:4071578/1‑150 (MQ=255)
                                                                                        ccAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGcct                                                            >  1:3786636/1‑150 (MQ=255)
                                                                                         cAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGcctc                                                           >  2:3345526/1‑150 (MQ=255)
                                                                                         cAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGcctc                                                           >  2:2420931/1‑150 (MQ=255)
                                                                                          aaCGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGcctcc                                                          >  2:3566268/1‑150 (MQ=255)
                                                                                          aaCGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGcctcc                                                          <  1:1051459/150‑1 (MQ=255)
                                                                                                  cGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCtt                                                  >  1:4204631/1‑150 (MQ=255)
                                                                                                   ggCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTc                                                 <  2:1016184/150‑1 (MQ=255)
                                                                                                       gATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATTCAGCAGATGCTGTGCCTCCTGCGGCTTCCACt                                             >  2:3442216/1‑150 (MQ=255)
                                                                                                       gATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACt                                             <  2:1211723/150‑1 (MQ=255)
                                                                                                       gATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACt                                             >  2:2015501/1‑150 (MQ=255)
                                                                                                       gATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACt                                             <  2:2021391/150‑1 (MQ=255)
                                                                                                        aTAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTg                                            >  1:2496895/1‑150 (MQ=255)
                                                                                                          aaTCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGAc                                          >  1:1902137/1‑150 (MQ=255)
                                                                                                          aaTCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGAc                                          >  2:1715514/1‑150 (MQ=255)
                                                                                                            tCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGAcgg                                        >  2:3227393/1‑150 (MQ=255)
                                                                                                                 aaTTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGAcggtcggt                                   <  2:4490933/150‑1 (MQ=255)
                                                                                                                      tttCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAAcc                              >  1:2918679/1‑150 (MQ=255)
                                                                                                                      tttCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAAcc                              >  1:3440838/1‑150 (MQ=255)
                                                                                                                      attCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAAcc                              >  1:3845534/2‑150 (MQ=255)
                                                                                                                      attCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAAcc                              >  1:3110566/2‑150 (MQ=255)
                                                                                                                       ttCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAAccc                             <  1:820333/150‑1 (MQ=255)
                                                                                                                       ttCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAAccc                             <  2:2205593/150‑1 (MQ=255)
                                                                                                                       ttCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAAccc                             <  2:776981/150‑1 (MQ=255)
                                                                                                                        tCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCg                            <  2:3190540/150‑1 (MQ=255)
                                                                                                                        tCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCg                            >  1:4440416/1‑150 (MQ=255)
                                                                                                                        tCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCg                            >  1:2145998/1‑150 (MQ=255)
                                                                                                                        tCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCGGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCg                            >  2:3458120/1‑150 (MQ=255)
                                                                                                                          aaTCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGc                          >  1:3157523/1‑150 (MQ=255)
                                                                                                                             cAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTg                       >  2:1537932/1‑150 (MQ=255)
                                                                                                                             cAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTg                       <  2:1479940/150‑1 (MQ=255)
                                                                                                                             cAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTg                       >  2:4274780/1‑150 (MQ=255)
                                                                                                                             cAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTg                       >  1:2331347/1‑150 (MQ=255)
                                                                                                                               gCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGat                     <  2:576537/150‑1 (MQ=255)
                                                                                                                                cTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGata                    >  2:2189480/1‑150 (MQ=255)
                                                                                                                                 ttACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGatat                   >  2:2483660/1‑150 (MQ=255)
                                                                                                                                 ttACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGatat                   >  1:1473612/1‑150 (MQ=255)
                                                                                                                                    cgcccgccTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGATATATc                <  1:376878/150‑1 (MQ=255)
                                                                                                                                     gcccgccTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGATATATCa               >  2:565246/1‑150 (MQ=255)
                                                                                                                                       ccgccTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGATATATCAcc             <  1:498645/150‑1 (MQ=255)
                                                                                                                                          ccTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGATATATCACCTTc          <  2:766816/150‑1 (MQ=255)
                                                                                                                                              cGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGATATATCACCTTCCAGa      >  1:541513/1‑150 (MQ=255)
                                                                                                                                              cGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGATATATCACCTTCCAGa      >  1:1067787/1‑150 (MQ=255)
                                                                                                                                                 gtAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGATATATCACCTTCCAGAAGc   >  1:4143562/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGATATATCACCTTCCAGAAGCG  >  CP009273/3208432‑3208727

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: ‑