Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 110,651 C→T 18.3% D170D (GAC→GAT USA300HOU_0109 → staphylococcal tandem lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730110,6510CT18.3% 139.4 / 26.0 71D170D (GAC→GATUSA300HOU_0109staphylococcal tandem lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (25/32);  new base T (8/5);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 3.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-02

AAGATATCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAA  >  CP000730/110513‑110785
                                                                                                                                          |                                                                                                                                      
aaGATATCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACg                                                                                                                                       <  1:370257/140‑1 (MQ=255)
  gATATCCGGTTAAAATGGTCTATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGga                                                                                                                                     <  1:101542/140‑1 (MQ=255)
    tatCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTAGAAAAATTATAAAAACGgaga                                                                                                                                   >  1:11713/1‑140 (MQ=255)
          gTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAAAATAGCGACTTTAAAAATTTGACAAAATATAAAGACGGAGATATTTc                                                                                                                             >  2:131413/1‑140 (MQ=255)
            tAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCat                                                                                                                           >  1:28343/1‑140 (MQ=255)
             aaaaTGGTCGATAATAAAATCATTCCAACAAAATAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCata                                                                                                                          <  1:598781/140‑1 (MQ=255)
              aactGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATCAAAAAAGAAATCGAAAACTCTAAGTTCTTCGTTCAATATCGCGACTTTAACAATTTGAAAACTTATAAAACCGGAGATATTTCctat                                                                                                                         >  2:141774/4‑140 (MQ=17)
              aaaTGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCatat                                                                                                                         <  1:645112/140‑1 (MQ=255)
                aTGGTCGATAATAAAATCATTCCAACAAAAGAAATCAGATATGAATAAATAAAAAAAGAAATCGAAAACTTTACGTGCGTTGTTCAATATGCAGACTTTGAAAATGGCAAAAATGATAAAAAAGGGGACAtttaaaaaaa                                                                                                                       >  2:546563/1‑133 (MQ=9)
                 tGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGATATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTCTAAAGACGGAGATATTTCATATAAt                                                                                                                      >  1:439229/1‑140 (MQ=255)
                 tGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAAt                                                                                                                      >  2:610214/1‑140 (MQ=255)
                  ggTCGATAATAAAATCATTCCAACAAAAGNAATTAAAGATGNAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGCCGACTTTTAAAATTTGAAAAAGTACAAAAACGGAGATATGTCANATaagc                                                                                                                     >  2:17522/1‑138 (MQ=17)
                      aataataAAATCATTCCAACACGAGAAAAAAAAGATGAAAAAATAAAAAAAGAAATCTAAAACTTGAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCaga                                                                                                                 <  1:194774/139‑1 (MQ=17)
                       ataataAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATCTGGCGCCTATAAACGTTTGCCCAATGATAAAGACGGAGATATTTCATATAATCCagag                                                                                                                <  1:350002/140‑1 (MQ=14)
                       ataataAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCagag                                                                                                                <  1:469277/140‑1 (MQ=255)
                        taataaAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAgg                                                                                                               <  1:611363/140‑1 (MQ=255)
                            aaaaTCATTCCAACAAAAGAAATTAAAGATGAAAAAATACAAAAAGAAATCGAACACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGcc                                                                                                           >  1:141113/1‑140 (MQ=255)
                               aTCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAAATATAAAGACGGAGATATTTCATATAATcacgagtagaccaa                                                                                                        >  2:610984/1‑127 (MQ=17)
                                     ccAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAAAGGTGCACAGGTATTc                                                                                                  >  2:432240/1‑140 (MQ=255)
                                       aacaaAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGCAAACTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCgg                                                                                                <  2:476326/140‑1 (MQ=255)
                                             aGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAat                                                                                          <  1:436556/140‑1 (MQ=255)
                                               aaaTTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAatat                                                                                        <  2:353231/140‑1 (MQ=255)
                                                 aTTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAAACGGAGATATTTCATATAATCCAGAGATGCCGAGTTATTCGGCTAAATAACa                                                                                      >  2:260746/1‑140 (MQ=255)
                                                 aTTAAAGATAAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTATTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGATGCCGAGTTATTCGGCTAAATATCa                                                                                      >  2:338750/1‑140 (MQ=255)
                                                    aaaGATGAAAAAATAACAAGAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAtt                                                                                   >  1:405838/1‑140 (MQ=255)
                                                    aaaGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAtt                                                                                   >  1:552465/1‑140 (MQ=255)
                                                        aTGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACt                                                                               <  1:703166/140‑1 (MQ=255)
                                                           aaaaaaTAAAAAAAGAAATCGAAAAATTCAAGTTCTTCATTCAATAGTGAGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAat                                                                            <  1:561128/140‑1 (MQ=14)
                                                             aaaaaaaaaaATNAAAAGAAAAACTTTAANTAATTTGTTCAATATTGCNACATTAAAAATTTAAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCGGAGTTATTCGGCTAAATATCAATTAACTAatga                                                                          <  1:17522/140‑1 (MQ=14)
                                                                aTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAAtgatga                                                                       <  2:631985/140‑1 (MQ=255)
                                                                     aaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAATTTTTAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtata                                                                  <  1:258139/140‑1 (MQ=255)
                                                                           aaTCGAAAACTTTAAGTTCTTTGTTCAATATTGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaa                                                            <  2:433052/140‑1 (MQ=255)
                                                                                  aaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATAAATTAACTAATGATGGTTATAATGTAAAACAAtt                                                     >  2:116802/1‑140 (MQ=255)
                                                                                    cTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTAc                                                   <  2:28343/140‑1 (MQ=255)
                                                                                     ctttaGTTCTTTGTACAACATGGCGACTTTAAACATTTGAAAAATGATAAAGACGGAGATAGTTCATATAATCCAGAGGTGCCCACCTATTCGGCTAAATATCAACTACCTATTGATGATTATAATGTAAAACAATTACg                                                  <  2:442362/136‑1 (MQ=17)
                                                                                      ttAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGc                                                 >  1:661565/1‑140 (MQ=255)
                                                                                      ttAAGTGCTGTGTGGAAAGTGGCCAGTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGc                                                 <  2:604233/140‑1 (MQ=17)
                                                                                       tAAGTTCTTTGTTCAATATTGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCa                                                <  2:664099/140‑1 (MQ=255)
                                                                                       tAAGTTCTTCGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCa                                                <  1:10695/140‑1 (MQ=255)
                                                                                       tAAGTACTTTGTTCAATATAGCGACTTTAAAAATTTTAAAAATTATAAATATGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATGCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCa                                                <  1:10763/140‑1 (MQ=17)
                                                                                        aaGTTCTTTGTTCAATATGGCTACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaa                                               <  2:557728/140‑1 (MQ=255)
                                                                                         aGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTAAGCaaa                                              >  2:722130/1‑140 (MQ=255)
                                                                                         aGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATAACAATTAACTAATGATGATTATAATGTAAAACAATTACGCaaa                                              >  2:256915/1‑140 (MQ=255)
                                                                                             cTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgat                                          <  1:603487/140‑1 (MQ=255)
                                                                                             cTTTGTTCAATATGGCGACGTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgat                                          <  1:216365/140‑1 (MQ=255)
                                                                                                   tCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgatatgata                                    <  1:187970/140‑1 (MQ=255)
                                                                                                   tCAATATGGAGATTTCAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGAAGAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgatatgata                                    <  1:20988/140‑1 (MQ=17)
                                                                                                              aCTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATAAACCGACGAGt                         >  2:575040/1‑140 (MQ=255)
                                                                                                                tttAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTaa                       <  2:595658/140‑1 (MQ=11)
                                                                                                                tttAAAAATTTGAAAAATTATAAAGACGGATATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTaa                       <  1:95532/140‑1 (MQ=255)
                                                                                                                  tAAAAATTTGAAAAATTATAAAGACGGAGATAGTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAAGTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAg                     <  2:405838/140‑1 (MQ=255)
                                                                                                                   aaagatttaatgaaTTATAAAGATGGAGATATATCATATAATCCAGGGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTAAAATGAAAAACAATTACGCAAAAGATATGATATACCGACGAGAAAAGc                    >  1:256829/13‑140 (MQ=11)
                                                                                                                   aaagatttaatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGc                    <  1:397488/128‑1 (MQ=11)
                                                                                                                   aaaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGc                    >  1:99895/1‑140 (MQ=255)
                                                                                                                   aaaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTATGCTAAATATCAATTAACTAATAATGATTATAATGTACAACAGTTACGCAAAAGATATGATCTACCGACAACTAAAGc                    >  2:171744/1‑140 (MQ=17)
                                                                                                                   aaaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATACAGAGGTGCCGAGTTATTTGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAAAAATTAAGCAAGAGGTAAGTAATACCGACGAGTAAAgt                    >  1:14434/1‑139 (MQ=18)
                                                                                                                    aagatttaatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCt                   >  2:227635/12‑140 (MQ=11)
                                                                                                                     agatttaatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAAAGTAAAACAATTAAGCAAAAAGTATGATATACAGACGAGAAAAGCTc                  >  1:188842/11‑140 (MQ=11)
                                                                                                                      gatttaatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCGCTTGACTGATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTcc                 >  1:168002/10‑140 (MQ=11)
                                                                                                                      gatttaatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTAGGAAAAAGATATGATATAACGAAGAGTAAAGCTcc                 >  1:216453/10‑140 (MQ=11)
                                                                                                                      gatttaatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGccac                 >  2:446556/10‑137 (MQ=11)
                                                                                                                      gatttaatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTcc                 <  1:434495/131‑1 (MQ=14)
                                                                                                                      aaTTTTAAAAAGTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTcc                 <  2:389372/140‑1 (MQ=255)
                                                                                                                      aaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGGGGTGCCGAGTTATTCGGCTAAATACCAAGTAACTAATGATGATTATAATGTAAAACAATTACGCAACAGATATGATATACCGACGAGTAAAGCTcc                 <  2:269228/140‑1 (MQ=255)
                                                                                                                      aaTTTAAAAAATTAAAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTACTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTcc                 <  1:131225/140‑1 (MQ=255)
                                                                                                                         ttGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCaaa              >  1:10908/1‑140 (MQ=255)
                                                                                                                           aatgaaTTATAAAGATGGAGACATATCATAAAATCCAGAGGGGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGt            <  2:243777/136‑1 (MQ=17)
                                                                                                                            atgaaTTATAAAGATGAAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGAAAAGCTCCAAAGtt           >  1:432976/4‑140 (MQ=14)
                                                                                                                             aaaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCTAAGTTa          >  2:524688/1‑140 (MQ=255)
                                                                                                                                aTTATAAAGACGGAGATATTTCATATAATCAAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATTCAAACGAGTAAAGCTCCaacaatctaa       >  2:724624/1‑132 (MQ=17)
                                                                                                                                     aaaGATGGAGCTATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaaa  >  2:492527/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                                      
AAGATATCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAA  >  CP000730/110513‑110785

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: