Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 323,222 G→A 6.5% intergenic (+34/+18) lytM → / ← USA300HOU_0290 M23 family lysostaphin/ABC superfamily ATP binding cassette transporter, ABC protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730323,2220GA 6.5% 288.1 / 4.2 95intergenic (+34/+18)lytM/USA300HOU_0290M23 family lysostaphin/ABC superfamily ATP binding cassette transporter, ABC protein
Reads supporting (aligned to +/- strand):  ref base G (44/44);  new base A (4/3);  total (48/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTACGGGTAATTCAACAGCGCCTCACGTACACTTCCAACGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAATTAACTTGATATCATCTGCAAGTGCTTCAACTTCATCTAAATGATGGGTAGTTAAT  >  CP000730/323084‑323360
                                                                                                                                           |                                                                                                                                          
agtaCGGGTAATTCAACAGCGCCTCACGTACACTTCCAACGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCCAAGTAGAGAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCt                                                                                                                                          <  1:248084/141‑1 (MQ=255)
 gtaCGGGTAATTCAACAGCGCCTCACGTACACTTCCAACGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCtt                                                                                                                                         <  2:211033/141‑1 (MQ=255)
   aCGGGTAATTCAACAGCGCCTCACGTACACTTCCAACGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCCTACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTgg                                                                                                                                       <  1:959023/141‑1 (MQ=255)
     gggTAATTCAACAGCGCCTCACGTACACTTCCAACGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGa                                                                                                                                     <  1:256779/141‑1 (MQ=255)
          aTTCAACAGCGCCTCACGTACACTTCCAACGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGAttttt                                                                                                                                <  1:2327/141‑1 (MQ=255)
          aTTCAACAGCGCCTCACGTACACTTCCAACGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGAttttt                                                                                                                                <  2:385308/141‑1 (MQ=255)
           ttCAACAGCGCCTCACGTACACTTCCAACGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTTATACTTGCAAAGTAGATAATACAGAATATCCCTAGCTGCG‑TATGATATCGCAACTTGGGAttttct                                                                                                                               <  2:179678/141‑3 (MQ=37)
            tCAACAGCGCCTCACGTACACTTCCAACGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTc                                                                                                                              >  1:53744/1‑141 (MQ=255)
              aaCAGCGCCTCACGTACACTTCCAACGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGt                                                                                                                            >  2:104888/1‑141 (MQ=255)
               aCAGCGCCTCACGTACACTTCCAACGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGGTAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGtt                                                                                                                           <  1:233094/141‑1 (MQ=255)
               aCAGCGCCTCACGTACACTTCCAACGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGtt                                                                                                                           <  1:347429/141‑1 (MQ=255)
               aCAGCGCCTCACGTACACTTCCAACGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGtt                                                                                                                           >  2:144297/1‑141 (MQ=255)
               aCAGCGCCTCACGTACACTTCCAACGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGtt                                                                                                                           >  1:27924/1‑141 (MQ=255)
                cAGCGCCTCACGTACACTTCCAACGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGttt                                                                                                                          <  2:992464/141‑1 (MQ=255)
                  gcgcCTCACGTACACTTCCAACGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTta                                                                                                                        <  2:622308/141‑1 (MQ=255)
                      cTCACGTACACTTCCAACGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGAGCCAACGTCATACTTGCAAAGTAGATAATACAGTAAATCCCAAGTTACAATATCATA‑CGCAACTTGGGATTTTTTCGTTTtaataa                                                                                                                    <  2:334576/141‑1 (MQ=255)
                      cTCACGTACACTTCCAACGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTtaataa                                                                                                                    <  2:76142/141‑1 (MQ=255)
                       tCACGTACACTTCCAACGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTaataat                                                                                                                   >  2:305028/1‑141 (MQ=255)
                         aCGTACACTTCCAACGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAg                                                                                                                 >  2:96263/1‑141 (MQ=255)
                           gTACACTTCCAACGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGt                                                                                                               >  1:875432/1‑141 (MQ=255)
                           gTACACTTCCAACGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGt                                                                                                               >  1:928979/1‑141 (MQ=255)
                             acacTTCCAACGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGtat                                                                                                             <  2:637229/141‑1 (MQ=255)
                              cacTTCCAACGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCATAGTTGCAAAGTAGATAATACAGAAAACCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGtata                                                                                                            <  1:246374/141‑1 (MQ=255)
                                cTTCCAACGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAg                                                                                                          <  2:955170/141‑1 (MQ=255)
                                 ttCCAACGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGt                                                                                                         <  1:530425/141‑1 (MQ=255)
                                  tCCAACGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTc                                                                                                        >  2:445890/1‑141 (MQ=255)
                                    cAACGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCa                                                                                                      <  1:163515/141‑1 (MQ=255)
                                    cAACGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCa                                                                                                      <  1:769326/141‑1 (MQ=255)
                                      aCGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCAtt                                                                                                    >  2:195395/1‑141 (MQ=255)
                                       cGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCAttg                                                                                                   <  2:828566/141‑1 (MQ=255)
                                          aTGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCAttgttt                                                                                                >  1:412178/1‑141 (MQ=255)
                                             tCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTc                                                                                             >  1:2150/1‑141 (MQ=255)
                                                  tggCATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTa                                                                                        <  1:324421/141‑1 (MQ=255)
                                                     cATCGGTAATCAATATGCAGTAGACCCACCGTCATACTTGCAAAGTAGATAATACAGAAAAGCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTaaaa                                                                                     <  2:547146/141‑1 (MQ=255)
                                                     cATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTaaaa                                                                                     >  1:750458/1‑141 (MQ=255)
                                                     cATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTaaaa                                                                                     >  1:817942/1‑141 (MQ=255)
                                                         ggTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATg                                                                                 >  1:313126/1‑141 (MQ=255)
                                                         gattaTCAATATGCGGTTGTCACAACGTCATACTTGCAATGTTGATAATATTGAAAATCCCAAATAGCGATATCATA‑CGCAGCTTGGCATTTTTTCGTTTTAATAATAGGTATACGTCCATTTTTTGTCCTCTAAAAAATg                                                                                 <  2:977531/137‑1 (MQ=255)
                                                          gTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGt                                                                                >  2:739224/1‑141 (MQ=255)
                                                          gTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGt                                                                                >  2:45260/1‑141 (MQ=255)
                                                            aaTCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGCTGCGAAATAATA‑CGCAACTTGGGATTTTTTTTTTTTTTTTATAGGTATAATTCCATTGTTTTTGTTCTTAAACATAttt                                                                              >  2:995612/1‑140 (MQ=255)
                                                             aTCAATATGCAGTAGCCCGGCCGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGa                                                                             <  1:283614/141‑1 (MQ=255)
                                                             aTCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCGTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGa                                                                             <  2:963654/141‑1 (MQ=255)
                                                              tCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGAt                                                                            <  1:581568/141‑1 (MQ=255)
                                                              tCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGAt                                                                            <  1:641760/141‑1 (MQ=255)
                                                               cAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATa                                                                           >  2:840659/1‑141 (MQ=255)
                                                               cAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATa                                                                           >  1:459485/1‑141 (MQ=255)
                                                               cAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATa                                                                           >  2:242038/1‑141 (MQ=255)
                                                                   atGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATaaaaa                                                                       <  1:129383/141‑1 (MQ=255)
                                                                      cAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATaaaccgca                                                                    >  1:817561/1‑136 (MQ=255)
                                                                      cAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGa                                                                    >  2:735954/1‑141 (MQ=255)
                                                                        ggaGACCCAACGTCATACGTGCAAAGTAGATAATACAGAAAAGCCCAAGTTGCGATATCATA‑CGCTGCTTGGGATTTTTGCGTTTTAATAATATTTATAAGTCCATTGTTTGTCCTCTAAAACATGGTGATAAAAAGGATc                                                                  <  2:974411/139‑1 (MQ=255)
                                                                        gTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAGATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATc                                                                  >  2:895076/1‑141 (MQ=255)
                                                                        gTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATc                                                                  >  1:661884/1‑141 (MQ=255)
                                                                           gACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTCTTAATAATAGGTATAAGTCCATTGTTTCTCCTCTAAAACATCTTGATAAAAAGGATCTTg                                                               >  1:937443/1‑141 (MQ=255)
                                                                              ccAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGcc                                                            >  1:90017/1‑141 (MQ=255)
                                                                              ccAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGcc                                                            >  1:89970/1‑141 (MQ=255)
                                                                              ccAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGcc                                                            <  2:866761/141‑1 (MQ=255)
                                                                                 aCGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCGTA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAAt                                                         <  2:27428/141‑1 (MQ=255)
                                                                                 aCGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAAt                                                         <  1:592469/141‑1 (MQ=255)
                                                                                 aCGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAAt                                                         >  1:1016897/1‑141 (MQ=255)
                                                                                     cATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAATTAAc                                                     <  1:656979/141‑1 (MQ=255)
                                                                                     cATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAATTAAc                                                     >  2:1078541/1‑141 (MQ=255)
                                                                                     cATACTTGCAAAGTAGATAATACAGAAAATCCAAAGTTGCGCTATCATA‑CGCGGCTTGGGATTTTTTCGTTTTAATGCTAGGTATAAGTCCATTGTTTGTCCTGTAAAACATGTTGATAAAAAGGTTCTTGGCCAGTTatc                                                     >  2:1078543/1‑139 (MQ=255)
                                                                                       tACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAATTAACtt                                                   >  1:534828/1‑141 (MQ=255)
                                                                                       tACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAATTAACtt                                                   >  1:458201/1‑141 (MQ=255)
                                                                                        aCTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAATTAACTTg                                                  >  2:516688/1‑141 (MQ=255)
                                                                                               aaGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAATTAACTTGATatcat                                           >  2:413178/1‑141 (MQ=255)
                                                                                               aaGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAATTAACTTGATatcat                                           >  1:480387/1‑141 (MQ=255)
                                                                                                 gTAGATAATACAGAAAATCTCAAGTTGCGATATCATA‑CGCAGCTAGGGAGTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAGAACATGTTGATAAAAAGGATCTTGGCCAATTAACTTGATATCATCt                                         >  2:1006681/1‑141 (MQ=255)
                                                                                                   aGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAATTAACTTGATATCATCTGc                                       <  1:751837/141‑1 (MQ=255)
                                                                                                   aGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAATTAACTTGATATCATCTGc                                       <  2:473545/141‑1 (MQ=255)
                                                                                                   aGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAATTAACTTGATATCATCTGc                                       <  1:823395/141‑1 (MQ=255)
                                                                                                   aGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAAAAAATTGATAAAAAGGATCTTGGCCAATTAACTTGATATCATCTGc                                       <  2:364846/141‑1 (MQ=255)
                                                                                                     ataataCAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAATTAACTTGATATCATCTGCaa                                     <  1:953404/141‑1 (MQ=255)
                                                                                                      taataCAGAATATCCCATGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAATTAACTTGATATCATCTGCAAg                                    <  2:1053610/141‑1 (MQ=255)
                                                                                                         taCAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAATTAACTTGATATCATCTGCAAGTGc                                 <  2:797988/141‑1 (MQ=255)
                                                                                                          aCAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAATTAACTTGATATCATCTGCAAGTGCt                                <  2:361335/141‑1 (MQ=255)
                                                                                                          aCAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAATTAACTGGATATCATCTGCAAGTGCt                                <  1:89695/141‑1 (MQ=255)
                                                                                                             gAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAATTAACTTGATATCATCTGCAAGTGCTTCa                             >  2:272958/1‑141 (MQ=255)
                                                                                                                   cccAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAATTAACTTGATATCATCTGCAAGTGCTTCAACTTCa                       >  1:84542/1‑141 (MQ=255)
                                                                                                                    ccAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAATTAACTTGATATCATCTGCAAGTGCTTCAACTTCAt                      <  1:681378/141‑1 (MQ=255)
                                                                                                                    ccAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGGTGATAAAAAGGATCTTGGCCAATTAACTTGATATCATCTGCAAGTGCTTCAACTTCAt                      >  1:514660/1‑141 (MQ=255)
                                                                                                                    ccAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACAAGTTAAAAAAAAGGATCTTGGCCAATTAACTTGATATCATCTGCAAGTGCTTCAACTTCAt                      <  2:140075/141‑1 (MQ=255)
                                                                                                                     cAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAATTAACTTGATATCATCTGCAAGTGCTTCAACTTCATc                     >  1:711894/1‑141 (MQ=255)
                                                                                                                      aaGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAATTAACTTGATATCACCTGCAAGTGCTTCAACTTCATCt                    >  2:137294/1‑141 (MQ=255)
                                                                                                                          ctgcgtatgatatcGCAACTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAATTAACTTGATATCATCTGCAAGTGCTTCAACTTCATCTaaa                 >  1:216786/14‑141 (MQ=255)
                                                                                                                           gCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAATTAACTTGATATCATCTGCAAGTGCTTCAACTTCATCTAAatg               <  1:242038/141‑1 (MQ=255)
                                                                                                                            gcgtatgatatcGCAACTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAATTAACTTGATATCATCTGCAAGTGCTTCAACTTCATCTAAatg               >  1:231306/12‑141 (MQ=255)
                                                                                                                              atatCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAATTAACTTGATATCATCTGCAAGTGCTTCAACTTCATCTAAatgatg            <  1:43873/141‑1 (MQ=255)
                                                                                                                               tatgatatcGCAACTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAATTAACTTGATATCATCTGCAAGTGCTTCAACTTCATCTAAatgatg            >  2:1074431/9‑141 (MQ=255)
                                                                                                                                  cATA‑CGCAACTTGGGATTTTTTCGTTTTAATAATAGGTATAAGCCCATTGCTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAATTAACTTGATATCATCTGCAAGTGCTTCCACTTCATCTAAATGATGGGTa        <  1:269349/141‑1 (MQ=255)
                                                                                                                                       cGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAATTAACTTGATATCATCTGCAAGTGCTTCAACTTCATCTAAATGACGGGAAGTTa    >  1:154119/1‑141 (MQ=255)
                                                                                                                                         cAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAATTAACTTGATATCATCTGCAAGTGCTTCAACTTCATCTAAATGATGGGTAGTTAat  >  1:867723/1‑141 (MQ=255)
                                                                                                                                         cAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGAACTTGGCCAATTAACTTGATATCATCTGCAAGTGCTTCAACTTCATCTAAATGATGGGTAGTTAat  >  2:448261/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                          
AGTACGGGTAATTCAACAGCGCCTCACGTACACTTCCAACGTATGTCTGGTGGCATCGGTAATCAATATGCAGTAGACCCAACGTCATACTTGCAAAGTAGATAATACAGAAAATCCCAAGTTGCGATATCATA‑CGCAGCTTGGGATTTTTTCGTTTTAATAATAGGTATAAGTCCATTGTTTGTCCTCTAAAACATGTTGATAAAAAGGATCTTGGCCAATTAACTTGATATCATCTGCAAGTGCTTCAACTTCATCTAAATGATGGGTAGTTAAT  >  CP000730/323084‑323360

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