Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,139,910 G→T 100% R1055L (CGA→CTA)  yegE → putative diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,139,9100GT98.8% 804.7 / 0.4 247R1055L (CGA→CTA) yegEputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base T (104/140);  total (104/143)
Fisher's exact test for biased strand distribution p-value = 2.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-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 247 total aligned reads displayed.

CCAGGTGGGCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTATTTCGCG  >  CP009273/2139764‑2140055
                                                                                                                                                  |                                                                                                                                                 
ccAGGTGGTCAGCGATCTTCAAATCTTCAACTCGCTGAAAGCGAATATTTCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTAc                                                                                                                                                 <  2:1374858/149‑1 (MQ=255)
ccAGGTGGGCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTAc                                                                                                                                                 <  2:507108/149‑1 (MQ=255)
  aGGTGGGCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATCATTCAGGGGCACGCTCAGCTACTc                                                                                                                                               <  1:3933134/149‑1 (MQ=255)
   ggTGGGCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCg                                                                                                                                              >  2:2663233/1‑149 (MQ=255)
      gggCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGa                                                                                                                                           >  2:3919935/1‑149 (MQ=255)
            cgATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGaaaa                                                                                                                                     >  1:3865954/1‑149 (MQ=255)
            cgATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGaaaa                                                                                                                                     <  2:1360256/149‑1 (MQ=255)
             gATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAAc                                                                                                                                    >  1:1468764/1‑149 (MQ=255)
              aTCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAAcc                                                                                                                                   <  1:162538/149‑1 (MQ=255)
               tCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCa                                                                                                                                  <  2:29822/149‑1 (MQ=255)
                cTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCAt                                                                                                                                 <  1:562525/149‑1 (MQ=255)
                 tGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATc                                                                                                                                >  1:1173308/1‑149 (MQ=255)
                   cAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGATTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACTATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGc                                                                                                                              <  2:3682297/149‑1 (MQ=255)
                    aaaTCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGcc                                                                                                                             >  1:4002365/1‑149 (MQ=255)
                     aaTCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCg                                                                                                                            >  1:195228/1‑149 (MQ=255)
                      aTCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCgg                                                                                                                           >  1:447090/1‑149 (MQ=255)
                      aTCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCgg                                                                                                                           >  2:1599525/1‑149 (MQ=255)
                         ttCCACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGcc                                                                                                                        <  2:3908927/149‑1 (MQ=255)
                           cAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAg                                                                                                                      >  1:1782101/1‑149 (MQ=255)
                             aCTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTc                                                                                                                    <  1:3475279/149‑1 (MQ=255)
                               tCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGt                                                                                                                  >  1:1331622/1‑149 (MQ=255)
                                cGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGtt                                                                                                                 <  2:2007780/149‑1 (MQ=255)
                                cGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGtt                                                                                                                 <  2:1906521/149‑1 (MQ=255)
                                cGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGtt                                                                                                                 <  2:3025726/149‑1 (MQ=255)
                                    gAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTAc                                                                                                             >  2:3809839/1‑149 (MQ=255)
                                     aaaGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTAcc                                                                                                            <  2:2354928/149‑1 (MQ=255)
                                      aaGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTGCTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTAccc                                                                                                           >  2:174942/1‑149 (MQ=255)
                                       aGCGAATATGGCAGATTACCTGCTACTTGATTGTGAGTTATGCGCCAACGTGCAGTGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGAACGATCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCt                                                                                                          <  2:839111/149‑1 (MQ=255)
                                       aGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCt                                                                                                          >  1:1437625/1‑149 (MQ=255)
                                           aaTATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGt                                                                                                      <  1:702309/149‑1 (MQ=255)
                                            atatGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTg                                                                                                     <  2:4049555/149‑1 (MQ=255)
                                              atGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGAt                                                                                                   <  2:2505468/149‑1 (MQ=255)
                                                ggCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATgg                                                                                                 >  1:4089862/1‑149 (MQ=255)
                                                ggCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATgg                                                                                                 >  2:2646417/1‑149 (MQ=255)
                                                 gCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCTGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGa                                                                                                <  2:847312/149‑1 (MQ=255)
                                                 gCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGa                                                                                                <  2:282662/149‑1 (MQ=255)
                                                 gCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGa                                                                                                >  2:3613731/1‑149 (MQ=255)
                                                 gCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGa                                                                                                >  1:321591/1‑149 (MQ=255)
                                                  cAGATTACCTGCTACTTTATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGAt                                                                                               <  2:2585552/149‑1 (MQ=255)
                                                   aGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATa                                                                                              <  2:4090466/149‑1 (MQ=255)
                                                   aGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATa                                                                                              >  2:3263129/1‑149 (MQ=255)
                                                     aTTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATaag                                                                                            >  1:1717852/1‑147 (MQ=255)
                                                     aTTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACg                                                                                            >  1:901978/1‑149 (MQ=255)
                                                       tACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCt                                                                                          >  1:1028864/1‑149 (MQ=255)
                                                       tACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCt                                                                                          <  2:3863019/149‑1 (MQ=255)
                                                        aCCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCtt                                                                                         >  2:3402265/1‑149 (MQ=255)
                                                         ccTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCttt                                                                                        <  1:2236409/149‑1 (MQ=255)
                                                           tGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCt                                                                                      <  1:173916/149‑1 (MQ=255)
                                                           tGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCt                                                                                      <  1:464110/149‑1 (MQ=255)
                                                            gCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTg                                                                                     <  1:3185905/149‑1 (MQ=255)
                                                                  tGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCg                                                                               <  1:2602488/149‑1 (MQ=255)
                                                                   gATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTAACCTTAGTGTTGGATACGCTTTCTGGCATCgg                                                                              >  1:1912195/1‑149 (MQ=255)
                                                                     tGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCg                                                                            >  1:1983973/1‑149 (MQ=255)
                                                                       gTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTc                                                                          <  1:1590708/149‑1 (MQ=255)
                                                                           gTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATc                                                                      <  2:2617783/149‑1 (MQ=255)
                                                                            ttATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCt                                                                     <  2:3423657/149‑1 (MQ=255)
                                                                             tATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTg                                                                    <  2:155229/149‑1 (MQ=255)
                                                                               tGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCGTAGTGATGGATACGCTTTCTGGCCTCGGCGTCGATCTGAt                                                                  >  2:2471970/1‑149 (MQ=255)
                                                                                gcgcCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGAtt                                                                 <  1:2973593/149‑1 (MQ=255)
                                                                                   ccAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATAAGCTTTCTGGCATCGGCGTCGATCTGATTTAt                                                              <  2:1717852/149‑1 (MQ=255)
                                                                                      aCGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGCTCTGATTTATggt                                                           >  2:1496223/1‑149 (MQ=255)
                                                                                      aCGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATggt                                                           <  1:2946040/149‑1 (MQ=255)
                                                                                        gTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGTTTTATggtga                                                         >  2:1083517/1‑149 (MQ=255)
                                                                                         tGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATggtgag                                                        >  2:3487863/1‑149 (MQ=255)
                                                                                           cAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATggtgaggt                                                      <  1:3307298/149‑1 (MQ=255)
                                                                                            aGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTTCCCGTAGTGATGGATACGCTTTCTGGCCTCGGCGGCGCTCTGATTTTTggtgaggtg                                                     >  1:1823879/1‑149 (MQ=255)
                                                                                            aGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATggtgaggtg                                                     <  1:2071091/149‑1 (MQ=255)
                                                                                            aGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATggtgaggtg                                                     <  1:4104672/149‑1 (MQ=255)
                                                                                            aGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCGTAGTGATGGATACGCTTTGTGGCAAAGGCGTCGATCGGTTTTATtgtgtggtg                                                     >  1:1824497/1‑149 (MQ=255)
                                                                                             gggTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTATTGATTGATAAGCTTTCTGGCATCGGCGTCGAGGTGGTTTTTtgtgtggtgt                                                    >  1:769587/1‑148 (MQ=255)
                                                                                              ggTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGAt                                                   >  1:3416069/1‑149 (MQ=255)
                                                                                              ggTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGAt                                                   <  1:497040/149‑1 (MQ=255)
                                                                                                 aaTTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGc                                                >  2:737434/1‑149 (MQ=255)
                                                                                                 aaTTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGc                                                <  1:268041/149‑1 (MQ=255)
                                                                                                    ttGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGa                                             <  2:2124213/149‑1 (MQ=255)
                                                                                                     tGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGAt                                            <  2:3799839/149‑1 (MQ=255)
                                                                                                      gatggatgAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATg                                           >  1:2289240/1‑149 (MQ=255)
                                                                                                             gagaTGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGTTGGATACGCTTTCTGGCATCGGCGGCGATCTGGTTTATGGGGAGGTGTTTGCCGATGCCCCAcc                                    >  1:129551/1‑149 (MQ=255)
                                                                                                             gagaTGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATAAGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGCTTCCCAAcc                                    >  1:2204941/1‑149 (MQ=255)
                                                                                                               gaTGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGc                                  <  1:3441961/149‑1 (MQ=255)
                                                                                                                aTGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCt                                 <  1:2646474/149‑1 (MQ=255)
                                                                                                                 tGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTg                                >  1:3377264/1‑149 (MQ=255)
                                                                                                                  gCTGATTACGATTATACAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTgg                               >  2:932281/1‑149 (MQ=255)
                                                                                                                    tGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGAt                             <  1:2478006/149‑1 (MQ=255)
                                                                                                                     gATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGGCGATCTGATTTATGGTGAGGTGATTTCCGATGCCCAACCGCTGGAtt                            >  2:3257497/1‑149 (MQ=255)
                                                                                                                     gATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCGTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGAtt                            >  1:2583249/1‑149 (MQ=255)
                                                                                                                      aTTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGAttt                           <  2:901978/149‑1 (MQ=255)
                                                                                                                        tACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGc                         <  1:1075366/149‑1 (MQ=255)
                                                                                                                        tACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGc                         <  2:3046988/149‑1 (MQ=255)
                                                                                                                         aCGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTATTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATTGTGAGGTGATTGCCGATGCCCAAACGCTGGATTTGCt                        >  2:2703871/1‑149 (MQ=255)
                                                                                                                         aCGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCt                        <  1:2953506/149‑1 (MQ=255)
                                                                                                                           gATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTgg                      <  1:3724342/149‑1 (MQ=255)
                                                                                                                            attattCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGt                     <  2:2289240/149‑1 (MQ=255)
                                                                                                                            attattCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGt                     >  2:3566750/1‑149 (MQ=255)
                                                                                                                                    ggggCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAAtagta             <  2:2805801/149‑1 (MQ=255)
                                                                                                                                    ggggCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAAtagta             <  1:1796315/149‑1 (MQ=255)
                                                                                                                                       gCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGtt          <  2:413460/149‑1 (MQ=255)
                                                                                                                                       gCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGtt          <  2:66623/149‑1 (MQ=255)
                                                                                                                                       gCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTAACCTTAGTGATGGATACGCGTTCTGGCATCGGCGTCGATGTGATTTATGGTGAGGTGTTTGGCGATTCCCCACCCCTGGGTTTGCTGGTGAAAATTAGtt          >  1:2175718/1‑149 (MQ=255)
                                                                                                                                          cGCTCAGCTACTCGTGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTAtt       <  1:3099981/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CCAGGTGGGCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTATTTCGCG  >  CP009273/2139764‑2140055

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: