Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 110,651 C→T 7.9% 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,6510CT 7.9% 186.5 / 5.3 68D170D (GAC→GATUSA300HOU_0109staphylococcal tandem lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (26/35);  new base T (2/3);  total (29/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.66e-01

AGATATCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTC  >  CP000730/110514‑110791
                                                                                                                                         |                                                                                                                                            
aGATATCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGATCAAAATGGCGACTTTAAAAATAAGAAAAATTATAAAGACGga                                                                                                                                           >  1:256354/1‑141 (MQ=255)
aGATATCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAAATTTAAGTCCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGAAGga                                                                                                                                           >  1:101213/1‑141 (MQ=255)
aGATATCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAAAAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTATAAAAATTTGAAAAATTATAAAGACGga                                                                                                                                           >  1:160162/1‑141 (MQ=255)
   tatCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGat                                                                                                                                        >  2:125397/1‑141 (MQ=255)
    atCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCGATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGata                                                                                                                                       <  1:522179/141‑1 (MQ=255)
    atCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGata                                                                                                                                       <  1:451438/141‑1 (MQ=255)
    atCCGATTAAAATGGTCGATAATAAAATAATTCAAACAAAAAAAAATAAAGATGAAAAAATAAAAAAGGAAATCGAAAATTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAACTTGAAAAATTATAAAGACGGAGata                                                                                                                                       <  1:606171/141‑1 (MQ=17)
            aaaaTGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAAATATACATACGGAGATATTTCataa                                                                                                                               >  1:680491/1‑140 (MQ=255)
                tGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATc                                                                                                                           >  2:373744/1‑141 (MQ=255)
                      ataataAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAgg                                                                                                                     >  1:761299/1‑141 (MQ=255)
                         ataaAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGc                                                                                                                  <  2:286260/141‑1 (MQ=255)
                          aaaaaTCATTCCAACAAAAGAAATGAAAGAAGAAAAAATAAAAAACGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGcc                                                                                                                 <  2:257394/140‑1 (MQ=255)
                           agaaTCATGCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAACCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTAGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCg                                                                                                                <  2:498549/139‑1 (MQ=255)
                            aaaTCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGa                                                                                                               <  1:616539/141‑1 (MQ=255)
                                cATTCCAACAAAAGAAATTAAAGAAGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTa                                                                                                           <  1:296303/141‑1 (MQ=255)
                                   tctaacaaAAGAAATTAAAGATGTACAGATAAAAAAGGAGATCGACCACTTTAATTTTTTTGTTCAATATGGCGACATTAAAACTTGTACAAATTAGACAGACGGACTAATTTCATAAAATCCGGGGGTGCCGAgatgttg                                                                                                        >  2:111221/4‑135 (MQ=9)
                                       aaaaaaGAAATTAAGGATGAAAAAATAAAAAAAGAAAATTATTACTTTAATTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCt                                                                                                    <  2:78692/139‑1 (MQ=17)
                                            aGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAata                                                                                               <  1:25302/141‑1 (MQ=255)
                                              aaaTTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTGATTCGGCTAAATATc                                                                                             <  1:565749/141‑1 (MQ=255)
                                                 ttAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAt                                                                                          >  2:511181/1‑141 (MQ=255)
                                                 atAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAt                                                                                          <  2:532423/140‑1 (MQ=255)
                                                  tAAAGATGAAAAAATAAAAAAAGAAGTCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAtt                                                                                         <  2:14672/141‑1 (MQ=255)
                                                  tAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAGAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAtt                                                                                         >  2:564932/1‑141 (MQ=255)
                                                   aaaGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTa                                                                                        >  2:478562/1‑141 (MQ=255)
                                                     aGATGAAAAAATAAAAAAAGAAATCGAAACATTTAAGTTATTTGTTCAATATGGAGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTAGGCTAAatctgacttcac                                                                                      >  2:434806/1‑131 (MQ=17)
                                                     aGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAAc                                                                                      <  1:8571/141‑1 (MQ=255)
                                                      gATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACt                                                                                     >  1:629170/1‑141 (MQ=255)
                                                          aaaaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatg                                                                                 <  1:412669/141‑1 (MQ=255)
                                                            aaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTCatgat                                                                               >  2:697782/1‑141 (MQ=255)
                                                             aaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgatg                                                                              <  1:187482/141‑1 (MQ=255)
                                                              aaTAAAAAAAGAAATCGAAAACTTTAAGTTCTGTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAAtgatga                                                                             <  1:256716/141‑1 (MQ=255)
                                                               aTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATgatgat                                                                            <  2:488889/141‑1 (MQ=255)
                                                                 aaaaaaaGAAATCGAAAACTCTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATta                                                                          <  1:9114/141‑1 (MQ=255)
                                                                  aaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                                                                         >  1:60834/1‑141 (MQ=255)
                                                                      aaGAAATCGAAATCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATg                                                                     >  1:536417/1‑141 (MQ=255)
                                                                      aaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATg                                                                     <  1:686057/141‑1 (MQ=255)
                                                                         aaaGCGAAAACTGTAAGTTCTTTGTTCAATATGGCGACTTTAAAAAGTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaa                                                                  <  1:316119/141‑1 (MQ=255)
                                                                             cggaaaCTTTGAGTTCTTTGTTCAATGTGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACa                                                              <  2:488030/138‑1 (MQ=255)
                                                                              gAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACaa                                                             >  1:696018/1‑141 (MQ=255)
                                                                               aaaaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                                                            <  1:235267/141‑1 (MQ=255)
                                                                                       aaGTTCTTTGTTCAATATGGCGACTTTCTAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaaa                                                    >  1:755421/1‑141 (MQ=255)
                                                                                         gTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaaaaa                                                  >  2:78275/1‑140 (MQ=255)
                                                                                         gTTCTTTGTTCAATAAGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAg                                                  >  1:564649/1‑141 (MQ=255)
                                                                                              ttGTTCCATATGCTGACGTTAAAAATGTGAAAAATTATAAAGACGGAAGTTTTTCATATAATCCAGAGGTGCCGAGTTATTCTGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTCCGCAAAAgatatg                                             <  2:266965/141‑1 (MQ=17)
                                                                                              ttGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACACAAAAgatatg                                             >  1:648340/1‑141 (MQ=255)
                                                                                                 ttCTATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgatatgata                                          <  2:340029/141‑1 (MQ=255)
                                                                                                  tCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgatatgatat                                         <  1:98928/141‑1 (MQ=255)
                                                                                                      tatGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACcg                                     >  1:63935/1‑141 (MQ=255)
                                                                                                      tatGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACcg                                     >  1:263298/1‑141 (MQ=255)
                                                                                                                    agatttaatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTcc                       <  1:410376/131‑1 (MQ=14)
                                                                                                                    aaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGATTAAAGCTcc                       <  2:464176/141‑1 (MQ=255)
                                                                                                                     aaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCa                      >  1:425181/1‑141 (MQ=255)
                                                                                                                      aTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATCTACCGACGAGTAAAGCTCCaa                     <  2:755421/141‑1 (MQ=255)
                                                                                                                      aTTTGAAAAATGATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCaa                     <  2:426068/141‑1 (MQ=255)
                                                                                                                        ttaatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAg                   <  1:643646/135‑1 (MQ=255)
                                                                                                                          taaaaaTTATAAAGACGGAAATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGtt                 <  2:497928/140‑1 (MQ=255)
                                                                                                                           atgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTa                >  2:482506/4‑141 (MQ=255)
                                                                                                                           aaaaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTGGGCTTAATAGCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTa                >  1:714796/1‑141 (MQ=255)
                                                                                                                           aaaaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATGCCGACGAGTAAAGCTCCAAAGTTa                <  1:485891/141‑1 (MQ=255)
                                                                                                                            cgcatcataAAGAGGGAGATATATCATATAATCCAGAGGTGCCGAGTTATCCCGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTAt               <  1:41116/135‑1 (MQ=17)
                                                                                                                            aaaaTTATAAAGCCGGAGATATTTCATATAATCCAGTGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTAt               <  2:619572/141‑1 (MQ=255)
                                                                                                                            aaaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTAt               <  1:373744/141‑1 (MQ=255)
                                                                                                                             gaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAGTTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAattatt              >  2:609137/2‑141 (MQ=14)
                                                                                                                                 tataAAGGCGGACATATTTCATATAATCCAGCGGCGCAGTGTTATCCGGCTGAACGTTTAATAACAAAATATGAGTATAAAATAAGACAATTACGCATAATATAGTATATACCGACGAGTAAAGCTCCAAAGTTATTGTTa          <  2:178479/141‑1 (MQ=255)
                                                                                                                                 tataAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTa          <  2:244336/141‑1 (MQ=255)
                                                                                                                                      aGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGCTATACCGACGAGTAAAGCTCCAAAGGTATTGTTAAAAgg     >  1:615666/1‑141 (MQ=255)
                                                                                                                                      aGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAGGTTATTGTTAAAAgg     >  1:226816/1‑141 (MQ=255)
                                                                                                                                      aGACGGAGACATTTCATAAAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGtaacaaag     <  2:315012/141‑8 (MQ=255)
                                                                                                                                         cGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGtta  >  1:331639/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                            
AGATATCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTC  >  CP000730/110514‑110791

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