Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,138,039 C→A 74.1% I431I (ATC→ATA yegE → putative diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,138,0390CA74.1% 450.6 / 277.4 335I431I (ATC→ATAyegEputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (43/44);  new base A (99/149);  total (142/193)
Fisher's exact test for biased strand distribution p-value = 1.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

Only 100 of 338 total aligned reads displayed.

AACCGCAATGGCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAGTGCGTGCTGCC  >  CP009273/2137890‑2138188
                                                                                                                                                     |                                                                                                                                                     
aaCCGCAATGGCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATc                                                                                                                                                       >  2:1787485/1‑150 (MQ=255)
 aCCGCAATGGCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCa                                                                                                                                                      <  2:3765714/150‑1 (MQ=255)
         ggCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCt                                                                                                                                              >  2:1248958/1‑150 (MQ=255)
            gATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAAc                                                                                                                                           >  1:1002906/1‑150 (MQ=255)
            gATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAAc                                                                                                                                           >  2:1261293/1‑150 (MQ=255)
             aTGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACg                                                                                                                                          <  2:4637461/150‑1 (MQ=255)
                  gtCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAgcg                                                                                                                                     >  1:4554200/1‑150 (MQ=255)
                     tGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAgcgggc                                                                                                                                  <  2:4395291/150‑1 (MQ=255)
                      gggCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAgcgggcg                                                                                                                                 >  2:2085902/1‑150 (MQ=255)
                      gggCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCACATAACTCTGGCTAACGAAgcgggcg                                                                                                                                 >  2:471086/1‑150 (MQ=255)
                        gCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGcgggcggg                                                                                                                               <  2:5305002/150‑1 (MQ=255)
                          gTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGAt                                                                                                                             <  2:1178996/150‑1 (MQ=255)
                          gTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGAt                                                                                                                             <  1:2234199/150‑1 (MQ=255)
                          gTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGAt                                                                                                                             <  1:4429710/150‑1 (MQ=255)
                            tgcttgcCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTg                                                                                                                           <  1:301038/150‑1 (MQ=255)
                              cttgcCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGc                                                                                                                         >  1:4866439/1‑150 (MQ=255)
                              cttgcCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGc                                                                                                                         <  2:2863229/150‑1 (MQ=255)
                                 gcagTCTCACTGGTGCGCCACTCGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATc                                                                                                                      <  1:941126/147‑1 (MQ=255)
                                  ccGTCTCTCTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCt                                                                                                                     <  2:5245634/150‑1 (MQ=255)
                                  ccGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCt                                                                                                                     <  1:4514430/150‑1 (MQ=255)
                                   cGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtg                                                                                                                    <  1:771882/150‑1 (MQ=255)
                                    gTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgg                                                                                                                   <  1:3364216/150‑1 (MQ=255)
                                     tctcACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtggg                                                                                                                  <  1:3612516/150‑1 (MQ=255)
                                     tctcACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtggg                                                                                                                  <  1:5161646/150‑1 (MQ=255)
                                      ctcACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtggga                                                                                                                 <  1:2793848/150‑1 (MQ=255)
                                      ctcACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtggga                                                                                                                 >  2:550642/1‑150 (MQ=255)
                                      ctcACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtggga                                                                                                                 >  1:2431523/1‑150 (MQ=255)
                                       tcACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgggag                                                                                                                <  1:1785685/150‑1 (MQ=255)
                                       tcACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAAAAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgggag                                                                                                                >  1:1684278/1‑150 (MQ=255)
                                          cTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgggagtgg                                                                                                             <  1:319333/150‑1 (MQ=255)
                                          cTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgggagtgg                                                                                                             <  1:1212921/150‑1 (MQ=255)
                                            ggTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgggagtggga                                                                                                           <  1:3280097/150‑1 (MQ=255)
                                            ggTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCggggagtggga                                                                                                           >  1:3498716/1‑150 (MQ=255)
                                                cgcCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTg                                                                                                       <  1:1734932/150‑1 (MQ=255)
                                                   cacaCGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCTCATCACTCTGGCTAACTAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAg                                                                                                    >  1:210170/1‑150 (MQ=255)
                                                   cacaCGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAg                                                                                                    <  1:4184633/150‑1 (MQ=255)
                                                   cacaCGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAg                                                                                                    <  2:4489776/150‑1 (MQ=255)
                                                   cacaCGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAg                                                                                                    <  2:3817568/150‑1 (MQ=255)
                                                     cacGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGcc                                                                                                  <  2:4965808/150‑1 (MQ=255)
                                                        ggATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGaa                                                                                               >  1:2730731/1‑150 (MQ=255)
                                                           tGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAatat                                                                                            >  2:4838371/1‑150 (MQ=255)
                                                           tGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAatat                                                                                            >  2:3107278/1‑150 (MQ=255)
                                                            ggCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTACCGCGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATAtt                                                                                           <  2:133049/150‑1 (MQ=255)
                                                            ggCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATAtt                                                                                           >  2:5182648/1‑150 (MQ=255)
                                                            ggCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATAtt                                                                                           <  1:4793346/150‑1 (MQ=255)
                                                               aCGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATAttttt                                                                                        <  1:2434173/150‑1 (MQ=255)
                                                                   cgcTCTATTTTATCGCTCAGATCGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCt                                                                                    >  1:3251800/1‑150 (MQ=255)
                                                                    gcTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTg                                                                                   >  1:2929473/1‑150 (MQ=255)
                                                                      tctATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTggg                                                                                 >  2:1455271/1‑150 (MQ=255)
                                                                         aTTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATa                                                                              >  1:4047866/1‑150 (MQ=255)
                                                                         aTTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATa                                                                              <  2:1507560/150‑1 (MQ=255)
                                                                         aTTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATa                                                                              <  1:468599/150‑1 (MQ=255)
                                                                          ttttATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATaa                                                                             <  2:4547496/150‑1 (MQ=255)
                                                                          ttttATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATaa                                                                             >  2:1953881/1‑150 (MQ=255)
                                                                          ttttATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACAGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATaa                                                                             <  2:274314/150‑1 (MQ=255)
                                                                            ttATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGc                                                                           <  2:2101017/150‑1 (MQ=255)
                                                                             tATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCg                                                                          <  2:4249410/150‑1 (MQ=255)
                                                                                cGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGAt                                                                       <  1:3761311/150‑1 (MQ=255)
                                                                                cGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGAt                                                                       >  2:2331190/1‑150 (MQ=255)
                                                                                cGCTCAGATTGAAGACATTAACGAGATAAAACGCACCGAACAGTTCAATCAGCAACTGATGGAGCGCATAACTCTTGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGAt                                                                       <  1:268063/150‑1 (MQ=255)
                                                                                 gCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATg                                                                      >  2:3747358/1‑150 (MQ=255)
                                                                                  cTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGt                                                                     >  1:4341787/1‑150 (MQ=255)
                                                                                  cTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGt                                                                     >  1:5065059/1‑150 (MQ=255)
                                                                                    cAGATTGAAGACATTAACGAGCTAAAACTCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTTGTAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTc                                                                   <  2:4126081/150‑1 (MQ=255)
                                                                                    cAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCCGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTc                                                                   <  2:1675300/150‑1 (MQ=255)
                                                                                    cAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTc                                                                   >  2:654075/1‑150 (MQ=255)
                                                                                    cAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTc                                                                   <  1:4627043/150‑1 (MQ=255)
                                                                                     aGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCg                                                                  >  2:2284264/1‑150 (MQ=255)
                                                                                       aTTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAg                                                                >  1:150838/1‑150 (MQ=255)
                                                                                            aGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTa                                                           >  1:1302428/1‑150 (MQ=255)
                                                                                             gACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTAt                                                          <  1:4979900/150‑1 (MQ=255)
                                                                                              aCATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATg                                                         >  2:100008/1‑150 (MQ=255)
                                                                                               cATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGa                                                        >  2:4773005/1‑150 (MQ=255)
                                                                                                  tAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAAt                                                     <  1:1984237/150‑1 (MQ=255)
                                                                                                   aaCGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGACGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAAtt                                                    <  1:2539517/150‑1 (MQ=255)
                                                                                                   aaCGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAAtt                                                    <  1:4247292/150‑1 (MQ=255)
                                                                                                    aCGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATtc                                                   <  2:3572189/150‑1 (MQ=255)
                                                                                                      gAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATtcct                                                 >  1:1166915/1‑150 (MQ=255)
                                                                                                       aGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATtcctc                                                <  1:847664/150‑1 (MQ=255)
                                                                                                       aGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCACATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATtcctc                                                >  2:2918107/1‑150 (MQ=255)
                                                                                                           aaaaCGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCcaca                                            >  1:637923/1‑150 (MQ=255)
                                                                                                           aaaaCGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCcaca                                            <  1:669293/150‑1 (MQ=255)
                                                                                                            aaaCGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACat                                           >  1:1615950/1‑150 (MQ=255)
                                                                                                                   ccGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAAc                                    <  1:2316684/150‑1 (MQ=255)
                                                                                                                   ccGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAAc                                    <  1:3657007/150‑1 (MQ=255)
                                                                                                                   ccGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAAc                                    <  1:1415097/150‑1 (MQ=255)
                                                                                                                   ccGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAAc                                    <  2:3060792/150‑1 (MQ=255)
                                                                                                                        cAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAAc                               <  2:883792/150‑1 (MQ=255)
                                                                                                                        cAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAAc                               >  2:2642061/1‑150 (MQ=255)
                                                                                                                         aGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACt                              <  1:509153/150‑1 (MQ=255)
                                                                                                                          ggTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTg                             >  2:847807/1‑150 (MQ=255)
                                                                                                                          ggTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTg                             <  2:3498716/150‑1 (MQ=255)
                                                                                                                                 cAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGgtg                      <  1:862155/150‑1 (MQ=255)
                                                                                                                                      aCTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTa                 <  1:3913299/150‑1 (MQ=255)
                                                                                                                                          aTGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAg             >  2:4021333/1‑150 (MQ=255)
                                                                                                                                            ggAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAgtg           <  1:2826122/150‑1 (MQ=255)
                                                                                                                                            ggAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAgtg           <  1:2613243/150‑1 (MQ=255)
                                                                                                                                            ggAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAgtg           <  2:3924694/150‑1 (MQ=255)
                                                                                                                                             gAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAgtgc          <  2:3234722/150‑1 (MQ=255)
                                                                                                                                                  cATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAGTGCGtgct     >  2:3853776/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AACCGCAATGGCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAGTGCGTGCTGCC  >  CP009273/2137890‑2138188

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: