Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 110,651 C→T 15.7% 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,6510CT15.7% 152.3 / 14.8 71D170D (GAC→GATUSA300HOU_0109staphylococcal tandem lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (24/35);  new base T (4/8);  total (28/43)
Fisher's exact test for biased strand distribution p-value = 7.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

TATCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTT  >  CP000730/110517‑110790
                                                                                                                                      |                                                                                                                                           
tatCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGat                                                                                                                                       >  2:444125/1‑141 (MQ=255)
tatCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGat                                                                                                                                       >  2:759245/1‑141 (MQ=255)
  tCCGGTTAAAATGGTCGATCCATTTTGCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGatat                                                                                                                                     >  1:769101/1‑141 (MQ=18)
  tCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGTCGgagctct                                                                                                                                     >  2:575317/1‑137 (MQ=255)
  tCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGatat                                                                                                                                     >  2:766676/1‑141 (MQ=255)
    cGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATAttt                                                                                                                                   <  1:318086/141‑1 (MQ=255)
       ttAAAATGGTCGAAAATAAAATCATACCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCat                                                                                                                                <  2:177719/141‑1 (MQ=255)
        tAAAACGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGAAGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCata                                                                                                                               <  2:602119/141‑1 (MQ=255)
        tAAAACAGTCAATAATAAAACTATTCCAAAGAATAAAATTAAAGATTAACAAATAATAAAAGAAATCTAAAACTTTAAGTTCGTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCata                                                                                                                               <  2:704983/141‑1 (MQ=17)
           aaTGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGACGCAAAAATAAACACAGACATTGAAAAATAGAACCTCTTTGTTCAAGATAGCGACTTGAAAGATTGGAAAAGTTATAAAGACGGAGATATTTCATATaa                                                                                                                            >  2:765191/1‑141 (MQ=14)
             tGGTCGATAATAAAATCATTCCAATAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATc                                                                                                                          >  1:494431/1‑141 (MQ=255)
                tcaataataAAAACATTCCAACAAAAGACATTAAAGATAAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAAGATTGCGACTTTAAAAATTTGAAAAATTATAAAGACGAAGATATATTATATAATCCag                                                                                                                       <  2:51342/138‑1 (MQ=17)
                    taataaAATCATTCCAACAAAAGAAGTAAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCATTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAAACAGAGATATTTAATATAATCCaaaaaa                                                                                                                   >  1:394735/1‑136 (MQ=17)
                            ccATTCCAACAAAAAAAATCAAAAAACAAAAAATAAAAAAAGAAATCGAAAACCTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGtt                                                                                                           <  2:599684/140‑1 (MQ=17)
                                   aacaaAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGc                                                                                                    <  2:147373/141‑1 (MQ=255)
                                    gcacaaGAAATTAAAGATGAAAAAATAAAAAAAGAAAACGAAAACTTTCAGTTCTTTATTAAATAAGGGTACTTGAAAAATTAGAAAAATTATAAAGACGGAGATATTTCATAAAATCAGGAGGTGACGAGGTAttatgtt                                                                                                   >  2:193058/5‑136 (MQ=17)
                                    acgaaaGAAATGAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCt                                                                                                   <  1:247110/138‑1 (MQ=255)
                                    acaaAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCt                                                                                                   <  2:232181/141‑1 (MQ=255)
                                    acaaAAGAAAACAAATATGTAAAAAGAAAAAAAGAGATCGAAAACTTTAAGGTCTGTGTTCAATATGGCGACTTTAAAAATATGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGACATATACGGCt                                                                                                   >  2:184260/1‑141 (MQ=17)
                                     caaAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTa                                                                                                  >  2:24268/1‑141 (MQ=255)
                                      aaaaGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTaa                                                                                                 <  1:465006/141‑1 (MQ=255)
                                        aaGAAATTAAAGAGGAAAAAATAAAAAAACAAATCCAAAACTTTAAATTCTTTGTTCAATATGGCGAATAAAAAAAAAAAAAAAATTATAAAGACGGAGATATTTCATATAATCCAGTGGTGCCGAGTTATTCGGCTAAat                                                                                               <  2:748688/141‑1 (MQ=17)
                                                  aGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAAc                                                                                     <  1:26001/141‑1 (MQ=255)
                                                  aGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAAc                                                                                     <  2:78005/141‑1 (MQ=255)
                                                  aGATAAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAAc                                                                                     <  2:220608/141‑1 (MQ=255)
                                                   gATGAAAAAATAAAAAAAGAAATCGAAAAATTAGAGTTCTTTGATAAATAGGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAAACAAAAGGTGCAGAGTTATTCGGCTAAATATAAATTAAca                                                                                    >  2:336073/1‑140 (MQ=17)
                                                    aagcaacaaTAAAAAATGCAACCGAAAATATGAAGTTCTTTGTGCAATATGGCGACTGTAAAAATTTGAAAAATTATAAATACAGAGATATTTCAGATAATCCATAGGGGCCGAGTTATTCGGCTAAATATCAATTAACTa                                                                                   <  2:640528/134‑1 (MQ=17)
                                                      gAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAat                                                                                 <  1:311789/141‑1 (MQ=255)
                                                       aaaaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGCGGTGCCGAGTTGTTCGGCTAAATATCAATTCACTAatg                                                                                >  2:650898/1‑141 (MQ=255)
                                                          aaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgatg                                                                             <  1:759245/141‑1 (MQ=255)
                                                          aaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgatg                                                                             <  1:444125/141‑1 (MQ=255)
                                                            aTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGACAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATTGCAATCAACTAATgaggat                                                                           >  2:172766/1‑141 (MQ=14)
                                                                 aaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTATCTAATGATGATTATaa                                                                      <  1:383059/141‑1 (MQ=255)
                                                                     gAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTa                                                                  >  2:758350/1‑141 (MQ=255)
                                                                     aaaaTCTAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTa                                                                  <  1:390983/140‑1 (MQ=255)
                                                                             aaaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAtt                                                          >  1:140765/1‑141 (MQ=255)
                                                                             aaaCTTTAAGTTCTTAGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAtt                                                          >  1:575525/1‑141 (MQ=255)
                                                                                cTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTAAAAAGACGAAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATGATAATGTAAAACAATTACg                                                       >  1:644359/1‑141 (MQ=255)
                                                                                 tttAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGc                                                      <  2:587219/141‑1 (MQ=255)
                                                                                 tttAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGc                                                      <  1:701788/141‑1 (MQ=255)
                                                                                        tCTGTGTTCAATATGGCGACTTTAAAAATTCACTAGATTTGAAAGACTTGTATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATCATTATAATCTAACCCCATTACGCACAAcat                                               <  2:542759/141‑1 (MQ=9)
                                                                                                 aaTATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATAc                                      >  1:314700/1‑141 (MQ=255)
                                                                                                       gCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACcgacga                                >  1:68929/1‑141 (MQ=255)
                                                                                                        cGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAg                               <  1:766676/141‑1 (MQ=255)
                                                                                                        cGACTTTAAAAATTTGAAAAAGTATAAAGACGTAGATATTTCATATAATCCAGAGGGGCCGAGTTATTCTGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAg                               <  1:738530/141‑1 (MQ=255)
                                                                                                           cTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTaa                            <  1:732286/141‑1 (MQ=255)
                                                                                                            tttAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTaaa                           <  2:633037/141‑1 (MQ=11)
                                                                                                              tAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGc                         <  2:394735/141‑1 (MQ=255)
                                                                                                                aagatttaatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTc                       >  1:289029/12‑141 (MQ=11)
                                                                                                                aaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTc                       >  2:361640/1‑141 (MQ=255)
                                                                                                                 agatttaatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGCATTACCCACGAGTAAAGCTcc                      >  2:310156/11‑141 (MQ=11)
                                                                                                                 aaaaTTGAAAAATAATAAAAACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTcc                      <  2:328481/141‑1 (MQ=255)
                                                                                                                 aaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTcc                      <  2:769101/141‑1 (MQ=255)
                                                                                                                   atttaattaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCaa                    <  2:104720/133‑1 (MQ=14)
                                                                                                                   atttaatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCaa                    <  2:2006/133‑1 (MQ=14)
                                                                                                                   atttaatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCaa                    <  2:56807/133‑1 (MQ=255)
                                                                                                                    attaatgaaTTAAAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAAAGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCaaa                   <  2:48219/134‑1 (MQ=14)
                                                                                                                     ttGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATCTACCGACGAGTAAAGCTCCAAAg                  >  2:191653/1‑141 (MQ=255)
                                                                                                                        aaaaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTa               <  2:597467/141‑1 (MQ=255)
                                                                                                                        aaaaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTa               <  2:212178/141‑1 (MQ=255)
                                                                                                                        aaaaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTa               <  2:42796/141‑1 (MQ=255)
                                                                                                                         tgaaTTATAAAGATGGAGCTCTATCATATAATCAAGAGGTGCCGAGGTATTCGGCGAAATTTCAATTAACTAATGAGGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTAt              >  2:534705/3‑141 (MQ=17)
                                                                                                                         tgaaTTATAAAGATGGAGATATATCATATAATCAAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGAAAAAGATATGATATACCGACGAGTAAAGCTACAAAGTTAt              >  2:708309/3‑141 (MQ=14)
                                                                                                                         aaaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTTTTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTAt              <  1:59804/141‑1 (MQ=255)
                                                                                                                            aTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCAACGAGTAAAGCTCCAAAGGTaaggt           >  2:53968/1‑137 (MQ=255)
                                                                                                                             ttATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGtt          <  2:289029/141‑1 (MQ=255)
                                                                                                                                   aGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAgg    <  2:428757/141‑1 (MQ=255)
                                                                                                                                   aGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAgg    <  1:260646/141‑1 (MQ=255)
                                                                                                                                   aGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAgg    <  2:33456/141‑1 (MQ=255)
                                                                                                                                    gACGGAGATATTTCATATAATCCAGAGGTGCCGAGTAATTCGGCTAAATATCAATTAATTAATGCTGATAAAAAGGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGATCCAAAGTTATTGTTAAAAGGt   >  1:534965/1‑141 (MQ=18)
                                                                                                                                     aCGGAGATATTTCATATAATCCAGCGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGtt  <  2:179290/141‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                           
TATCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTT  >  CP000730/110517‑110790

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: