Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,295,911 C→A 100% Q22K (CAG→AAG)  yhbV → U32 peptidase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,295,9110CA98.2% 535.2 / 1.1 169Q22K (CAG→AAG) yhbVU32 peptidase family protein
Reads supporting (aligned to +/- strand):  ref base C (3/0);  new base A (94/72);  total (97/72)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 173 total aligned reads displayed.

GGATGGAGACGCTCGGGTCGATGTCCGAAGGCACGCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATCAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAATCGCTCGCCGGGAGTGGTAAGCAGATTGTGCTCTCCACGCTGG  >  CP009273/3295766‑3296055
                                                                                                                                                 |                                                                                                                                                
ggatggaGACGCTCGGGTCGATGTCCGAAGGCACGCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATaag                                                                                                                                                <  2:713710/148‑1 (MQ=255)
ggatggaGACGCTCGGGTCGATGTCCGAAGGCACGCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATaag                                                                                                                                                >  2:2193658/1‑148 (MQ=255)
  atggaGACGCTCGGGTCGATGTCCGAAGGCACGCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATaagca                                                                                                                                              >  2:1648101/1‑148 (MQ=255)
    ggaGACGCTCGGGTCGATGTCCGAAGGCACGCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAgg                                                                                                                                            <  1:3252342/148‑1 (MQ=255)
    ggaGACGCTCGGGTCGATGTCCGAAGGCACGCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAgg                                                                                                                                            >  2:2919495/1‑148 (MQ=255)
    ggaGACGCTCGGGTCGATGTCCGAAGGCACGCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAgg                                                                                                                                            >  2:2919508/1‑148 (MQ=255)
    ggaGACGCTCGGGTCGATGTCCGAAGGCACGCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAgg                                                                                                                                            <  2:3423345/148‑1 (MQ=255)
      agaCGCTCGGGTCGATGTCCGAAGGCACGCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGgcc                                                                                                                                          >  2:2853094/1‑148 (MQ=255)
        aCGCTCGGGTCGATGTCCGAAGGCACGCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGgccgc                                                                                                                                        >  2:455760/1‑148 (MQ=255)
        aCGCTCGGGTCGATGTCCGAAGGCACGCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGgccgc                                                                                                                                        >  2:1345241/1‑148 (MQ=255)
        aCGCTCGGGTCGATGTCCGAAGGCACGCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGgccgc                                                                                                                                        >  1:1107619/1‑148 (MQ=255)
           cTCGGGTCGATGTCCGAAGGCACGCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGccac                                                                                                                                     >  1:2544432/1‑148 (MQ=255)
           cTCGGGTCGATGTCCGAAGGCACGCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGccac                                                                                                                                     >  1:2518553/1‑148 (MQ=255)
              gggTCGATGTCCGAAGGCACGCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACcag                                                                                                                                  <  2:1389090/148‑1 (MQ=255)
               ggTCGATGTCCGAAGGCACGCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACcagc                                                                                                                                 >  1:414914/1‑148 (MQ=255)
                    aTGTCCGAAGGCACGCAAACCACTCTTGGCGCATATCGCCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAgcgc                                                                                                                            <  2:2648327/148‑1 (MQ=255)
                    aTGTCCGAAGGCACGCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAgcgc                                                                                                                            >  1:640148/1‑148 (MQ=255)
                    aTGTCCGAAGGCACGCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAgcgc                                                                                                                            <  1:2515434/148‑1 (MQ=255)
                      gTCCGAAGGCACGCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCg                                                                                                                          >  1:1685220/1‑148 (MQ=255)
                      gTCCGAAGGCACGCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCg                                                                                                                          >  2:3015881/1‑148 (MQ=255)
                          gAAGGCACGCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGt                                                                                                                      <  1:560543/148‑1 (MQ=255)
                          gAAGGCACGCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGt                                                                                                                      <  1:538613/148‑1 (MQ=255)
                            aGGCACGCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGa                                                                                                                    >  2:353400/1‑148 (MQ=255)
                            aGGCACGCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGa                                                                                                                    >  1:131363/1‑148 (MQ=255)
                             ggCACGCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGAt                                                                                                                   >  1:504674/1‑148 (MQ=255)
                                 cGCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTAt                                                                                                               <  2:2726883/148‑1 (MQ=255)
                                 cGCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTAt                                                                                                               <  2:2598706/148‑1 (MQ=255)
                                  gCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATc                                                                                                              <  1:857990/148‑1 (MQ=255)
                                      aCCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTgg                                                                                                          >  2:1376952/1‑148 (MQ=255)
                                       ccACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGt                                                                                                         <  1:2146350/148‑1 (MQ=255)
                                             ttGGCGCATATCACCGTAAATGGCAGTTATTAAAGCAATGTAATATTCCTTTGGGGCTTTGCTGTGGTACTGGCCAAAAGAGACTCTGGAAGTATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCg                                                                                                   <  1:2087698/148‑1 (MQ=255)
                                                gcgcATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCGGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTa                                                                                                >  1:2658068/1‑148 (MQ=255)
                                                gcgcATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTa                                                                                                >  1:2792818/1‑148 (MQ=255)
                                                    atatCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATgca                                                                                            <  2:2502800/148‑1 (MQ=255)
                                                     tatCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATgcag                                                                                           >  1:2376271/1‑148 (MQ=255)
                                                      atCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATgcagc                                                                                          <  1:2894777/148‑1 (MQ=255)
                                                       tCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATgcagca                                                                                         >  1:935887/1‑148 (MQ=255)
                                                             tAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTc                                                                                   >  2:2405802/1‑148 (MQ=255)
                                                             tAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTc                                                                                   >  2:855419/1‑148 (MQ=255)
                                                              aaaTGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCg                                                                                  <  2:2524832/148‑1 (MQ=255)
                                                                aTGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACTTGATTTATCTTGGTGAAGCTGTATGCAGCAAGCGTCtgg                                                                                >  1:3141416/1‑148 (MQ=255)
                                                                 tGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGc                                                                               >  2:1288680/1‑148 (MQ=255)
                                                                   gCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCaa                                                                             <  2:1905071/148‑1 (MQ=255)
                                                                        gagaAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCaaa                                                                        >  2:1720474/1‑148 (MQ=255)
                                                                         agaAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAg                                                                       <  1:732618/148‑1 (MQ=255)
                                                                          gaACAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGt                                                                      <  2:2055605/148‑1 (MQ=255)
                                                                          gaAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATCAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGt                                                                      >  2:3095737/1‑148 (MQ=255)
                                                                              aGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGc                                                                  >  1:379548/1‑148 (MQ=255)
                                                                              aGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGc                                                                  >  1:1231710/1‑148 (MQ=255)
                                                                               gCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCg                                                                 >  2:610659/1‑148 (MQ=255)
                                                                               gCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTGGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCg                                                                 >  2:2002178/1‑148 (MQ=255)
                                                                                  aTGAAATATTCCTTAGGGCCAGTGCTGTGGTCCTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGAct                                                              >  1:3019745/1‑148 (MQ=255)
                                                                                  aTGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGAct                                                              >  1:3019770/1‑148 (MQ=255)
                                                                                    gAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGActgg                                                            >  1:1551326/1‑148 (MQ=255)
                                                                                      aaTATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGActggct                                                          <  1:243318/148‑1 (MQ=255)
                                                                                         atTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGga                                                       <  2:1354653/148‑1 (MQ=255)
                                                                                         atTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGga                                                       <  2:2885639/148‑1 (MQ=255)
                                                                                          ttCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGTAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGgag                                                      >  2:3349301/1‑148 (MQ=255)
                                                                                           tCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGgaga                                                     >  1:467979/1‑148 (MQ=255)
                                                                                           tCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGgaga                                                     >  1:2742249/1‑148 (MQ=255)
                                                                                             cTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATg                                                   <  1:1617096/148‑1 (MQ=255)
                                                                                              ttAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATgg                                                  >  1:173684/1‑148 (MQ=255)
                                                                                                 gggCTAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCaa                                               >  2:3045938/1‑148 (MQ=255)
                                                                                                     cAGTGCTGTGGTACTTGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGAGTTATCGTGGTGAGGGGGTATGCAGCACGCGTCGGGCACCCAAAGTTGGCGACTGGCTGGAGATGGCAAAAtc                                           >  2:297322/1‑148 (MQ=255)
                                                                                                     cAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAAtc                                           >  2:981686/1‑148 (MQ=255)
                                                                                                      aGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAAtcg                                          <  2:2098186/148‑1 (MQ=255)
                                                                                                      aGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAAtcg                                          <  2:793005/148‑1 (MQ=255)
                                                                                                        tGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAAtcgct                                        >  2:2237675/1‑148 (MQ=255)
                                                                                                           tgtgGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAAtcgctcgc                                     <  2:392608/148‑1 (MQ=255)
                                                                                                             tgGTACTGGACAAAAGAGACGCTTGAAGAATTTTATAAGCAGGCCGCCACCAGAAGCGCCGACTTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATTGCAAAATCTCTCGACg                                   >  1:3504991/1‑148 (MQ=255)
                                                                                                                tACTGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAATCGCTCGCCGGGa                                >  1:1847737/1‑148 (MQ=255)
                                                                                                                   tGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAATCGCTCGCCGGGAGTg                             <  1:2245022/148‑1 (MQ=255)
                                                                                                                   tGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAATCGCTCGCCGGGAGTg                             <  1:1345241/148‑1 (MQ=255)
                                                                                                                   tGGCCAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAATCGCTCGCCGGGAGTg                             <  1:455760/148‑1 (MQ=255)
                                                                                                                      ccAAAAGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAATCGCTCGCCGGGAGTGGTa                          <  1:1376891/148‑1 (MQ=255)
                                                                                                                         aaaGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAATCGCTCGCCGGGAGTGGTAAGc                       >  2:1153617/1‑148 (MQ=255)
                                                                                                                          aaGAGACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAATCGCTCGCCGGGAGTGGTAAGCa                      <  2:1278210/148‑1 (MQ=255)
                                                                                                                           agagACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAATCGCTCGCCGGGAGTGGTAAGCAg                     >  1:1923817/1‑148 (MQ=255)
                                                                                                                           agagACGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAATCGCTCGCCGGGAGTGGTAAGCAg                     <  1:1281861/148‑1 (MQ=255)
                                                                                                                             agaCGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAATCGCTCGCCGGGAGTGGTAAGCAGAt                   <  1:1363090/148‑1 (MQ=255)
                                                                                                                             agaCGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAATCGCTCGCCGGGAGTGGTAAGCAGAt                   >  2:1502111/1‑148 (MQ=255)
                                                                                                                              gaCGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAATCGCTCGCCGGGAGTGGTAAGCAGAtt                  <  2:559073/148‑1 (MQ=255)
                                                                                                                               aCGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAATCGCTCGCCGGGAGTGGTAAGCAGATtg                 >  1:224800/1‑148 (MQ=255)
                                                                                                                               aCGCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAATCGCTCGCCGGGAGTGGTAAGCAGATtg                 <  1:1974281/148‑1 (MQ=255)
                                                                                                                                 gCTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCTAAATCGCTCGCCGGGAGTGGTCAGCAGATtgtg               >  2:1807795/1‑148 (MQ=255)
                                                                                                                                  cTGGAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAATCGCTCGCCGGGAGTGGTAAGCAGATTGTGc              <  1:3015881/148‑1 (MQ=255)
                                                                                                                                    ggAAGAATTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAATCGCTCGCCGGGAGTGGTAAGCAGATTGTGctc            <  2:414914/148‑1 (MQ=255)
                                                                                                                                      aagaaTTTTATCAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAATCGCTCGCCGGGAGTGGTAAGCAGATTGTGCtctc          >  1:805509/1‑148 (MQ=255)
                                                                                                                                      aagaaTTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAATCGCTCGCCGGGAGTGGTAAGCAGATTGTGCtctc          >  1:1398341/1‑148 (MQ=255)
                                                                                                                                      aagaaTTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAATCGCTCGCCGGGAGTGGTAAGCAGATTGTGCtctc          <  1:2919508/148‑1 (MQ=255)
                                                                                                                                       agaaTTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAATCGCTCGCCGGGAGTGGTAAGCAGATTGTGCTCTcc         >  1:1263022/1‑148 (MQ=255)
                                                                                                                                       agaaTTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAATCGCTCGCCGGGAGTGGTAAGCAGATTGTGCTCTcc         >  2:2311389/1‑148 (MQ=255)
                                                                                                                                       agaaTTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGAGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAATCGCTCGCCGGGAGTGGTAAGCAGATTGTGCTCTcc         <  2:1743511/148‑1 (MQ=255)
                                                                                                                                        gaaTTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAATCGCTCGCCGGGAGTGGTAAGCAGATTGTGCTCTCCa        <  2:499947/148‑1 (MQ=255)
                                                                                                                                        gaaTTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAATCGCTCGCCGGGAGTGGTAAGCAGATTGTGCTCTCCa        >  1:1051100/1‑148 (MQ=255)
                                                                                                                                        gaaTTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAATCGCTCGCCGGGAGTGGTAAGCAGATTGTGCTCTCCa        >  1:1038437/1‑148 (MQ=255)
                                                                                                                                        gaaTTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACCGGCTGGAGATGGCAAAATCGCTCGCCGGGAGTGGTAAGCAGATTGTGCTCTCCa        >  1:1038489/1‑148 (MQ=255)
                                                                                                                                          aTTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAATCGCTCGCCGGGAGTGGTAAGCAGATTGTGCTCTCCAcg      >  1:1963301/1‑148 (MQ=255)
                                                                                                                                          aTTTTATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAATCGCTCGCCGGGAGTGGTAAGCAGATTGTGCTCTCCAcg      <  1:2877541/148‑1 (MQ=255)
                                                                                                                                           ttttATAAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAATCGCTCGCCGGGAGTGGTAAGCAGATTGTGCTCTCCAcgc     >  1:2908320/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GGATGGAGACGCTCGGGTCGATGTCCGAAGGCACGCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATCAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGTGAAGCGGTATGCAGCAAGCGTCGGGCAACCAAAGTTGGCGACTGGCTGGAGATGGCAAAATCGCTCGCCGGGAGTGGTAAGCAGATTGTGCTCTCCACGCTGG  >  CP009273/3295766‑3296055

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: