Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,619,105 T→A 19.1% A353A (GCA→GCT yhhJ ← inner membrane putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,619,1050TA19.1% 679.7 / 153.0 266A353A (GCA→GCTyhhJinner membrane putative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (117/98);  new base A (24/27);  total (141/125)
Fisher's exact test for biased strand distribution p-value = 3.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 274 total aligned reads displayed.

TTATCCGGCCTAAAAACAACTAAAATTCAATAAATTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGCATAATGTCCTGCACCATCTGCGGCATACTTTCGCGCGGCGTGGAACC  >  CP009273/3618957‑3619254
                                                                                                                                                    |                                                                                                                                                     
ttATCCGGCCTAAAAACAACTAAAATTCAATAAATTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTg                                                                                                                                                      <  2:3445422/150‑1 (MQ=255)
  aTCCGGCCTAAAAACAACTAAAATTCAATAAATTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGcc                                                                                                                                                    <  2:1206950/150‑1 (MQ=255)
     cGGCCTAAAAACAACTAAAATTCAATAAATTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATc                                                                                                                                                 >  2:16879/1‑150 (MQ=255)
     cGGCCTAAAAACAACTAAAATTCAATAAATTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTAGTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATc                                                                                                                                                 >  2:15629/1‑150 (MQ=255)
      ggCCTAAAAACAACTAAAATTCAATAAATTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCa                                                                                                                                                >  2:4174220/1‑150 (MQ=255)
       gCCTAAAAACAACTAAAATTCAATAAATTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAg                                                                                                                                               >  2:4357499/1‑150 (MQ=255)
         cTAAAAACAACTAAAATTCAATAAATTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCg                                                                                                                                             >  2:228710/1‑150 (MQ=255)
          tAAAAACAACTAAAATTCAATAAATTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGt                                                                                                                                            >  1:3717133/1‑150 (MQ=255)
           aaaaaCAACTAAAATTCGATAAATTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATAGCCATCAGCGTc                                                                                                                                           >  2:764078/1‑150 (MQ=255)
           aaaaaCAACTAAAATTCAATAAATTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTCCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTc                                                                                                                                           <  1:613811/150‑1 (MQ=255)
            aaaaCAACTAAAATTCAATAAATTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCa                                                                                                                                          <  2:3852348/150‑1 (MQ=255)
               acaacTAAAATTCAATAAATTGCACCGATGCGTAGTCCCTATAAGATTACGCATCGGGCAATTGTATTAACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATAGCCATCAGCGTCAGaa                                                                                                                                       >  1:2669656/1‑150 (MQ=255)
                caacTAAAATTCAATAAATTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATAGCCATCAGCGTCAGaaa                                                                                                                                      >  1:4380209/1‑150 (MQ=255)
                 aacTAAAATTCAATAAATTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAAc                                                                                                                                     >  2:3587241/1‑150 (MQ=255)
                 aacTAAAATTCAATAAATTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATAGCCATCAGCGTCAGAAAc                                                                                                                                     >  1:4254495/1‑150 (MQ=255)
                    tAAAATTCAATAAATTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGc                                                                                                                                  >  1:1468127/1‑150 (MQ=255)
                    tAAAATTCAATAAATTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGc                                                                                                                                  >  2:4199498/1‑150 (MQ=255)
                     aaaaTTCAATAAATTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCg                                                                                                                                 >  1:991803/1‑150 (MQ=255)
                     aaaaTTCAATAAATTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATAGCCATCAGCGTCAGAAACTGCg                                                                                                                                 >  2:3140918/1‑150 (MQ=255)
                       aaTTCAATAAATTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGc                                                                                                                               >  2:718806/1‑150 (MQ=255)
                        aTTCAATAAATTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATAGCCATCAGCGTCAGAAACTGCGGcc                                                                                                                              <  2:2853479/150‑1 (MQ=255)
                         ttCAATAAATTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATAGCCATCAGCGTCAGAAACTGCGGCCa                                                                                                                             <  2:4133273/150‑1 (MQ=255)
                          tCAATAAATTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAg                                                                                                                            >  2:684347/1‑150 (MQ=255)
                          tCAATAAATTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATAGCCATCAGCGTCAGAAACTGCGGCCAg                                                                                                                            >  1:4296486/1‑150 (MQ=255)
                           cAATAAATTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGa                                                                                                                           >  1:1075120/1‑150 (MQ=255)
                           cAATAAATTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGa                                                                                                                           >  2:1019770/1‑150 (MQ=255)
                            aaTAAATTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAACAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGAc                                                                                                                          >  1:2629862/1‑150 (MQ=255)
                            aaTAAATTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCgg                                                                                                                                <  2:2136369/144‑1 (MQ=255)
                             aTAAATTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACg                                                                                                                         <  1:1401024/150‑1 (MQ=255)
                               aaaTTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGAt                                                                                                                       <  2:3337703/150‑1 (MQ=255)
                               aaaTTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGAt                                                                                                                       <  2:4193852/150‑1 (MQ=255)
                                aaTTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGAtt                                                                                                                      <  2:3261821/150‑1 (MQ=255)
                                aaTTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGAtt                                                                                                                      <  2:1022677/150‑1 (MQ=255)
                                 aTTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATAGCCATCAGCGTCAGAAACTGCGGCCAGACGAttt                                                                                                                     >  1:1905885/1‑150 (MQ=255)
                                   tGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCg                                                                                                                   <  1:3431830/150‑1 (MQ=255)
                                      aCCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAAt                                                                                                                <  1:3565377/150‑1 (MQ=255)
                                         gATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCg                                                                                                             >  1:1295073/1‑150 (MQ=255)
                                              gTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCAcc                                                                                                        >  2:1713796/1‑150 (MQ=255)
                                              gTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGAGTTCGAATCCGGCAcc                                                                                                        >  1:81989/1‑150 (MQ=255)
                                                   cccTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGt                                                                                                   <  2:2454378/150‑1 (MQ=255)
                                                    ccTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTa                                                                                                  <  1:2483538/150‑1 (MQ=255)
                                                     cTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTag                                                                                                 <  1:895421/150‑1 (MQ=255)
                                                       ataAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTagag                                                                                               <  2:921913/150‑1 (MQ=255)
                                                              tACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATAGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGcc                                                                                        <  1:2929737/150‑1 (MQ=255)
                                                               aCGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCt                                                                                       <  1:1346203/150‑1 (MQ=255)
                                                               aCGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCt                                                                                       <  1:3630184/150‑1 (MQ=255)
                                                                cGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTg                                                                                      >  2:2416023/1‑150 (MQ=255)
                                                                  cATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATAGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTgcg                                                                                    >  2:1321363/1‑150 (MQ=255)
                                                                      gggCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAg                                                                                >  1:10297/1‑150 (MQ=255)
                                                                        gCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGc                                                                              <  2:2785145/150‑1 (MQ=255)
                                                                              gTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACaa                                                                        >  2:841325/1‑150 (MQ=255)
                                                                                aTTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAg                                                                      >  2:3748119/1‑150 (MQ=255)
                                                                                  ttACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAgtg                                                                    <  2:4171270/150‑1 (MQ=255)
                                                                                     cGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAgtgtgt                                                                 <  2:404020/150‑1 (MQ=255)
                                                                                      gCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTgtc                                                                <  2:1526328/150‑1 (MQ=255)
                                                                                       ccATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTgtcg                                                               <  2:2926814/150‑1 (MQ=255)
                                                                                       ccATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTgtcg                                                               <  1:4510210/150‑1 (MQ=255)
                                                                                         cttGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATAGCCATCAGCGTCAGAAACTGCGGCCAGACGAGTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTgtcgtc                                                             <  2:3926018/149‑1 (MQ=255)
                                                                                         aTTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTgtcgtc                                                             <  2:2558325/150‑1 (MQ=255)
                                                                                         aTTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATAGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTgtcgtc                                                             <  2:3926076/150‑1 (MQ=255)
                                                                                           tGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCgg                                                           <  2:4042514/150‑1 (MQ=255)
                                                                                             tCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCa                                                         <  2:510645/150‑1 (MQ=255)
                                                                                                 aaTCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGt                                                     >  2:3700623/1‑150 (MQ=255)
                                                                                                  aTCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTc                                                    >  2:1944553/1‑150 (MQ=255)
                                                                                                  aTCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTc                                                    >  2:3364812/1‑150 (MQ=255)
                                                                                                  aTCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCa                                                   >  2:443606/1‑150 (MQ=255)
                                                                                                   tCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACACAGTGTGTCGTCGGCATGGTCa                                                   <  2:4274688/150‑1 (MQ=255)
                                                                                                    cGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAg                                                  >  1:1622041/1‑150 (MQ=255)
                                                                                                    cGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAg                                                  >  2:4158780/1‑150 (MQ=255)
                                                                                                     gTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGGGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGc                                                 >  2:3247216/1‑150 (MQ=255)
                                                                                                        ttCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGCATa                                              <  2:1808023/150‑1 (MQ=255)
                                                                                                          acTGAATCGCAGCAGCGTTATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGTAACTGTGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGATGCTATCAAAGTGTGTCTTCGGCATGGTCAGCTTAAt                                            <  2:3120824/149‑1 (MQ=255)
                                                                                                            tGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGCATAATGt                                          <  2:2457488/150‑1 (MQ=255)
                                                                                                            tGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGCATAATGt                                          >  1:2814702/1‑150 (MQ=255)
                                                                                                              aaTCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGCATAATGTcc                                        >  1:2440260/1‑150 (MQ=255)
                                                                                                              aaTCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGCATAATGTcc                                        >  1:652097/1‑150 (MQ=255)
                                                                                                                tCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGCATAATGTCCTg                                      <  1:3010161/150‑1 (MQ=255)
                                                                                                                tCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGCATAATGTCCTg                                      <  1:3910947/150‑1 (MQ=255)
                                                                                                                 cGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGCATAATGTCCTGc                                     >  1:1520278/1‑150 (MQ=255)
                                                                                                                   cagcagCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGCATAATGTCCTGCAc                                   >  1:3293595/1‑150 (MQ=255)
                                                                                                                   cagcagCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGCATAATGTCCTGCAc                                   >  1:1014538/1‑150 (MQ=255)
                                                                                                                      cagcGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGCATAATGTCCTGCACCAt                                <  1:4030047/150‑1 (MQ=255)
                                                                                                                       agcGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGCATAATGTCCTGCACCATc                               <  1:3123847/150‑1 (MQ=255)
                                                                                                                        gcgcAATGGTGAAAAATGCGCCGCCAAGTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGCATAATGTCCTGCACCATCt                              <  2:4450/150‑1 (MQ=255)
                                                                                                                         cgcAATGGTGAAAAATGCGCCGCCAATAGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGCATAATGTCCTGCACCATCTg                             >  1:1711996/1‑150 (MQ=255)
                                                                                                                          gcAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGCATAATGTCCTGCACCATCTGc                            >  1:2020561/1‑150 (MQ=255)
                                                                                                                          gcAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGCATAATGTCCTGCACCATCTGc                            >  1:2019842/1‑150 (MQ=255)
                                                                                                                           cAATGGTGAAAAATGCGCCGCCAATAGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGCATAATGTCCTGCACCATCTGCg                           <  1:225480/150‑1 (MQ=255)
                                                                                                                                      aaTGCGCCGCCAATAGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGGAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGCATAATGTCCTGCACCATCTGCGGCATACTTTcg                >  2:1610276/1‑150 (MQ=255)
                                                                                                                                       aTGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGCATAATGTCCTGCACCATCTGCGGCATACTTTcgc               >  1:4154470/1‑150 (MQ=255)
                                                                                                                                          cgccgcCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGCATAATGTCCTGCACCATCTGCGGCATACTTTcgcgcg            <  1:3397734/150‑1 (MQ=255)
                                                                                                                                              gccAATAGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGCATAATGTCCTGCACCATCTGCGGCATACTTTCGCGCGGCGt        <  1:4341147/150‑1 (MQ=255)
                                                                                                                                               ccAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGCATAATGTCCTGCACCATCTGCGGCATACTTTCGCGCGGCGTg       >  1:4122621/1‑150 (MQ=255)
                                                                                                                                                cAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGCATAATGTCCTGCACCATCTGCGGCATACTTTCGCGCGGCGTgg      >  1:4127193/1‑150 (MQ=255)
                                                                                                                                                 aaTTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGCATAATGTCCTGCACCATCTGCGGCATACTTTCGCGCGGCGTGGa     >  2:2011687/1‑150 (MQ=255)
                                                                                                                                                 aaTTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGCATAATGTCCTGCACCATCTGCGGCATACTTTCGCGCGGCGTGGa     <  2:2094781/150‑1 (MQ=255)
                                                                                                                                                    tGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGCATAATGTCCTGCACCATCTGCGGCATACTTTCGCGCGGCGTGGAacc  >  2:2736428/1‑150 (MQ=255)
                                                                                                                                                    tGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGCATAATGTCCTGCACCATCTGCGGCATACTTTCGCGCGGCGTGGAacc  >  2:1096872/1‑150 (MQ=255)
                                                                                                                                                    tGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGCATAATGTCCTGCACCATCTGCGGCATACTTTCGCGCGGCGTGGAacc  >  1:792671/1‑150 (MQ=255)
                                                                                                                                                    tGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGCATAATGTCCTGCACCATCTGCGGCATACTTTCGCGCGGCGTGGAacc  >  2:3033323/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTATCCGGCCTAAAAACAACTAAAATTCAATAAATTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGAAAAATGCGCCGCCAATTGCCATCAGCGTCAGAAACTGCGGCCAGACGATTTCGAATCCGGCACCCCGGTAGAGGATGGCCTGCGCGAGGCTAACAAAGTGTGTCGTCGGCATGGTCAGCATAATGTCCTGCACCATCTGCGGCATACTTTCGCGCGGCGTGGAACC  >  CP009273/3618957‑3619254

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: