Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 36,315 G→A 12.1% Y79Y (TAC→TAT)  USA300HOU_0026 ← IS431mec transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073036,3150GA12.1% 152.1 / 14.3 73Y79Y (TAC→TAT) USA300HOU_0026IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base G (37/27);  new base A (3/6);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 2.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.29e-01

AAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTC  >  CP000730/36180‑36455
                                                                                                                                       |                                                                                                                                            
aaaTTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGt                                                                                                                                         >  2:528858/1‑141 (MQ=32)
 aaTTGTTTAATGAGACGTTTGATAAACGCATATGCTGGATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCCCCATTTTCCTTTTATTTTGATGTACGtc                                                                                                                                        >  1:496123/1‑141 (MQ=17)
 aaTTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGtc                                                                                                                                        <  2:391854/141‑1 (MQ=32)
 aaTTGGTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCCATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGtc                                                                                                                                        <  2:60829/141‑1 (MQ=17)
  aTTGTTTAATGAGACGATTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGtcc                                                                                                                                       >  1:304844/1‑140 (MQ=17)
        ttaTGAGAAGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCaa                                                                                                                                 <  1:28464/139‑1 (MQ=255)
        tAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCaa                                                                                                                                 >  1:232590/1‑141 (MQ=32)
         aaTGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAAt                                                                                                                                >  1:379554/1‑141 (MQ=32)
          aTGAGACGTTTTATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATa                                                                                                                               <  2:24372/141‑1 (MQ=17)
          aTGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGGCTCATCACTa                                                                                                                               >  1:494296/1‑141 (MQ=17)
             agaCGTTTAATAAACGGATATGCTGAATGATTATCTCGTTGCCTACGCAACCAAAGATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGc                                                                                                                            <  1:232474/141‑1 (MQ=11)
               aCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTCCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGTAATCAATGAAAGAGCTACATTTTTCCTGAATTTTGATATACGTCTCATCAATACGCCa                                                                                                                          >  1:536405/1‑141 (MQ=11)
                  tttGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCTTCAATACGCCAttt                                                                                                                       >  1:70300/1‑141 (MQ=17)
                   ttGATAAACGCATAAGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCAtttt                                                                                                                      >  1:336644/1‑140 (MQ=17)
                     taTAAAGCCATATTCTTGATGTTTATGCCGTTCCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGta                                                                                                                    <  2:541692/140‑1 (MQ=11)
                     gATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGtg                                                                                                                    >  2:113425/1‑140 (MQ=17)
                                atGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTCCGTTTCAACAATACGCCATTTGTAATAAGCttttt                                                                                                         <  1:145630/141‑1 (MQ=255)
                                  gCTGAATGATTATCTCGTTGCTTACGCAATCCATTATCTTATGTATGACCCTCTGCATCTATGGCACCATATAAATACCTCCATTTTCCTTTTATTTTTATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTa                                                                                                       <  2:205022/141‑1 (MQ=11)
                                     gAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGc                                                                                                    <  1:421508/141‑1 (MQ=40)
                                     gAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGc                                                                                                    <  2:302756/141‑1 (MQ=38)
                                     gAATGATTATCGTGTTGCTTACGCAACCAAATATCTTATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGc                                                                                                    <  2:517391/141‑1 (MQ=255)
                                         gATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTGGTACGAAGCTTTTTTTTGCtttt                                                                                                >  1:464354/1‑141 (MQ=255)
                                                    tGCTTACGCAACCAAATATCTAATGTATGTCCCTCTACCTCAATGGAACGATATAAATAGCTCCATTTTCCTTTTCTTTTGATGAACGTTTAATCACTACGCCATTTGTAATAGAATTTTTTGTGTTTTTTGTTCTaacgt                                                                                     >  1:269360/1‑138 (MQ=11)
                                                       ttACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGa                                                                                  >  2:576262/1‑141 (MQ=40)
                                                       ttACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTAAGTCTTATGAAGACGCCTTTTGTAATCAGCTTTTTTATGGTTTTTCTTCCAAATTTGa                                                                                  >  1:1681/1‑141 (MQ=11)
                                                       gtACGCAACCAAAGATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGa                                                                                  <  1:171576/140‑1 (MQ=255)
                                                        tACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGATTTTTTCTGCGTTTTCTTCCCACTTTGAt                                                                                 >  2:426853/1‑141 (MQ=14)
                                                         aCGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTTATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATa                                                                                >  1:431712/1‑141 (MQ=255)
                                                         aCGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATa                                                                                >  1:527987/1‑141 (MQ=40)
                                                          cGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATAc                                                                               <  1:312544/141‑1 (MQ=40)
                                                          agCAACCAAATATCTAATGTATGTCTCTCTGCATCAATTGCACGATATAAATCGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATAc                                                                               <  1:355394/140‑1 (MQ=255)
                                                            cAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACaa                                                                             <  2:246764/141‑1 (MQ=255)
                                                            cAACCAAATATCTAATGTATGGCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACaa                                                                             <  2:435338/141‑1 (MQ=255)
                                                              aCCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAATCTTTTTTATGGTTTTTCTTCCAAATTTGATACaaaa                                                                           >  1:540058/1‑141 (MQ=255)
                                                              aCCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAACTTTTTTATGCTTTTTCCTCCAAATTTGATATaaaa                                                                           >  1:243191/1‑141 (MQ=255)
                                                                cAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAAtt                                                                         >  2:74738/1‑141 (MQ=255)
                                                                    tatCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTgggg                                                                     <  2:232590/141‑1 (MQ=255)
                                                                    tatCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGATCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCTTTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTTATACAAAATTgggg                                                                     >  2:336857/1‑141 (MQ=255)
                                                                     atCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGc                                                                    >  2:174817/1‑141 (MQ=255)
                                                                             tatgtCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCtt                                                            >  2:486745/1‑141 (MQ=255)
                                                                              ttgtCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTg                                                           <  1:61055/140‑1 (MQ=255)
                                                                                gtCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGac                                                         >  2:332576/1‑140 (MQ=255)
                                                                                 tCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAAc                                                        >  2:182778/1‑141 (MQ=21)
                                                                                 tCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAACTGGGGGCATCTTGtttcac                                                        >  1:380277/1‑137 (MQ=14)
                                                                                 tCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGTGGCATATTCTTGAAc                                                        >  1:142263/1‑141 (MQ=255)
                                                                                 tCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAAc                                                        <  1:123951/141‑1 (MQ=40)
                                                                                   ccTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAAccc                                                      <  1:377296/141‑1 (MQ=40)
                                                                                         cATCAATGGTACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATATGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCTTATCCTTGAACCCAACGGt                                                <  1:455784/141‑1 (MQ=255)
                                                                                             aaTGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCAATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACGAAATTGGGGCATATCCTTGAACCCaaagtgagcc                                            >  2:483970/1‑133 (MQ=14)
                                                                                                 gCACGTTATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTCTACCGtt                                        <  1:216200/141‑1 (MQ=255)
                                                                                                    cGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTTaa                                     >  1:68961/1‑141 (MQ=255)
                                                                                                    cGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGaa                                     >  1:96042/1‑141 (MQ=40)
                                                                                                      atatAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAatg                                   <  2:379554/141‑1 (MQ=40)
                                                                                                          aattAGCTCCATTTTCCTCTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAAtgatga                               <  1:483970/138‑1 (MQ=255)
                                                                                                           aaTAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGATATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGaa                              >  1:541872/1‑141 (MQ=11)
                                                                                                           aaTACCGCCATTTTCCTTTTATTTTGATGTACGTCGCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATCAGAAAAAACACTGGGGCATATTCTTGAACCCAACGGTAGACCGTTCAATGATGaa                              <  2:585779/141‑1 (MQ=11)
                                                                                                              aGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCACCGGTAGACCGTTGGATGATGAACGt                           >  1:544553/1‑141 (MQ=255)
                                                                                                              aGCTCCATTTTCCTGTGATTTTGATGTACGCCTCATCAATACGCCATTTGTAATAAGCGTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATGGAGGAACATTCTTCAACCCAACCGACAACCCTGGGATGATGAACGt                           >  2:559397/1‑141 (MQ=11)
                                                                                                                cTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGttt                         >  1:392754/1‑141 (MQ=255)
                                                                                                                   cATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTaca                      >  1:481481/1‑141 (MQ=40)
                                                                                                                    aTTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTTATACAACATTGGGGCATATTCTTGAACCCAACGGTAGACCGCTGCATGGTCACCGtttttac                     >  2:247622/1‑137 (MQ=14)
                                                                                                                    aTTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTacac                     >  1:525587/1‑141 (MQ=40)
                                                                                                                     ttttCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTAcacc                    <  1:38464/141‑1 (MQ=40)
                                                                                                                         ccTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGt                <  2:96042/141‑1 (MQ=40)
                                                                                                                         ccTTTTATTTTGACAAACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGTACGTTTACACCACGt                <  2:467641/141‑1 (MQ=11)
                                                                                                                           ttttATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTc              <  1:455602/141‑1 (MQ=255)
                                                                                                                            tttATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTcc             >  1:32683/1‑141 (MQ=40)
                                                                                                                             ttATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTccc            <  1:174817/141‑1 (MQ=40)
                                                                                                                             ttATTGTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTccc            <  1:300952/141‑1 (MQ=255)
                                                                                                                               tttttGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTTTGGTGTGTCTGCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCtt          <  1:17232/140‑1 (MQ=14)
                                                                                                                               aTTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCtt          <  2:461496/141‑1 (MQ=39)
                                                                                                                                    gATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGATTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGGAGACCGTTGAATGATGAACGTTTAAACAACGTGCCCTTGatat     >  1:507434/1‑141 (MQ=14)
                                                                                                                                    gATGTACGTCTCATCAATACGCAATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAAATGGGAAATATAATTACACCCAACGGTAACACGTTGAATGATGAACGTTTACACCACGTTCCCTTAatat     >  1:580722/1‑141 (MQ=14)
                                                                                                                                       gTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAGATTTGATACAAAATTGGGGAATATTCGTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTACCTTAATATTTc  >  2:336583/1‑141 (MQ=255)
                                                                                                                                       |                                                                                                                                            
AAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTC  >  CP000730/36180‑36455

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