Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 35,850 C→T 31.0% intergenic (+81/+27) rlmH → / ← tpn hypothetical protein/IS431mec transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073035,8500CT31.0% 144.5 / 52.6 71intergenic (+81/+27)rlmH/tpnhypothetical protein/IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base C (24/24);  new base T (10/12);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 8.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGC  >  CP000730/35713‑35960
                                                                                                                                         |                                                                                                              
aTTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCAc                                                                                                             <  2:255112/141‑1 (MQ=255)
aTTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCAc                                                                                                             >  1:97949/1‑141 (MQ=255)
 ttGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACt                                                                                                            <  1:81718/141‑1 (MQ=255)
         gTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATc                                                                                                    >  1:30971/1‑141 (MQ=255)
         gTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATc                                                                                                    >  1:335287/1‑141 (MQ=255)
           aTTTAATGTAGAAGCTATTCAGTATAAAGTTTTTTTTTTTTTTTGATTACTAGAACAAGAGATAAAATAAAATTGTGTTCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCAt                                                                                                  >  2:183847/1‑141 (MQ=25)
           aTTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCAt                                                                                                  >  1:418786/1‑141 (MQ=255)
           aTTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCAt                                                                                                  <  1:381431/141‑1 (MQ=255)
            tttAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATg                                                                                                 >  2:160546/1‑141 (MQ=255)
              tAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTc                                                                                               <  1:410352/141‑1 (MQ=255)
                           aTTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTaa                                                                                  >  2:271947/1‑141 (MQ=255)
                           aTTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTaa                                                                                  >  2:62521/1‑141 (MQ=255)
                                 tataAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCt                                                                            <  1:214523/141‑1 (MQ=255)
                                       gTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATCGTGGGACACTCTGTTTGCATGTCAGTGTTCGCGTAACTTGCTAGCatg                                                                      >  1:205252/1‑141 (MQ=37)
                                       gTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCatg                                                                      >  2:348203/1‑141 (MQ=255)
                                       gTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCatg                                                                      >  1:14571/1‑141 (MQ=255)
                                       gTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCatg                                                                      >  1:15322/1‑141 (MQ=255)
                                       gTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCatg                                                                      >  2:407848/1‑141 (MQ=255)
                                        gttttGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCatga                                                                     >  1:410010/2‑141 (MQ=255)
                                           ttGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGc                                                                  <  2:443251/141‑1 (MQ=255)
                                             ttttaattaTGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTa                                                                >  2:119648/9‑141 (MQ=255)
                                             ttctatttttGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTa                                                                >  2:115869/10‑141 (MQ=255)
                                             ggCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTa                                                                >  2:293968/1‑141 (MQ=255)
                                             ggCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTa                                                                <  2:264902/141‑1 (MQ=255)
                                             ggCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTa                                                                >  2:378104/1‑141 (MQ=255)
                                              gCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTaa                                                               >  2:388464/1‑141 (MQ=255)
                                              gCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTaa                                                               >  2:354678/1‑141 (MQ=255)
                                               ttttttttGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAAt                                                              >  2:396431/8‑141 (MQ=255)
                                                        aTTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCa                                                     >  1:332871/1‑141 (MQ=255)
                                                         ttACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCAt                                                    >  2:396499/1‑141 (MQ=255)
                                                          tACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATg                                                   >  2:447869/1‑141 (MQ=255)
                                                                     gagCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCg                                        <  2:66418/141‑1 (MQ=255)
                                                                     gagCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGAAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCg                                        <  2:400051/141‑1 (MQ=255)
                                                                               aaaatTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGaa                              >  1:27870/1‑141 (MQ=255)
                                                                                aaatTGTGATCTAATAACAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTGATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAag                             <  1:57848/141‑1 (MQ=255)
                                                                                  atTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAagag                           <  2:365240/141‑1 (MQ=255)
                                                                                  agagtattttaataaAAATGATGTTTCTTTTTTAAAGTAAATAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTGACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                         >  2:459699/10‑141 (MQ=11)
                                                                                    tgtgATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                         <  1:268911/141‑1 (MQ=255)
                                                                                      tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                     >  2:378866/6‑141 (MQ=11)
                                                                                      tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                     <  2:237572/136‑1 (MQ=11)
                                                                                      tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                     <  2:48009/136‑1 (MQ=11)
                                                                                      tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                     <  2:22792/136‑1 (MQ=11)
                                                                                      tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                     <  2:381418/136‑1 (MQ=11)
                                                                                      tagtttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTGATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCCTGGCTAAAAGCCGTAGAGCTGAAGAGACCTGc                     <  2:319034/136‑1 (MQ=11)
                                                                                       attttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTCGCAGGTTGCTAGTTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg                    <  2:323534/137‑1 (MQ=11)
                                                                                       attttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg                    <  2:228964/137‑1 (MQ=11)
                                                                                       attttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg                    <  2:6572/137‑1 (MQ=11)
                                                                                       attttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg                    <  2:407987/137‑1 (MQ=11)
                                                                                        ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                   >  2:257345/4‑141 (MQ=11)
                                                                                         tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAGCTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                  <  1:455586/139‑1 (MQ=11)
                                                                                         tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                  >  1:269923/3‑141 (MQ=11)
                                                                                         tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                  >  2:382213/3‑141 (MQ=11)
                                                                                         tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                  >  2:456938/3‑141 (MQ=11)
                                                                                         tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                  >  2:395541/3‑141 (MQ=11)
                                                                                         tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                  >  1:116691/3‑141 (MQ=11)
                                                                                         tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAAATTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTCGATCTGAAGAGACCTGCGGt                  >  1:92603/3‑141 (MQ=11)
                                                                                          cTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                   <  2:79064/141‑1 (MQ=40)
                                                                                          cTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                   <  1:396431/141‑1 (MQ=40)
                                                                                          cTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                   <  2:168743/141‑1 (MQ=40)
                                                                                           taataaAAATAGAAGGTTCTGTTGTAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                  <  2:138603/141‑1 (MQ=255)
                                                                                           taataaAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                  <  2:332871/141‑1 (MQ=40)
                                                                                           taataaAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                  <  2:215480/141‑1 (MQ=40)
                                                                                           taataaAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                  <  1:429799/141‑1 (MQ=40)
                                                                                           taataaAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                  <  2:101320/141‑1 (MQ=40)
                                                                                           taataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTc                <  2:457404/141‑1 (MQ=11)
                                                                                           taataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTc                <  1:224968/141‑1 (MQ=11)
                                                                                             ataaAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTc                <  2:320463/141‑1 (MQ=38)
                                                                                               aaaaaTAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCtt              <  2:433373/141‑1 (MQ=37)
                                                                                                    tAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCTCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTtat         <  2:244208/141‑1 (MQ=17)
                                                                                                    tAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTtat         <  1:300467/141‑1 (MQ=32)
                                                                                                          gtacTGTTCCAAAGTAAAAAATCTAGCGAATCACTAATATATCACGTCAGTGTTCGGTCAACTTGCTAGCATGATTCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGATCTTTTTATACAGAGc  >  2:391649/4‑141 (MQ=1)
                                                                                                                                         |                                                                                                              
ATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGC  >  CP000730/35713‑35960

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