Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,564,984 C→A 33.8% E627* (GAA→TAA)  glgB ← 1,4‑alpha‑glucan branching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,564,9840CA33.8% 374.0 / 191.0 215E627* (GAA→TAA) glgB1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base C (75/67);  new base A (44/29);  total (119/96)
Fisher's exact test for biased strand distribution p-value = 3.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 219 total aligned reads displayed.

ACGCCATTTGCCCGGCTGGTTTATGCCGAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACT  >  CP009273/3564835‑3565133
                                                                                                                                                     |                                                                                                                                                     
aCGCCATTTGCCCGGCTGGTTTATGCCGAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTc                                                                                                                                                       >  1:1975021/1‑150 (MQ=255)
  gCCATTTGCCCGGCTGGTTTATGCCGAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCaa                                                                                                                                                     >  1:2582280/1‑150 (MQ=255)
   ccATTTGCCCGGCTGGTTTATGCCGAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCaaa                                                                                                                                                    <  2:1780500/150‑1 (MQ=255)
    cATTTGCCCGGCTGGTTTATGCCGAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAg                                                                                                                                                   <  2:444580/150‑1 (MQ=255)
    cATTTGCCCGGCTGGTTTATGCCGAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAg                                                                                                                                                   >  2:2066659/1‑150 (MQ=255)
       ttGCCCGGCTGGTTTATGCCGAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTAAAAGccc                                                                                                                                                >  2:3622821/1‑149 (MQ=255)
        tGCCCGGCTGGTTTATGCCGAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTAAAAGCCGt                                                                                                                                               >  2:889000/1‑150 (MQ=255)
         gCCCGGCTGGTTTATGCCGAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTAAAAGCCGTa                                                                                                                                              >  2:899535/1‑150 (MQ=255)
                tGGTTTATGCCGAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCa                                                                                                                                       <  1:2802813/150‑1 (MQ=255)
                 ggTTTATGCCGAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCaa                                                                                                                                      >  1:4152044/1‑150 (MQ=255)
                 ggTTTATGCCGAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTAAAAGCCGTACGGGTCaa                                                                                                                                      <  2:938749/150‑1 (MQ=255)
                   tttATGCCGAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCaaaa                                                                                                                                    <  2:3240038/150‑1 (MQ=255)
                     tATGCCGAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATc                                                                                                                                  <  1:1084615/150‑1 (MQ=255)
                      aTGCCGAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTAAAAGCCGTACGGGTCAAAATcc                                                                                                                                 >  1:650821/1‑150 (MQ=255)
                        gCCGAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAg                                                                                                                               <  2:1306153/150‑1 (MQ=255)
                          cGAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGtt                                                                                                                             <  2:3054442/150‑1 (MQ=255)
                            aaGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCa                                                                                                                           >  1:4273656/1‑150 (MQ=255)
                             aGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCAt                                                                                                                          >  2:1460375/1‑150 (MQ=255)
                              gCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATg                                                                                                                         >  1:4386644/1‑150 (MQ=255)
                               cGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGc                                                                                                                        >  2:2504320/1‑150 (MQ=255)
                                gATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCa                                                                                                                       >  1:3003517/1‑150 (MQ=255)
                                 aTAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCAt                                                                                                                      >  2:1385018/1‑150 (MQ=255)
                                    aTCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGc                                                                                                                   <  1:2080785/150‑1 (MQ=255)
                                    aTCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGc                                                                                                                   >  1:175422/1‑150 (MQ=255)
                                      cATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCtt                                                                                                                 >  2:1003797/1‑150 (MQ=255)
                                      cATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCtt                                                                                                                 <  1:1909281/150‑1 (MQ=255)
                                      cATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCtt                                                                                                                 <  1:3892527/150‑1 (MQ=255)
                                        tGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTa                                                                                                               <  2:362877/150‑1 (MQ=255)
                                        tGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTa                                                                                                               <  1:837674/150‑1 (MQ=255)
                                         gACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTAt                                                                                                              <  1:1225097/150‑1 (MQ=255)
                                           cgcgGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCGTTAtgg                                                                                                            >  2:4319190/1‑150 (MQ=255)
                                            gcgGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTAtggt                                                                                                           <  2:2451308/150‑1 (MQ=255)
                                               gTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGc                                                                                                        <  2:3008827/150‑1 (MQ=255)
                                               gTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGc                                                                                                        <  1:2963808/150‑1 (MQ=255)
                                                tACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCg                                                                                                       <  1:2142576/150‑1 (MQ=255)
                                                 aCCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCgg                                                                                                      <  1:2778251/150‑1 (MQ=255)
                                                 aCCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCgg                                                                                                      >  2:3441122/1‑150 (MQ=255)
                                                   cGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGACAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTa                                                                                                    >  1:507605/1‑150 (MQ=255)
                                                    ggCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAg                                                                                                   <  1:4206310/150‑1 (MQ=255)
                                                       gTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTg                                                                                                >  2:3191396/1‑150 (MQ=255)
                                                        tAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGa                                                                                               >  2:3000630/1‑150 (MQ=255)
                                                         aaaGTTTCTGGCAACGATGATTTCGTTACCCTCTTTATCGAGAAGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTAAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAg                                                                                              <  1:1008809/150‑1 (MQ=255)
                                                         aaaGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAg                                                                                              <  1:397330/150‑1 (MQ=255)
                                                           aGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGt                                                                                            <  2:3109987/150‑1 (MQ=255)
                                                            gTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGtt                                                                                           <  1:1091081/150‑1 (MQ=255)
                                                                cTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGa                                                                                       <  2:3551479/150‑1 (MQ=255)
                                                                cTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGATCTTTGTCATCCACCACCAGCCATAAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGa                                                                                       >  2:4112857/1‑150 (MQ=255)
                                                                 tGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGAt                                                                                      >  2:1516296/1‑150 (MQ=255)
                                                                      aCGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCCTTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCgcgc                                                                                 >  2:3134298/1‑150 (MQ=255)
                                                                        gatgatTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCAc                                                                               >  2:1457259/1‑150 (MQ=255)
                                                                           gatTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAg                                                                            <  2:1091973/150‑1 (MQ=255)
                                                                           gatTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAg                                                                            >  2:1831151/1‑150 (MQ=255)
                                                                              ttCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACg                                                                         <  1:246672/150‑1 (MQ=255)
                                                                                 gTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTg                                                                      <  2:2694099/150‑1 (MQ=255)
                                                                                 gTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTg                                                                      <  2:3694620/150‑1 (MQ=255)
                                                                                 gTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTg                                                                      >  2:511737/1‑150 (MQ=255)
                                                                                     cccTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGaca                                                                  <  2:1615745/150‑1 (MQ=255)
                                                                                      ccTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGATTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGacac                                                                 <  2:1901748/150‑1 (MQ=255)
                                                                                          tttATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCgtg                                                             >  1:3718908/1‑150 (MQ=255)
                                                                                           atATCGCGACGCACAAAGATCAGCACCGAGCGTACTTTGTCATCAACCAACAGCCATACAAAGCCGTACGGGTCAAAATCCAGTTCAAGCATTTCTTTATGGTGGCTGAAGGTGAGGTTCAGATCGCGCACCAGACGCAGGACACCgagg                                                            >  1:3296614/2‑150 (MQ=255)
                                                                                               cgcgACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGGGGTGcc                                                        >  2:4178973/1‑150 (MQ=255)
                                                                                                  gACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGt                                                     >  2:609390/1‑150 (MQ=255)
                                                                                                  gACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGt                                                     >  2:1211307/1‑150 (MQ=255)
                                                                                                     gCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTAt                                                  >  2:3391283/1‑150 (MQ=255)
                                                                                                     gCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGAGCACCAGACGCTGGACACCGTGGTGCCAGTTAt                                                  >  2:3647456/1‑150 (MQ=255)
                                                                                                        caAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATcgc                                               >  2:1117913/1‑150 (MQ=255)
                                                                                                         aaaGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATcgcc                                              <  2:2823561/150‑1 (MQ=255)
                                                                                                         aaaGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATcgcc                                              >  2:2766462/1‑150 (MQ=255)
                                                                                                         aaaGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATcgcc                                              <  2:2958663/150‑1 (MQ=255)
                                                                                                             aTCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCt                                          <  2:4415225/150‑1 (MQ=255)
                                                                                                               cAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTc                                        <  2:3787190/150‑1 (MQ=255)
                                                                                                               cAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTc                                        <  1:1264731/150‑1 (MQ=255)
                                                                                                               cAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTc                                        >  2:3644711/1‑150 (MQ=255)
                                                                                                                aGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTcc                                       <  1:1376387/150‑1 (MQ=255)
                                                                                                                  cACCGAGCGTTCTTTGTCATCCACCACCAGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCaa                                     <  2:4040481/150‑1 (MQ=255)
                                                                                                                   aCCGAGCGTTCTTTGTCATCCACCACCAGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAAc                                    >  1:1831483/1‑150 (MQ=255)
                                                                                                                   aCCGAGCGTTCTTTGTCATCCACCACCAGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAAc                                    >  1:1831534/1‑150 (MQ=255)
                                                                                                                       aGCGTTCTTTGTCATCCACCACCAGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGAt                                <  1:1196876/150‑1 (MQ=255)
                                                                                                                        gCGTTCTTTGTCATCCACCACCAGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATg                               <  2:1494128/150‑1 (MQ=255)
                                                                                                                          gTTCTTTGTCATCCACCACCAGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGcc                             >  1:2678823/1‑150 (MQ=255)
                                                                                                                            tCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAg                           <  2:1651195/150‑1 (MQ=255)
                                                                                                                            tCTTTGTCATCCACCACCAGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAg                           >  2:4560671/1‑150 (MQ=255)
                                                                                                                              tttGTCATCCACCCCCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTc                         <  2:2577406/150‑1 (MQ=255)
                                                                                                                               ttGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCg                        >  2:133700/1‑150 (MQ=255)
                                                                                                                                 gTCATCCACCACCAGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAg                      <  2:3904374/150‑1 (MQ=255)
                                                                                                                                    aTCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCt                   <  2:2164637/150‑1 (MQ=255)
                                                                                                                                     tCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTg                  >  1:146552/1‑150 (MQ=255)
                                                                                                                                     tCCACCACCAGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTg                  <  1:2323147/150‑1 (MQ=255)
                                                                                                                                          caccaGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTc             >  1:3191409/1‑150 (MQ=255)
                                                                                                                                           accaGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCa            <  2:2609072/150‑1 (MQ=255)
                                                                                                                                            ccaGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCAt           >  2:1947878/1‑150 (MQ=255)
                                                                                                                                             caGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATg          >  2:479366/1‑150 (MQ=255)
                                                                                                                                             caGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATg          >  2:3716622/1‑150 (MQ=255)
                                                                                                                                               gCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGt        >  1:4076250/1‑150 (MQ=255)
                                                                                                                                                 cATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTc      >  1:3133102/1‑150 (MQ=255)
                                                                                                                                                 cATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTc      >  1:340122/1‑150 (MQ=255)
                                                                                                                                                  aTTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTcc     >  2:3408485/1‑150 (MQ=255)
                                                                                                                                                   ttCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCa    <  1:3616191/150‑1 (MQ=255)
                                                                                                                                                    tCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCAc   <  2:2411068/150‑1 (MQ=255)
                                                                                                                                                     cAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACt  >  2:820878/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACGCCATTTGCCCGGCTGGTTTATGCCGAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACT  >  CP009273/3564835‑3565133

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: