Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,732,549 G→T 5.7% D565Y (GAT→TAT)  malS → alpha‑amylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,732,5490GT 5.7% 748.9 / 13.8 199D565Y (GAT→TAT) malSalpha‑amylase
Reads supporting (aligned to +/- strand):  ref base G (89/98);  new base T (3/9);  total (92/107)
Fisher's exact test for biased strand distribution p-value = 1.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 200 total aligned reads displayed.

CACGGCTTCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTACACGTTCGGATA  >  CP009273/3732405‑3732696
                                                                                                                                                |                                                                                                                                                   
  cGGCTTCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATAcc                                                                                                                                                >  2:1149985/1‑148 (MQ=255)
       tCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCt                                                                                                                                           >  2:2230586/1‑148 (MQ=255)
        cgatgcgatgATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTg                                                                                                                                          <  2:1739007/148‑1 (MQ=255)
        cgatgcgatgATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTg                                                                                                                                          <  2:887952/148‑1 (MQ=255)
        cgatgcgatgATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTg                                                                                                                                          >  1:1979035/1‑148 (MQ=255)
          atgcgatgATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGtt                                                                                                                                        >  2:1021772/1‑148 (MQ=255)
             cgatgATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCg                                                                                                                                     <  2:930302/148‑1 (MQ=255)
                tgatCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGa                                                                                                                                  <  2:1837162/148‑1 (MQ=255)
                tgatCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCAGTTCCGTGa                                                                                                                                  >  2:520307/1‑148 (MQ=255)
                 gatCAATTTCGATTATCAGGAGCAGGCGTCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGaa                                                                                                                                 <  2:3297842/148‑1 (MQ=255)
                 gatCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGaa                                                                                                                                 <  1:2842067/148‑1 (MQ=255)
                  atCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACTACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAg                                                                                                                                >  1:658538/1‑148 (MQ=255)
                  atCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAg                                                                                                                                <  1:2725490/148‑1 (MQ=255)
                       tttCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGc                                                                                                                           >  2:1124081/1‑148 (MQ=255)
                        ttCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCg                                                                                                                          <  1:1222544/148‑1 (MQ=255)
                        ttCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCg                                                                                                                          >  1:3218087/1‑148 (MQ=255)
                         tCGATTATCAGTAGCAGGCTGCGAATGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGa                                                                                                                         <  2:1363176/148‑1 (MQ=255)
                         tCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGa                                                                                                                         <  2:1983211/148‑1 (MQ=255)
                         tCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGa                                                                                                                         <  2:832700/148‑1 (MQ=255)
                         tCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGa                                                                                                                         <  1:426903/148‑1 (MQ=255)
                              tATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAag                                                                                                                    >  2:3188028/1‑148 (MQ=255)
                               aTCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAagc                                                                                                                   >  2:1407164/1‑148 (MQ=255)
                                tCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAagca                                                                                                                  <  1:2335807/148‑1 (MQ=255)
                                tCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAagca                                                                                                                  >  2:2061239/1‑148 (MQ=255)
                                    gAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGcagcag                                                                                                              >  2:610446/1‑148 (MQ=255)
                                      gCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCagag                                                                                                            <  1:1179961/148‑1 (MQ=255)
                                       ctggcggcGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGt                                                                                                           <  2:2983112/146‑1 (MQ=255)
                                       cAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGt                                                                                                           <  2:269672/148‑1 (MQ=255)
                                         ggcggcGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGtta                                                                                                         >  1:346794/1‑148 (MQ=255)
                                         ggcggcGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGtta                                                                                                         <  1:487916/148‑1 (MQ=255)
                                                aaaGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTAt                                                                                                  <  1:1321440/148‑1 (MQ=255)
                                                 aaGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTAtt                                                                                                 >  2:2741867/1‑148 (MQ=255)
                                                 aaGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTAtt                                                                                                 >  2:1084875/1‑148 (MQ=255)
                                                    cAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAgc                                                                                              >  2:1104481/1‑148 (MQ=255)
                                                    cAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAgc                                                                                              >  1:1045326/1‑148 (MQ=255)
                                                      gTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAgcgc                                                                                            <  2:1057314/148‑1 (MQ=255)
                                                      gTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTAATACTATTAgcgc                                                                                            >  1:2619409/1‑148 (MQ=255)
                                                        cGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCa                                                                                          <  1:2533284/148‑1 (MQ=255)
                                                        cGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTCTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCa                                                                                          <  1:97448/148‑1 (MQ=255)
                                                           cTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGgc                                                                                       >  1:94883/1‑148 (MQ=255)
                                                                tGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGt                                                                                  <  2:1519741/148‑1 (MQ=255)
                                                                    gcAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACaa                                                                              <  2:1291360/148‑1 (MQ=255)
                                                                    gcAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACaa                                                                              <  1:825470/148‑1 (MQ=255)
                                                                        aTGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCt                                                                          <  2:486162/148‑1 (MQ=255)
                                                                        aTGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCt                                                                          >  2:1234798/1‑148 (MQ=255)
                                                                        aTGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCt                                                                          <  1:2547584/148‑1 (MQ=255)
                                                                         tGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCtt                                                                         <  1:2290794/148‑1 (MQ=255)
                                                                          ggATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCttt                                                                        >  2:3055395/1‑148 (MQ=255)
                                                                           gATTCGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGTGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCtttt                                                                       <  1:1785291/148‑1 (MQ=255)
                                                                           gATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCtttt                                                                       <  2:642033/148‑1 (MQ=255)
                                                                           gATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCtttt                                                                       <  1:109841/148‑1 (MQ=255)
                                                                           gATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCtttt                                                                       >  1:691502/1‑148 (MQ=255)
                                                                           gATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCtttt                                                                       >  1:1928189/1‑148 (MQ=255)
                                                                           gATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCGACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCtttt                                                                       >  1:660583/1‑148 (MQ=255)
                                                                            aTACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTa                                                                      <  1:2402503/148‑1 (MQ=255)
                                                                             tACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTAt                                                                     <  2:570217/148‑1 (MQ=255)
                                                                              aCGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATg                                                                    <  1:2334728/148‑1 (MQ=255)
                                                                               cGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATgg                                                                   <  1:2517793/148‑1 (MQ=255)
                                                                                   cTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGtgat                                                               <  1:650886/148‑1 (MQ=255)
                                                                                      gcagcaAATGGCGGATAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGaa                                                            <  2:1926056/148‑1 (MQ=255)
                                                                                      gcagcaAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGaa                                                            <  2:706238/148‑1 (MQ=255)
                                                                                      gcagcaAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGaa                                                            >  2:593036/1‑148 (MQ=255)
                                                                                      gcagcaAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGaa                                                            <  2:1658893/148‑1 (MQ=255)
                                                                                      gcagcaAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGaa                                                            >  1:3061247/1‑148 (MQ=255)
                                                                                      gcagcaAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGaa                                                            >  1:2342269/1‑148 (MQ=255)
                                                                                          caAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCt                                                        <  2:2185963/148‑1 (MQ=255)
                                                                                           aaaTGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTc                                                       >  1:2678584/1‑148 (MQ=255)
                                                                                                   gagtttttGCAGTGTTTATACTTGTTGAGCTACATCTCGTCTCATTATTCCCGCCTGTTCCGTTTTGGGTTCGACAAAGAAGCAGAGTTATTACTATTAGCGCCAGGCGCTGTTCAAATGTTTTATTGTGATTAATACTCGCGTCAGt                                               <  2:1577551/142‑1 (MQ=255)
                                                                                                   gagaAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGt                                               <  2:1169000/148‑1 (MQ=255)
                                                                                                    agaAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGtt                                              >  2:38911/1‑148 (MQ=255)
                                                                                                        aTTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGt                                          >  2:2489868/1‑148 (MQ=255)
                                                                                                         ttGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTc                                         >  1:1559429/1‑148 (MQ=255)
                                                                                                              gggTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCATCGCGTCCGTTCGGTCCTACa                                    >  1:2180755/1‑148 (MQ=255)
                                                                                                               ggTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAg                                   >  1:1375699/1‑148 (MQ=255)
                                                                                                                 tttCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGt                                 >  1:2120895/1‑148 (MQ=255)
                                                                                                                 tttCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGt                                 >  1:2034599/1‑148 (MQ=255)
                                                                                                                  ttCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGtt                                >  1:3060362/1‑148 (MQ=255)
                                                                                                                  ttCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGtt                                <  2:2396973/148‑1 (MQ=255)
                                                                                                                  ttCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGtt                                <  2:658518/148‑1 (MQ=255)
                                                                                                                     aaCGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTg                             <  2:1800168/148‑1 (MQ=255)
                                                                                                                      aCGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGa                            >  1:677720/1‑148 (MQ=255)
                                                                                                                      aCGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGa                            <  2:3145406/148‑1 (MQ=255)
                                                                                                                       cGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACACATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGAt                           >  2:2164823/1‑148 (MQ=255)
                                                                                                                       cGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGAt                           >  2:2850833/1‑148 (MQ=255)
                                                                                                                       cGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGAt                           >  2:2916198/1‑148 (MQ=255)
                                                                                                                       cGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGAt                           >  1:3298201/1‑148 (MQ=255)
                                                                                                                       cGTGTTGAGCAACCTCTCGTCGCATGATACACGCCTTTACCGTGAAGGGGGAGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTAGGGTGATCAATCCTAGCGACACATAGGTCCTACAGGttatcat                           >  2:2850858/1‑143 (MQ=255)
                                                                                                                          gtTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCg                        <  1:172639/148‑1 (MQ=255)
                                                                                                                           ttGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGc                       <  1:728897/148‑1 (MQ=255)
                                                                                                                           ttGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGc                       <  1:1407164/148‑1 (MQ=255)
                                                                                                                                cTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCaa                  >  1:29165/1‑148 (MQ=255)
                                                                                                                                        cgtcgCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTACACg          <  2:692683/148‑1 (MQ=255)
                                                                                                                                        cgtcgCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTACACg          <  1:2934111/148‑1 (MQ=255)
                                                                                                                                          tcgCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTACACGtt        <  1:3129604/148‑1 (MQ=255)
                                                                                                                                           cgCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTACACGTTc       >  2:1335767/1‑148 (MQ=255)
                                                                                                                                           cgCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTACACGTTc       >  2:1335768/1‑148 (MQ=255)
                                                                                                                                           cgCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTACACGTTc       >  2:1637359/1‑148 (MQ=255)
                                                                                                                                             cATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTACACGTTCgg     <  2:941767/148‑1 (MQ=255)
                                                                                                                                             cATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTACACGTTCgg     <  2:2595238/148‑1 (MQ=255)
                                                                                                                                              aTGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTACACGTTCGGa    <  1:1149985/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CACGGCTTCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTACACGTTCGGATA  >  CP009273/3732405‑3732696

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