Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,872,054 A→T 100% E28V (GAG→GTG)  yeaP → diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,872,0540AT99.1% 746.4 / ‑2.5 227E28V (GAG→GTG) yeaPdiguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base T (102/122);  total (104/123)
Fisher's exact test for biased strand distribution p-value = 5.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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 228 total aligned reads displayed.

GTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTT  >  CP009273/1871909‑1872198
                                                                                                                                                 |                                                                                                                                                
  tttATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATg                                                                                                                                              <  1:1178979/148‑1 (MQ=255)
    tATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCt                                                                                                                                            <  2:96794/148‑1 (MQ=255)
    tATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCt                                                                                                                                            <  2:414735/148‑1 (MQ=255)
    tATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATNATAGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCt                                                                                                                                            >  1:1783504/1‑148 (MQ=255)
       cATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAACGTCTGGTCCGACAGCTTCTGGTGATGCTGGa                                                                                                                                         >  2:526777/1‑148 (MQ=255)
        aTGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGaa                                                                                                                                        <  2:3264212/148‑1 (MQ=255)
        aTGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGaa                                                                                                                                        <  2:2790154/148‑1 (MQ=255)
          gCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAAt                                                                                                                                      >  1:1712408/1‑148 (MQ=255)
             tGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGt                                                                                                                                   <  1:3090258/148‑1 (MQ=255)
                aaaGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCAc                                                                                                                                <  1:1576915/148‑1 (MQ=255)
                  aGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTg                                                                                                                              >  2:1809417/1‑148 (MQ=255)
                   gATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGa                                                                                                                             >  1:1328830/1‑148 (MQ=255)
                   gATTAAATAATCAGCATTTACCCGCCGTATCCAGGAGTTGTTCCGTGTCAGAACAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAACTTCTGGTGATGCTGGAAATGGTCACTGa                                                                                                                             >  2:1487546/1‑148 (MQ=255)
                    aTTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGat                                                                                                                            <  1:1654095/148‑1 (MQ=255)
                      tAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGatat                                                                                                                          >  2:1178934/1‑148 (MQ=255)
                           aatCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAAt                                                                                                                     <  1:2690216/148‑1 (MQ=255)
                           aatCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAAt                                                                                                                     <  2:706351/148‑1 (MQ=255)
                             tCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCa                                                                                                                   >  2:1156439/1‑148 (MQ=255)
                              cAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCaa                                                                                                                  >  2:353231/1‑148 (MQ=255)
                                gCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAacc                                                                                                                >  1:2980808/1‑148 (MQ=255)
                                gCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAacc                                                                                                                >  2:3309004/1‑148 (MQ=255)
                                gCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAacc                                                                                                                >  1:2219723/1‑148 (MQ=255)
                                   tttACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAacctac                                                                                                             >  1:194506/1‑148 (MQ=255)
                                   tttACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAacctac                                                                                                             >  1:2717298/1‑148 (MQ=255)
                                    ttACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGACAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAacctacc                                                                                                            >  2:3211222/1‑148 (MQ=255)
                                    ttACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAacctacc                                                                                                            >  2:1843160/1‑148 (MQ=255)
                                      aCCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTg                                                                                                          <  2:2356496/148‑1 (MQ=255)
                                        ccgccgTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGAc                                                                                                        >  1:2771272/1‑148 (MQ=255)
                                           ccgTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCaa                                                                                                     >  2:2595783/1‑148 (MQ=255)
                                             gTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAg                                                                                                   >  2:1068907/1‑148 (MQ=255)
                                              tATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGt                                                                                                  >  2:1792817/1‑148 (MQ=255)
                                               aTCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTg                                                                                                 <  2:2996703/148‑1 (MQ=255)
                                                tCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTgg                                                                                                >  1:223045/1‑148 (MQ=255)
                                                tCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTgg                                                                                                <  2:3432904/148‑1 (MQ=255)
                                                tCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGAGTCTCGCAATCCCTTGCCATAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTgg                                                                                                <  2:1986522/148‑1 (MQ=255)
                                                  cTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCTCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGAt                                                                                              <  2:2939472/148‑1 (MQ=255)
                                                    ggAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGt                                                                                            <  2:196552/148‑1 (MQ=255)
                                                    ggAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGt                                                                                            <  1:964838/148‑1 (MQ=255)
                                                    ggAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGt                                                                                            <  1:2559262/148‑1 (MQ=255)
                                                      aGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCg                                                                                          <  1:1518649/148‑1 (MQ=255)
                                                      aGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCg                                                                                          >  2:148429/1‑148 (MQ=255)
                                                           ttCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAgcg                                                                                     <  1:670088/148‑1 (MQ=255)
                                                              cGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAgcgcgc                                                                                  >  2:2556414/1‑148 (MQ=255)
                                                              cGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAgcgcgc                                                                                  <  1:753056/148‑1 (MQ=255)
                                                                 gtCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTg                                                                               >  1:239672/1‑148 (MQ=255)
                                                                  tcagatcagaTTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTgc                                                                              >  2:524917/1‑148 (MQ=255)
                                                                   cagatcagatTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTgca                                                                             >  2:1597075/1‑148 (MQ=255)
                                                                        cagatTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCata                                                                        >  2:1833233/1‑148 (MQ=255)
                                                                        cagatTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCata                                                                        <  2:2697459/148‑1 (MQ=255)
                                                                         agatTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCatat                                                                       <  1:2838671/148‑1 (MQ=255)
                                                                           atTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCAtatta                                                                     >  1:2473849/1‑148 (MQ=255)
                                                                              aTCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGt                                                                  <  1:422394/148‑1 (MQ=255)
                                                                               tCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGCTGTCGAAGCGCGCCTGCAGCATATTATGtt                                                                 >  2:503251/1‑148 (MQ=255)
                                                                               tCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGtt                                                                 <  2:796050/148‑1 (MQ=255)
                                                                                cGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGttt                                                                <  2:338961/148‑1 (MQ=255)
                                                                                  cccGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGc                                                              <  1:1801146/148‑1 (MQ=255)
                                                                                   ccGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGcc                                                             >  2:64065/1‑148 (MQ=255)
                                                                                        tctcGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTaa                                                        >  2:3100623/1‑148 (MQ=255)
                                                                                         ctcGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAAc                                                       <  2:217325/148‑1 (MQ=255)
                                                                                         ctcGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAAc                                                       >  1:3429653/1‑148 (MQ=255)
                                                                                         ctcGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAAc                                                       >  2:2276682/1‑148 (MQ=255)
                                                                                          tcGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACa                                                      <  2:836660/148‑1 (MQ=255)
                                                                                             cAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGcc                                                   >  1:3269396/1‑148 (MQ=255)
                                                                                               aTCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAg                                                 >  2:1542209/1‑148 (MQ=255)
                                                                                               aTCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAg                                                 >  1:1264619/1‑148 (MQ=255)
                                                                                                tCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGa                                                >  2:1730821/1‑148 (MQ=255)
                                                                                                tCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGa                                                <  2:334350/148‑1 (MQ=255)
                                                                                                tCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGa                                                <  1:3491111/148‑1 (MQ=255)
                                                                                                  ccTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGaaa                                              <  2:167202/148‑1 (MQ=255)
                                                                                                    ttGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAAt                                            >  1:75143/1‑148 (MQ=255)
                                                                                                      gCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGt                                          <  1:2187398/148‑1 (MQ=255)
                                                                                                          aaGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACAt                                      <  1:1949296/148‑1 (MQ=255)
                                                                                                          aaGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACAt                                      <  1:1430926/148‑1 (MQ=255)
                                                                                                          aaGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACAt                                      >  2:2753795/1‑148 (MQ=255)
                                                                                                           aGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATc                                     <  1:326241/148‑1 (MQ=255)
                                                                                                           aGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATc                                     >  1:1121408/1‑148 (MQ=255)
                                                                                                                   cACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGaa                             >  2:133331/1‑148 (MQ=255)
                                                                                                                     cTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAAtt                           >  2:2666079/1‑148 (MQ=255)
                                                                                                                     cTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAAtt                           >  2:1094071/1‑148 (MQ=255)
                                                                                                                      tGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAAttt                          >  2:3166977/1‑148 (MQ=255)
                                                                                                                      tGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAAttt                          <  2:1540313/148‑1 (MQ=255)
                                                                                                                       ggAACGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAAtttt                         <  1:526197/148‑1 (MQ=255)
                                                                                                                       ggAACGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAAtttt                         <  1:3179043/148‑1 (MQ=255)
                                                                                                                       ggAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAAtttt                         <  2:680722/148‑1 (MQ=255)
                                                                                                                       ggAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAAtttt                         <  1:265836/148‑1 (MQ=255)
                                                                                                                        gAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTa                        >  2:3340082/1‑148 (MQ=255)
                                                                                                                        gAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTa                        <  1:148436/148‑1 (MQ=255)
                                                                                                                        gAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTa                        <  1:1818996/148‑1 (MQ=255)
                                                                                                                         aaaGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTAc                       <  1:1987043/148‑1 (MQ=255)
                                                                                                                         aaaGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTAc                       <  2:426467/148‑1 (MQ=255)
                                                                                                                          aaGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTAcc                      >  1:1008331/1‑148 (MQ=255)
                                                                                                                          aaGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTAcc                      >  1:1008768/1‑148 (MQ=255)
                                                                                                                            gTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGt                    >  1:2825953/1‑148 (MQ=255)
                                                                                                                            gTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGt                    >  2:2844966/1‑148 (MQ=255)
                                                                                                                             tCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGtc                   <  1:960732/148‑1 (MQ=255)
                                                                                                                                    cGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTggg            <  2:1712408/148‑1 (MQ=255)
                                                                                                                                     gACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTCTGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGa           <  2:2471185/148‑1 (MQ=255)
                                                                                                                                      aCAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGAt          >  2:1207900/1‑148 (MQ=255)
                                                                                                                                        aGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGATTa        <  1:3246385/148‑1 (MQ=255)
                                                                                                                                          cTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGATTACt      <  2:919657/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTT  >  CP009273/1871909‑1872198

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: