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,9100GT99.3% 1041.6 / ‑1.8 303R1055L (CGA→CTA) yegEputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (165/136);  total (166/137)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
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 303 total aligned reads displayed.

AGGTGGGCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTATTTCGC  >  CP009273/2139766‑2140054
                                                                                                                                                |                                                                                                                                                
aGGTGGGCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTc                                                                                                                                              >  2:2820554/1‑149 (MQ=255)
aGGTGGGCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTc                                                                                                                                              >  2:262733/1‑149 (MQ=255)
      gCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATg                                                                                                                                        <  2:4581304/149‑1 (MQ=255)
                gCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTATGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCg                                                                                                                              >  2:1731526/1‑149 (MQ=255)
                 cAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGc                                                                                                                             <  1:3600073/149‑1 (MQ=255)
                   aaTCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCg                                                                                                                           >  1:4749091/1‑149 (MQ=255)
                   aaTCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCg                                                                                                                           <  2:1232996/149‑1 (MQ=255)
                   aaTCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCg                                                                                                                           <  2:1498032/149‑1 (MQ=255)
                   aaTCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCg                                                                                                                           >  2:3876586/1‑149 (MQ=255)
                   aaTCTTCAACTCGCTGAAAGCGAAAATGGCAGATTACCAGCTACTTGAAGGTGAGATATGCGCCAACGTGCAGGGTAATATGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCAACGCCg                                                                                                                           >  1:4430624/1‑149 (MQ=255)
                    aTCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCgg                                                                                                                          <  1:545556/149‑1 (MQ=255)
                    aTCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCgg                                                                                                                          <  1:3861276/149‑1 (MQ=255)
                     tCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCggg                                                                                                                         >  2:2704006/1‑149 (MQ=255)
                      cTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGc                                                                                                                        >  2:3513032/1‑149 (MQ=255)
                      cTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGc                                                                                                                        <  1:2732189/149‑1 (MQ=255)
                         cAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAg                                                                                                                     <  2:3225022/149‑1 (MQ=255)
                           aCTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTc                                                                                                                   >  1:3326046/1‑149 (MQ=255)
                             tCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGt                                                                                                                 <  1:924793/149‑1 (MQ=255)
                              cGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGtt                                                                                                                <  2:3973112/149‑1 (MQ=255)
                                cTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGtttt                                                                                                              >  2:2402283/1‑149 (MQ=255)
                                cTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGtttt                                                                                                              <  2:533064/149‑1 (MQ=255)
                                 tGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTa                                                                                                             >  2:2986775/1‑149 (MQ=255)
                                  gAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTAc                                                                                                            <  2:2220967/149‑1 (MQ=255)
                                  gAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTAc                                                                                                            >  2:1841742/1‑149 (MQ=255)
                                  gAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTAc                                                                                                            >  1:3013195/1‑149 (MQ=255)
                                         aaTATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGt                                                                                                     <  2:1030209/149‑1 (MQ=255)
                                          atatGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTg                                                                                                    >  1:208538/1‑149 (MQ=255)
                                          atatGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTg                                                                                                    <  2:1024660/149‑1 (MQ=255)
                                            atGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGAt                                                                                                  >  2:3749825/1‑149 (MQ=255)
                                                cAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGAt                                                                                              <  2:1300501/149‑1 (MQ=255)
                                                  gATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATAc                                                                                            >  2:3941575/1‑149 (MQ=255)
                                                   aTTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACg                                                                                           <  2:2900852/149‑1 (MQ=255)
                                                    ttACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGc                                                                                          <  2:2367404/149‑1 (MQ=255)
                                                        cTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTc                                                                                      <  1:1939716/149‑1 (MQ=255)
                                                         tGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCt                                                                                     <  2:3097215/149‑1 (MQ=255)
                                                          gCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATAAGCTTTCTg                                                                                    >  1:1578506/1‑149 (MQ=255)
                                                             aCTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCa                                                                                 >  2:1984473/1‑149 (MQ=255)
                                                              cTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCAt                                                                                >  2:334011/1‑149 (MQ=255)
                                                               ttGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTAACCTTAGTGATGGATAAGCTTTATTGCATc                                                                               >  1:899527/1‑149 (MQ=255)
                                                               ttGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTGTTAACCTTAGTGTTGGATAAGCTTTATTGCATc                                                                               >  1:1718788/1‑149 (MQ=255)
                                                                    ggTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTGTTTCCCGTAGTGATGGATACGCTTTCTGGCATCGGCGt                                                                          >  1:813509/1‑149 (MQ=255)
                                                                             tGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCGTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGAt                                                                 >  2:2027910/1‑149 (MQ=255)
                                                                              gcgcCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGAtt                                                                <  1:4465259/149‑1 (MQ=255)
                                                                              gcgcCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGAtt                                                                <  1:2268821/149‑1 (MQ=255)
                                                                              gcgcCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGAtt                                                                <  2:1435193/149‑1 (MQ=255)
                                                                               cgcCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGAttt                                                               <  1:67668/149‑1 (MQ=255)
                                                                                 ccAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTAt                                                             <  1:1999398/149‑1 (MQ=255)
                                                                                   aaCGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTGTTACCCTTAGTTATGGATACGCTTTCTGGCATCGGCGTCGATCTTATTTATgg                                                           >  1:1342334/1‑149 (MQ=255)
                                                                                    aCGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCGTTCTGGCATCGGCGTCGATATGATTTATggt                                                          >  2:1391267/1‑149 (MQ=255)
                                                                                      gTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATggtga                                                        <  2:4432840/149‑1 (MQ=255)
                                                                                         cAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATggtgaggt                                                     >  2:4300362/1‑149 (MQ=255)
                                                                                         cAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATggtgaggt                                                     <  1:1250025/149‑1 (MQ=255)
                                                                                           gggTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATggtgaggtgt                                                   >  1:2417056/1‑148 (MQ=255)
                                                                                                   tGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGAt                                           <  2:2594456/149‑1 (MQ=255)
                                                                                                    gatggatgAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTTGCATCGGCGTCGATATGATTTATTGTGAGGTGATTTCCCATg                                          >  1:3842491/1‑149 (MQ=255)
                                                                                                      tggatgAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGcc                                        <  2:1014378/149‑1 (MQ=255)
                                                                                                       ggatgAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGccc                                       >  2:698502/1‑149 (MQ=255)
                                                                                                       ggatgAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGccc                                       <  2:3218518/149‑1 (MQ=255)
                                                                                                        gatgAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCa                                      >  2:2513066/1‑149 (MQ=255)
                                                                                                        gatgAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGAAGCCCa                                      >  2:1528696/1‑149 (MQ=255)
                                                                                                        gatgAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTGTTACCCGTAGTGCTGGATACGCTTTCTGGCCTCGGCGTCGATCGGATTGATGGTGAGGGGATTGCCGATGcccc                                      >  1:870430/1‑148 (MQ=255)
                                                                                                            agaTGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCg                                  >  1:3506118/1‑149 (MQ=255)
                                                                                                            agaTGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCATAGTGATGGATAAGCTTTAATGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATTCCCaaacc                                  >  1:4564709/1‑146 (MQ=255)
                                                                                                            agaTGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTAACCATAGTGATGGATACGCTTTCTGGCATAGGCGTCGATATTATTTATGGTGAGGTGATTGCCGATGCCCAACCg                                  >  1:1116817/1‑149 (MQ=255)
                                                                                                              aTGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCt                                <  2:3623700/149‑1 (MQ=255)
                                                                                                               tGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTGTTACCCTTAGTGATGGATACGCTTTATGGCATCGGCGTCGATCTGATTTATTGTGAGGTGATTGCCGATGCCCAACCGCTg                               >  1:3878321/1‑149 (MQ=255)
                                                                                                                 cTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGa                             >  2:3140425/1‑149 (MQ=255)
                                                                                                                  tGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGAt                            <  1:3183566/149‑1 (MQ=255)
                                                                                                                    aTTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTTCCCGTAGTGATGGATAAGCGTTATGGCATCGGCGTCGGTCTGTTTTTTGGTGGGGTGGTTTCCGATGCCCAACCCCAGGAttt                          >  1:4403732/1‑149 (MQ=255)
                                                                                                                     ttACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCGTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGCTGCCCCACCGCTGGATTTg                         >  2:1802861/1‑149 (MQ=255)
                                                                                                                      tACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGc                        <  1:3346830/149‑1 (MQ=255)
                                                                                                                      tACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGc                        <  1:2529814/149‑1 (MQ=255)
                                                                                                                      tACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGc                        <  2:2319477/149‑1 (MQ=255)
                                                                                                                       aCGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCt                       <  1:2986449/149‑1 (MQ=255)
                                                                                                                        cGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGTTTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTg                      >  1:4305938/1‑149 (MQ=255)
                                                                                                                        cGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTg                      >  1:3542646/1‑149 (MQ=255)
                                                                                                                         gATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTgg                     <  2:2937906/149‑1 (MQ=255)
                                                                                                                         gATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTgg                     >  1:4369673/1‑149 (MQ=255)
                                                                                                                          attattCGGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGt                    <  1:4043777/149‑1 (MQ=255)
                                                                                                                          ataattCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTGTTACCCGTAGTGATGGATACGCTTTCTGGCAGCGGCGTCGATCGGATTTATGGTGTGGTGATTTGCGATGCCCCAACCCTGGAGTTGCTggg                    >  1:2249411/4‑148 (MQ=255)
                                                                                                                           ttattCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTTCCGATGCCCAACCGCTGGATTTGCTgggg                   >  1:3556101/1‑147 (MQ=255)
                                                                                                                           ttattCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTg                   >  1:2383396/1‑149 (MQ=255)
                                                                                                                           ttattCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTAGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTg                   <  2:4488934/149‑1 (MQ=255)
                                                                                                                           ttattCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCGTAGTGATGGATACGCTTTCTGGCATCGGCGTCCATCAGATTTAAGGTGAGGTGATTGCCGATGACCAACCGCTGGATTTTCTGGTg                   >  2:2155233/1‑149 (MQ=255)
                                                                                                                            tattCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGCTCTGATTTATGGGGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGa                  >  2:4661809/1‑149 (MQ=255)
                                                                                                                            tattCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCGTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGa                  >  1:3950945/1‑149 (MQ=255)
                                                                                                                             attCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTGTTACCCGTAGTGATGGATACGCGTTATGGCATCGGCGTCGATCTGTTTTATTGTGTGGTGTTTTCCGATGCCCAACCCCAGGATTTGCTggggta                 >  1:3255097/1‑145 (MQ=255)
                                                                                                                                cAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAAtag              <  2:3609394/149‑1 (MQ=255)
                                                                                                                                 aGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAAtagt             >  2:1873031/1‑149 (MQ=255)
                                                                                                                                 aGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAAtagt             <  1:1654781/149‑1 (MQ=255)
                                                                                                                                   gggCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAAtagtag           <  2:476731/149‑1 (MQ=255)
                                                                                                                                    ggCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATagtagt          <  2:3389480/149‑1 (MQ=255)
                                                                                                                                       aCGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTAt       <  1:2820554/149‑1 (MQ=255)
                                                                                                                                       aCGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTAt       <  1:3905953/149‑1 (MQ=255)
                                                                                                                                       aCGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCCGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTAt       <  2:1951249/149‑1 (MQ=255)
                                                                                                                                        cGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTAtt      <  2:590030/149‑1 (MQ=255)
                                                                                                                                          cTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTATTTc    <  1:3073869/149‑1 (MQ=255)
                                                                                                                                          cTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCGTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATTGTGTGGGGATTGCCGATGCCCAACCGCTGGATTTTCGGGGGAATATTAGTTAttta    >  1:3879109/1‑148 (MQ=255)
                                                                                                                                           tCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTATTTcg   <  2:4064440/149‑1 (MQ=255)
                                                                                                                                           tCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTATTTcg   >  2:1083961/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                
AGGTGGGCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTATTTCGC  >  CP009273/2139766‑2140054

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