Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,552,629 C→A 94.5% 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,6290CA94.5% 899.6 / 40.0 293D216Y (GAT→TAT) ddpDD,D‑dipeptide permease system, ATP‑binding component
Reads supporting (aligned to +/- strand):  ref base C (7/9);  new base A (161/116);  total (168/125)
Fisher's exact test for biased strand distribution p-value = 3.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 293 total aligned reads displayed.

CTCCATGTTCTGGTGCGCATTGCAGCAAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGGAGAATGCCAGCGCA  >  CP009273/1552482‑1552774
                                                                                                                                                   |                                                                                                                                                 
cTCCATGTTCTGGTGCGCATTGCAGCAAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAAt                                                                                                                                                 >  1:3770901/1‑150 (MQ=255)
 tCCATGTTCTGGTGCGCATTGCAGCAAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATg                                                                                                                                                >  1:4035270/1‑150 (MQ=255)
   cATGTTCTGGTGCGCATTGCAGCAAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGc                                                                                                                                              <  2:237881/150‑1 (MQ=255)
    aTGTTCTGGTGCGCATTGCAGCAAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCt                                                                                                                                             >  2:1183155/1‑150 (MQ=255)
        tCTGGTGCGCATTGCAGCAAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCtgatg                                                                                                                                         >  1:941513/1‑150 (MQ=255)
          tGGTGCGCATTGCAGCAAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGaa                                                                                                                                       <  2:342369/150‑1 (MQ=255)
             tGCGCATTGCAGCAAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAg                                                                                                                                    >  1:1947048/1‑150 (MQ=255)
             tGCGCATTGCAGCAAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAg                                                                                                                                    >  2:442767/1‑150 (MQ=255)
                gcATTGCAGCAAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAAGGGCTGATGAACAGAAc                                                                                                                                 <  2:2534967/150‑1 (MQ=255)
                 cATTGCAGCAAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAAcc                                                                                                                                >  1:3325190/1‑150 (MQ=255)
                    tGCAGCAAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCa                                                                                                                             >  1:1565759/1‑150 (MQ=255)
                      cagcaAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGt                                                                                                                           <  1:3928937/150‑1 (MQ=255)
                      cagcaAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGt                                                                                                                           >  2:1327208/1‑150 (MQ=255)
                       agcaAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGtt                                                                                                                          >  2:2047765/1‑150 (MQ=255)
                         caATCCAATGGTATACGGATGCCGGGGATTATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCTCAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTcc                                                                                                                        <  1:2703627/150‑1 (MQ=255)
                         caAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCAGATGAACAGAACCGCAGTTcc                                                                                                                        >  2:2342837/1‑150 (MQ=255)
                           aaCCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCAc                                                                                                                      >  1:507389/1‑150 (MQ=255)
                           aaCCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCAc                                                                                                                      <  1:128082/150‑1 (MQ=255)
                           aaCCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCAc                                                                                                                      >  1:1500885/1‑150 (MQ=255)
                             ccAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTg                                                                                                                    <  2:225584/150‑1 (MQ=255)
                              cAATGGTATACGGATGTCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTgg                                                                                                                   <  2:4064966/150‑1 (MQ=255)
                              cAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTgg                                                                                                                   >  2:3110224/1‑150 (MQ=255)
                               aaTGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGgc                                                                                                                  <  2:154931/150‑1 (MQ=255)
                                 tGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGgcgc                                                                                                                >  2:3825662/1‑150 (MQ=255)
                                   gTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCgg                                                                                                              >  1:2988630/1‑150 (MQ=255)
                                        cGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCttt                                                                                                         <  1:876873/150‑1 (MQ=255)
                                         ggATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTa                                                                                                        >  2:1910266/1‑150 (MQ=255)
                                          gATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTAt                                                                                                       <  1:3426105/150‑1 (MQ=255)
                                          gATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTAt                                                                                                       >  1:4171196/1‑150 (MQ=255)
                                          gATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTAt                                                                                                       >  1:1986044/1‑150 (MQ=255)
                                            tGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGt                                                                                                     >  2:1121513/1‑150 (MQ=255)
                                               cGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACAACACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAAAGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGttaa                                                                                                  >  1:249451/1‑148 (MQ=255)
                                                ggggATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTaa                                                                                                 >  2:758267/1‑150 (MQ=255)
                                                     tgatgGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCaaa                                                                                            <  2:2936391/150‑1 (MQ=255)
                                                      gatgGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAAc                                                                                           <  1:2190497/150‑1 (MQ=255)
                                                       atgGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACg                                                                                          <  2:3239269/150‑1 (MQ=255)
                                                        tgGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGc                                                                                         <  2:4346414/150‑1 (MQ=255)
                                                           aTAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGt                                                                                      >  2:3890028/1‑150 (MQ=255)
                                                             aaCGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTAc                                                                                    >  2:958735/1‑150 (MQ=255)
                                                             aaCGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTAc                                                                                    <  1:2016572/150‑1 (MQ=255)
                                                               cGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCt                                                                                  <  1:2588826/150‑1 (MQ=255)
                                                                gTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACTGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTg                                                                                 >  2:128945/1‑150 (MQ=255)
                                                                gTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTg                                                                                 <  1:1345358/150‑1 (MQ=255)
                                                                  cTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCa                                                                               <  2:3529902/150‑1 (MQ=255)
                                                                    gCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAAc                                                                             <  2:671696/150‑1 (MQ=255)
                                                                     ccGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACt                                                                            <  2:2096182/150‑1 (MQ=255)
                                                                      cGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTg                                                                           <  1:113084/150‑1 (MQ=255)
                                                                           cgccgcTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCg                                                                      <  1:367329/150‑1 (MQ=255)
                                                                            gccgcTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGt                                                                     <  2:698699/150‑1 (MQ=255)
                                                                              cgcTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGa                                                                   >  1:1385218/1‑150 (MQ=255)
                                                                               gcTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGAc                                                                  <  2:2788761/150‑1 (MQ=255)
                                                                                 tttCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCACAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGt                                                                >  2:1831575/1‑150 (MQ=255)
                                                                                  ttCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTc                                                               >  2:1236212/1‑150 (MQ=255)
                                                                                    cAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCa                                                             <  1:3630134/150‑1 (MQ=255)
                                                                                    cAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCa                                                             >  1:4094238/1‑150 (MQ=255)
                                                                                     aaTTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAg                                                            >  1:4392000/1‑150 (MQ=255)
                                                                                      aTTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAgc                                                           >  1:4629872/1‑150 (MQ=255)
                                                                                       ttACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAgcg                                                          <  2:3873306/150‑1 (MQ=255)
                                                                                        tACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAgcgc                                                         >  1:4437282/1‑150 (MQ=255)
                                                                                        tACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAgcgc                                                         >  1:4266222/1‑150 (MQ=255)
                                                                                          cGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCg                                                       >  1:2385328/1‑150 (MQ=255)
                                                                                          cGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCg                                                       >  1:1210319/1‑150 (MQ=255)
                                                                                               ccagCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTc                                                  <  1:2262316/147‑1 (MQ=255)
                                                                                                cGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCg                                                 >  1:3388787/1‑150 (MQ=255)
                                                                                                cGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCg                                                 <  1:1632761/150‑1 (MQ=255)
                                                                                                cGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCg                                                 <  1:3956579/150‑1 (MQ=255)
                                                                                                 ggCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCgg                                                >  1:2871615/1‑150 (MQ=255)
                                                                                                 ggCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCgg                                                >  2:4633853/1‑150 (MQ=255)
                                                                                                   cATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGAAGTCCAGCGCCGTAGTCGGtt                                              >  2:1606763/1‑150 (MQ=255)
                                                                                                    aTACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTtc                                             <  1:999645/150‑1 (MQ=255)
                                                                                                    aTACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTtc                                             <  1:79877/150‑1 (MQ=255)
                                                                                                    aTACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTtc                                             <  1:4373964/150‑1 (MQ=255)
                                                                                                       cATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTtcgtc                                          <  2:2306576/150‑1 (MQ=255)
                                                                                                        aTCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTtcgtcg                                         <  2:937741/150‑1 (MQ=255)
                                                                                                        aTCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTtcgtcg                                         >  1:647714/1‑150 (MQ=255)
                                                                                                          cACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGc                                       >  1:1689258/1‑150 (MQ=255)
                                                                                                           aCGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCa                                      <  1:4083006/150‑1 (MQ=255)
                                                                                                           aCGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCa                                      >  2:1991041/1‑150 (MQ=255)
                                                                                                                aaCGCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACAGCACCGTGACGTCCAGCGCCGTAGTCGGATCGTCGGCaaaaat                                 >  1:3822208/1‑150 (MQ=255)
                                                                                                                   gCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCaa                              >  2:3789993/1‑150 (MQ=255)
                                                                                                                   gCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCaa                              <  2:2482592/150‑1 (MQ=255)
                                                                                                                   gCTATCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCaa                              >  1:319779/1‑150 (MQ=255)
                                                                                                                      aTCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTg                           <  2:867672/150‑1 (MQ=255)
                                                                                                                       tCGCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGc                          >  2:448356/1‑150 (MQ=255)
                                                                                                                         gCACAGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCTGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCgg                        >  1:2767669/1‑150 (MQ=255)
                                                                                                                          cacaGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGc                       >  1:1401584/1‑150 (MQ=255)
                                                                                                                           acaGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCt                      <  2:210820/150‑1 (MQ=255)
                                                                                                                           acaGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCt                      <  1:1254072/150‑1 (MQ=255)
                                                                                                                             aGTTGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCg                    >  1:4239970/1‑150 (MQ=255)
                                                                                                                               ttGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCa                  <  2:2170293/150‑1 (MQ=255)
                                                                                                                               ttGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCa                  <  2:1486700/150‑1 (MQ=255)
                                                                                                                                tGCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAg                 >  2:2581413/1‑150 (MQ=255)
                                                                                                                                 gCGACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAgg                <  1:1065160/150‑1 (MQ=255)
                                                                                                                                   gACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGgag              >  2:3755853/1‑150 (MQ=255)
                                                                                                                                   gACACCACGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGgag              <  2:1717426/150‑1 (MQ=255)
                                                                                                                                     caccacGGCCATATAATGGCTGATGACCAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGGAGaa            <  1:1744639/150‑1 (MQ=255)
                                                                                                                                     caccacGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGGAGaa            >  2:3623441/1‑150 (MQ=255)
                                                                                                                                       ccacGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGGAGAATg          >  1:2513260/1‑150 (MQ=255)
                                                                                                                                         acGGCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGGAGAATGcc        >  2:2432263/1‑150 (MQ=255)
                                                                                                                                            gCCATATAATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGGAGAATGCCAgc     >  1:1704818/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CTCCATGTTCTGGTGCGCATTGCAGCAAACCAATGGTATACGGATGCCGGGGATGATGGATAACGTCTGCCGTCACGCCGCTTTCAATTACGCTTCCGGCATACATCACGTAAACGCTATCGCACAGTTGCGACACCACGGCCATATCATGGCTGATGAACAGAACCGCAGTTCCACTGGCGCGGGCTTTATGTTTAAGCAAACGCAGTACCTGCAACTGCACCGTGACGTCCAGCGCCGTAGTCGGTTCGTCGGCAATAATCAATTGCGGCTCGCAGGAGAATGCCAGCGCA  >  CP009273/1552482‑1552774

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: