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% 754.8 / ‑2.5 230E28V (GAG→GTG) yeaPdiguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base T (130/98);  total (130/100)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 235 total aligned reads displayed.

GTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTT  >  CP009273/1871909‑1872198
                                                                                                                                                 |                                                                                                                                                
gTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGgtga                                                                                                                                                <  2:3736983/148‑1 (MQ=255)
gTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGgtga                                                                                                                                                <  2:1165445/148‑1 (MQ=255)
gTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGgtga                                                                                                                                                >  1:1483020/1‑148 (MQ=255)
 ttttATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGAt                                                                                                                                               >  1:3239161/1‑148 (MQ=255)
   ttATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGc                                                                                                                                             <  1:1077983/148‑1 (MQ=255)
    tATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCt                                                                                                                                            >  2:2652452/1‑148 (MQ=255)
      tCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTgg                                                                                                                                          >  1:2764072/1‑148 (MQ=255)
           ccTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATg                                                                                                                                     <  2:121822/148‑1 (MQ=255)
              gCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTc                                                                                                                                  >  1:2791705/1‑148 (MQ=255)
               cAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCa                                                                                                                                 <  2:2612866/148‑1 (MQ=255)
               cAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCa                                                                                                                                 <  2:3276844/148‑1 (MQ=255)
               cAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCa                                                                                                                                 <  2:531280/148‑1 (MQ=255)
                  aGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTg                                                                                                                              <  1:3452202/148‑1 (MQ=255)
                  aGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTg                                                                                                                              >  2:2639354/1‑148 (MQ=255)
                   gATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGa                                                                                                                             >  1:2454178/1‑148 (MQ=255)
                            atCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATc                                                                                                                    >  1:1582481/1‑148 (MQ=255)
                             tCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCa                                                                                                                   >  1:3735003/1‑148 (MQ=255)
                              cAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCaa                                                                                                                  >  2:2936312/1‑148 (MQ=255)
                              cAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCaa                                                                                                                  >  1:3769629/1‑148 (MQ=255)
                              cAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCaa                                                                                                                  >  2:673517/1‑148 (MQ=255)
                               aGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAac                                                                                                                 >  1:3674510/1‑148 (MQ=255)
                               aGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAac                                                                                                                 >  1:556833/1‑148 (MQ=255)
                                gCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAacc                                                                                                                >  2:343875/1‑148 (MQ=255)
                                 cATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAacct                                                                                                               <  2:873874/148‑1 (MQ=255)
                                 cATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAacct                                                                                                               <  2:2787727/148‑1 (MQ=255)
                                   tttACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAacctac                                                                                                             >  2:3113626/1‑148 (MQ=255)
                                       cccGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGa                                                                                                         <  2:3725420/148‑1 (MQ=255)
                                       cccGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGa                                                                                                         <  2:3725230/148‑1 (MQ=255)
                                        ccgccgTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGAc                                                                                                        >  1:1990175/1‑148 (MQ=255)
                                         cgccgTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGACTCAACCTACCTGAcc                                                                                                       >  2:2718100/1‑148 (MQ=255)
                                             gTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAAACTACCTGACCAAAg                                                                                                   >  1:2952676/1‑148 (MQ=255)
                                                tCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACTAAAGTgg                                                                                                >  1:2495808/1‑148 (MQ=255)
                                                   tGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATg                                                                                             >  1:1436355/1‑148 (MQ=255)
                                                     gAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTc                                                                                           >  1:231529/1‑148 (MQ=255)
                                                     gAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTc                                                                                           >  2:2369128/1‑148 (MQ=255)
                                                        ttgttCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGaa                                                                                        >  1:2116687/1‑148 (MQ=255)
                                                             ccGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAgcgcg                                                                                   <  2:740892/148‑1 (MQ=255)
                                                              cGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAgcgcgc                                                                                  <  1:1916595/148‑1 (MQ=255)
                                                                  tcagatcagaTTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTgc                                                                              >  1:646692/1‑148 (MQ=255)
                                                                    agatcagatTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTgcag                                                                            >  2:140061/1‑148 (MQ=255)
                                                                     gatcagatTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTgcagc                                                                           >  2:196994/1‑148 (MQ=255)
                                                                     gatcagatTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTgcagc                                                                           <  1:2096542/148‑1 (MQ=255)
                                                                        cagatTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCata                                                                        >  1:975957/1‑148 (MQ=255)
                                                                               tCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGtt                                                                 <  1:1057330/148‑1 (MQ=255)
                                                                                   ccGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGcc                                                             <  2:1560344/148‑1 (MQ=255)
                                                                                   ccGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGcc                                                             <  2:2276027/148‑1 (MQ=255)
                                                                                      cgTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGt                                                          >  1:3049605/1‑148 (MQ=255)
                                                                                       gTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTa                                                         <  1:2203584/148‑1 (MQ=255)
                                                                                            gCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGc                                                    <  2:1855001/148‑1 (MQ=255)
                                                                                               aTCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAg                                                 >  2:3468014/1‑148 (MQ=255)
                                                                                               aTCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAg                                                 >  1:2550997/1‑148 (MQ=255)
                                                                                               aTCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAg                                                 >  1:1310778/1‑148 (MQ=255)
                                                                                               aTCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAg                                                 >  2:1997893/1‑148 (MQ=255)
                                                                                                 cccTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGaa                                               >  1:380527/1‑148 (MQ=255)
                                                                                                    ttGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAAt                                            <  1:1447854/148‑1 (MQ=255)
                                                                                                    ttGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAAt                                            >  2:3354149/1‑148 (MQ=255)
                                                                                                     tGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATg                                           >  1:1762416/1‑148 (MQ=255)
                                                                                                      gCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGt                                          >  2:1636485/1‑148 (MQ=255)
                                                                                                      gCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGt                                          >  2:1470945/1‑148 (MQ=255)
                                                                                                       ccAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTa                                         >  1:45990/1‑148 (MQ=255)
                                                                                                       ccAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTa                                         >  1:3466908/1‑148 (MQ=255)
                                                                                                        cAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTAc                                        <  1:315227/148‑1 (MQ=255)
                                                                                                         aaaGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACa                                       >  1:1628126/1‑148 (MQ=255)
                                                                                                          aaGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACAt                                      >  2:2563007/1‑148 (MQ=255)
                                                                                                           aGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATc                                     <  1:139795/148‑1 (MQ=255)
                                                                                                            gAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATcc                                    >  1:1728728/1‑148 (MQ=255)
                                                                                                             aaCAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATccc                                   >  1:1240593/1‑148 (MQ=255)
                                                                                                             aaCAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATccc                                   >  1:2911675/1‑148 (MQ=255)
                                                                                                             aaCAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATccc                                   >  1:1515610/1‑148 (MQ=255)
                                                                                                             aaCAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATccc                                   <  1:2362722/148‑1 (MQ=255)
                                                                                                              aCAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAACCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCg                                  >  1:893142/1‑148 (MQ=255)
                                                                                                                   cACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGaa                             <  2:231529/148‑1 (MQ=255)
                                                                                                                    aCTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAAt                            <  2:512430/148‑1 (MQ=255)
                                                                                                                    aCTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAAt                            <  2:2434590/148‑1 (MQ=255)
                                                                                                                      tGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAAttt                          <  1:1050465/148‑1 (MQ=255)
                                                                                                                       ggAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAAtttt                         >  1:1125813/1‑148 (MQ=255)
                                                                                                                        gAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTa                        >  2:689012/1‑148 (MQ=255)
                                                                                                                        gAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTa                        <  1:3390588/148‑1 (MQ=255)
                                                                                                                        gAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTa                        >  2:1267648/1‑148 (MQ=255)
                                                                                                                         aaaGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTAc                       >  2:3800959/1‑148 (MQ=255)
                                                                                                                           aGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCg                     >  2:1420892/1‑148 (MQ=255)
                                                                                                                             tCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGtc                   >  1:1651983/1‑148 (MQ=255)
                                                                                                                             tCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGtc                   >  1:1652287/1‑148 (MQ=255)
                                                                                                                                ggTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCg                <  2:131612/148‑1 (MQ=255)
                                                                                                                                ggTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCg                <  2:2339244/148‑1 (MQ=255)
                                                                                                                                 gTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGt               >  2:1311241/1‑148 (MQ=255)
                                                                                                                                 gTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGt               >  1:827056/1‑148 (MQ=255)
                                                                                                                                 gTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGt               >  1:718905/1‑148 (MQ=255)
                                                                                                                                 gTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGt               >  2:3585776/1‑148 (MQ=255)
                                                                                                                                  tCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTg              >  1:3367130/1‑148 (MQ=255)
                                                                                                                                  tCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTg              <  1:2080905/148‑1 (MQ=255)
                                                                                                                                     gACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGa           <  1:3611584/148‑1 (MQ=255)
                                                                                                                                     gACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAACTTTACCGTCTCGTGGGa           <  2:340685/148‑1 (MQ=255)
                                                                                                                                      aCAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGAt          >  2:928218/1‑148 (MQ=255)
                                                                                                                                      aCAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGAt          <  1:2050785/148‑1 (MQ=255)
                                                                                                                                      aCAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGAt          <  1:2307516/148‑1 (MQ=255)
                                                                                                                                          cTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGATTACt      <  2:2586448/148‑1 (MQ=255)
                                                                                                                                            tCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCGAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCg    <  1:3198205/148‑1 (MQ=255)
                                                                                                                                             cTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGt   <  2:3231286/148‑1 (MQ=255)
                                                                                                                                             cTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGt   <  1:3571811/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: ‑