Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,552,629 C→A 45.7% D216Y (GAT→TAT)  ddpD ← D,D‑dipeptide permease system, ATP‑binding component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,552,6290CA45.7% 257.1 / inf 356D216Y (GAT→TAT) ddpDD,D‑dipeptide permease system, ATP‑binding component
Reads supporting (aligned to +/- strand):  ref base C (115/78);  new base A (77/85);  total (192/164)
Fisher's exact test for biased strand distribution p-value = 2.51e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 361 total aligned reads displayed.

TACTCCATGTTCTGGTGCGCATTGCAGCAAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGGAGAATGCCAGCGCAATCA  >  CP009273/1552480‑1552778
                                                                                                                                                     |                                                                                                                                                     
 aCTCCATGTTCTGGTGCGCATTGCAGCAAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCa                                                                                                                                                      <  1:682527/150‑1 (MQ=255)
  cTCCATGTTCTGGTGCGCATTGCAGCAAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAAt                                                                                                                                                     <  2:1917575/150‑1 (MQ=255)
     cATGTTCTGGTGCGCATTGCAGCAAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGc                                                                                                                                                  <  1:1333693/150‑1 (MQ=255)
         ttCTGGTGCGCATTGCAGCAAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCtgat                                                                                                                                              <  2:190678/150‑1 (MQ=255)
         ttCTGGTGCGCATTGCAGCAAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCtgat                                                                                                                                              >  2:1458931/1‑150 (MQ=255)
             ggTGCGCATTGCAGCAAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAAc                                                                                                                                          <  2:4326226/150‑1 (MQ=255)
              gTGCGCATTGCAGCAAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACa                                                                                                                                         <  2:3128247/150‑1 (MQ=255)
              gTGCGCATTGCAGCAAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACa                                                                                                                                         <  2:2956154/150‑1 (MQ=255)
              gTGCGCATTGCAGCAAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACa                                                                                                                                         <  1:2203902/150‑1 (MQ=255)
               tGCGCATTGCAGCAAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCTGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAg                                                                                                                                        >  1:4177138/1‑150 (MQ=255)
               tGCGCATTGCAGCAAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAg                                                                                                                                        <  2:1576812/150‑1 (MQ=255)
                    aTTGCAGCAAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCg                                                                                                                                   >  1:1808500/1‑150 (MQ=255)
                          gcaAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTc                                                                                                                             >  1:4200593/1‑150 (MQ=255)
                           caAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTcc                                                                                                                            <  1:328126/150‑1 (MQ=255)
                           caAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTcc                                                                                                                            >  1:1490671/1‑150 (MQ=255)
                            aaaCCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCa                                                                                                                           >  2:801373/1‑150 (MQ=255)
                             aaCCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCAc                                                                                                                          <  2:3534431/150‑1 (MQ=255)
                              aCCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACt                                                                                                                         <  1:1925216/150‑1 (MQ=255)
                              aCCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACt                                                                                                                         <  1:603523/150‑1 (MQ=255)
                                cAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTgg                                                                                                                       >  2:1987774/1‑150 (MQ=255)
                                        taCGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCt                                                                                                               <  1:1044091/150‑1 (MQ=255)
                                           ggATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACCGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTa                                                                                                            >  2:4075502/1‑150 (MQ=255)
                                           ggATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTa                                                                                                            >  2:3169738/1‑150 (MQ=255)
                                           ggATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTa                                                                                                            >  2:936784/1‑150 (MQ=255)
                                           ggATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTa                                                                                                            >  1:732409/1‑150 (MQ=255)
                                             aTGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATg                                                                                                          <  1:79813/150‑1 (MQ=255)
                                                ccGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGttt                                                                                                       <  1:2968861/150‑1 (MQ=255)
                                                ccGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGttt                                                                                                       <  1:364724/150‑1 (MQ=255)
                                                 cGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTa                                                                                                      >  1:2415228/1‑150 (MQ=255)
                                                 cGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTa                                                                                                      >  1:709458/1‑150 (MQ=255)
                                                 cGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTa                                                                                                      >  1:4465610/1‑150 (MQ=255)
                                                  ggggATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTaa                                                                                                     >  2:3077602/1‑150 (MQ=255)
                                                  ggggATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTaa                                                                                                     >  1:1535490/1‑150 (MQ=255)
                                                  ggggATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTaa                                                                                                     <  2:2178799/150‑1 (MQ=255)
                                                  ggggATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTaa                                                                                                     >  2:3438045/1‑150 (MQ=255)
                                                    ggATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTGTCTGTTTAAGc                                                                                                   >  2:2535096/1‑150 (MQ=255)
                                                      atgatgGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCaa                                                                                                 <  2:3654514/150‑1 (MQ=255)
                                                       tgatgGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCaaa                                                                                                >  1:3312632/1‑150 (MQ=255)
                                                        gatgGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAAc                                                                                               >  2:3282971/1‑150 (MQ=255)
                                                         atgGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACg                                                                                              >  2:3026294/1‑150 (MQ=255)
                                                          tgGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGc                                                                                             >  2:703019/1‑150 (MQ=255)
                                                           ggATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCa                                                                                            <  2:3301068/150‑1 (MQ=255)
                                                             aTAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTACGGCATACAACACGTAAACGCTATCGCACAGTTGCGACACCACGCCCATATCAAGGATGATGAACAGAACCGCAGTTCCAATGGCGCGGGATTTATGTTTAAGCAAACGCAGt                                                                                          >  2:3684975/1‑150 (MQ=255)
                                                               aaCGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTAc                                                                                        <  2:1066660/150‑1 (MQ=255)
                                                               aaCGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTAc                                                                                        >  2:1871790/1‑150 (MQ=255)
                                                                aCGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTAcc                                                                                       >  1:942616/1‑150 (MQ=255)
                                                                aCGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTAcc                                                                                       <  2:3842590/150‑1 (MQ=255)
                                                                 cGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCt                                                                                      >  1:3442059/1‑150 (MQ=255)
                                                                 cGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCt                                                                                      >  2:4366924/1‑150 (MQ=255)
                                                                 cGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCt                                                                                      >  2:2509868/1‑150 (MQ=255)
                                                                   tCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGc                                                                                    >  1:2540108/1‑150 (MQ=255)
                                                                          tCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCa                                                                             >  2:1790578/1‑150 (MQ=255)
                                                                           cACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCAc                                                                            >  1:3190920/1‑150 (MQ=255)
                                                                           cACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCAc                                                                            >  2:1308772/1‑150 (MQ=255)
                                                                             cgccgcTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCg                                                                          >  2:2596083/1‑150 (MQ=255)
                                                                              gccgcTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGt                                                                         >  1:3115178/1‑150 (MQ=255)
                                                                                   tttCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGt                                                                    >  2:3485313/1‑150 (MQ=255)
                                                                                    ttCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTc                                                                   >  1:2645383/1‑150 (MQ=255)
                                                                                    ttCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTc                                                                   <  2:783015/150‑1 (MQ=255)
                                                                                      cAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCa                                                                 >  1:544884/1‑150 (MQ=255)
                                                                                       aaTTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAg                                                                >  2:2278404/1‑150 (MQ=255)
                                                                                        aTTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAgc                                                               <  1:2766575/150‑1 (MQ=255)
                                                                                              cTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTa                                                         <  2:1522610/150‑1 (MQ=255)
                                                                                                 ccGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTc                                                      <  2:1683738/150‑1 (MQ=255)
                                                                                                    gCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGt                                                   <  2:4521019/150‑1 (MQ=255)
                                                                                                     cATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGtt                                                  <  1:1701340/150‑1 (MQ=255)
                                                                                                      aTACATCACGTAAACGCGATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTtc                                                 <  1:1190047/150‑1 (MQ=255)
                                                                                                        aCATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTtcgt                                               <  1:4740596/150‑1 (MQ=255)
                                                                                                               gTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCaat                                        >  2:2848464/1‑150 (MQ=255)
                                                                                                                tAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCaata                                       <  1:1436208/150‑1 (MQ=255)
                                                                                                                 aaaCGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCaataa                                      <  1:4648146/150‑1 (MQ=255)
                                                                                                                           gCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCgg                            <  1:1516544/150‑1 (MQ=255)
                                                                                                                             acaGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCt                          >  2:4031647/1‑150 (MQ=255)
                                                                                                                             acaGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCt                          >  2:2711675/1‑150 (MQ=255)
                                                                                                                             acaGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCt                          >  1:4290182/1‑150 (MQ=255)
                                                                                                                              caGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTc                         <  1:4072859/150‑1 (MQ=255)
                                                                                                                               aGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCg                        >  2:1741952/1‑150 (MQ=255)
                                                                                                                               aGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCg                        <  2:4188675/150‑1 (MQ=255)
                                                                                                                                  tGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAg                     >  1:122769/1‑150 (MQ=255)
                                                                                                                                    cGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGga                   >  2:4727846/1‑150 (MQ=255)
                                                                                                                                    cGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGga                   >  2:2419593/1‑150 (MQ=255)
                                                                                                                                     gACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGgag                  >  1:4372268/1‑150 (MQ=255)
                                                                                                                                     gACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGgag                  >  1:2993647/1‑150 (MQ=255)
                                                                                                                                      acacCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGgaga                 >  2:476374/1‑150 (MQ=255)
                                                                                                                                      acacCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGgaga                 <  1:3566822/150‑1 (MQ=255)
                                                                                                                                        accacGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGGAGAAt               <  2:3673598/150‑1 (MQ=255)
                                                                                                                                         ccacGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGGAGAATg              >  1:3812322/1‑150 (MQ=255)
                                                                                                                                         ccacGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGGAGAATg              <  2:2120852/150‑1 (MQ=255)
                                                                                                                                         ccacGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGGAGAATg              <  1:3754770/150‑1 (MQ=255)
                                                                                                                                         ccacGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGGAGAATg              >  2:2398189/1‑150 (MQ=255)
                                                                                                                                         ccacGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGGAGAATg              <  1:3502262/150‑1 (MQ=255)
                                                                                                                                           acGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGGAGAATGcc            <  2:54493/150‑1 (MQ=255)
                                                                                                                                            cGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGGAGAATGCCa           <  2:1400316/150‑1 (MQ=255)
                                                                                                                                             ggCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGGAGAATGCCAg          >  2:4226991/1‑150 (MQ=255)
                                                                                                                                               ccATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGGAGAATGCCAgcg        <  2:873420/150‑1 (MQ=255)
                                                                                                                                                cATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGGAGAATGCCAgcgc       >  1:3456990/1‑150 (MQ=255)
                                                                                                                                                 atctCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGGAGAATGCCAGCGCa      <  1:400306/147‑1 (MQ=255)
                                                                                                                                                 atatCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGGAGAATGCCAGCGCa      >  1:1144897/1‑150 (MQ=255)
                                                                                                                                                   atCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGGAGAATGCCAGCGCAAt    <  2:1199305/150‑1 (MQ=255)
                                                                                                                                                     cATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGGAGAATGCCAGCGCAATCa  >  1:3982619/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TACTCCATGTTCTGGTGCGCATTGCAGCAAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGGAGAATGCCAGCGCAATCA  >  CP009273/1552480‑1552778

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: