Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,929,946 G→A 38.3% intergenic (‑1029/‑139) USA300HOU_1794 ← / → USA300HOU_1796 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,929,9460GA38.3% 78.9 / 88.5 100intergenic (‑1029/‑139)USA300HOU_1794/USA300HOU_1796hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (35/25);  new base A (20/18);  total (55/45)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

TAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAAATA  >  CP000730/1929807‑1930085
                                                                                                                                                |                                                                                                                                                  
taataaTAGATTCTCCTATCTCCCGGGCGTCAATT‑‑GACGAAGAGAGGACGTGCCTTACTTTGCTTCTCATTTTCGTTCACATCATAGCACCATTCATCAGTTGCTGTGCTGTTGCGTATTATTCTTATTGGCTTAGTAAAc                                                                                                                                                  <  2:115887/141‑1 (MQ=255)
 aagaatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACg                                                                                                                                                 <  2:227011/136‑1 (MQ=255)
 aaatatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACTAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACg                                                                                                                                                 <  2:816243/136‑1 (MQ=255)
  ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                                <  1:966583/141‑1 (MQ=255)
  ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                                <  2:339298/141‑1 (MQ=255)
  ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                                <  2:581846/141‑1 (MQ=255)
   tagtagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGc                                                                                                                                                <  1:1004251/141‑1 (MQ=255)
   tagtagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGc                                                                                                                                                <  1:1071881/141‑1 (MQ=255)
   tagtagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGc                                                                                                                                                <  1:216377/141‑1 (MQ=255)
    aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                                >  1:646883/3‑139 (MQ=255)
    aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                                <  2:646883/137‑1 (MQ=255)
     atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaac                                                                                                                                             >  2:1066255/2‑140 (MQ=255)
     atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGCAAACGCaac                                                                                                                                             >  2:851496/2‑140 (MQ=255)
      tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                                                                                                            <  1:280341/141‑3 (MQ=255)
      tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                                                                                                            <  1:571420/141‑3 (MQ=255)
      aagagTCGCCTATCTCTCATGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCacca                                                                                                                                            <  2:927093/140‑4 (MQ=255)
       agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                                                           >  2:595827/1‑138 (MQ=255)
       agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                                                           <  2:298789/141‑4 (MQ=255)
         agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTGTACTTCTTGGCTTAGTAAACGccacaaat                                                                                                                                         >  2:359396/1‑134 (MQ=255)
          gTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACaa                                                                                                                                       <  2:997173/141‑1 (MQ=255)
          gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                                                                                                        >  1:855455/1‑135 (MQ=255)
           tCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAGGTCGCCGCAAT‑‑ACaa‑‑a                                                                                                                                    >  1:315166/1‑141 (MQ=255)
           cagCCTATCTCTCAATCCTCAATTTATACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTCCTTATTGGCTAAGTAGACGCAAT‑‑ACaa‑‑a                                                                                                                                    <  1:467949/139‑1 (MQ=255)
            cGCCTATCTCTCAGTCGTCAACT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCAGCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTGGTGAACGCaacaaataga                                                                                                                                      >  2:758655/1‑133 (MQ=255)
            cGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACaa‑‑aa                                                                                                                                   <  2:201777/141‑1 (MQ=255)
             gCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGAATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACaa‑‑aaa                                                                                                                                  >  2:551636/1‑140 (MQ=255)
             gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                                                     <  1:741885/141‑10 (MQ=255)
             gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                                                     <  2:526243/141‑10 (MQ=255)
              ccTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGTCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCATT‑‑ACAA‑‑AATa                                                                                                                                 >  1:610897/1‑141 (MQ=255)
              ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacagctaatta                                                                                                                                    >  1:953996/1‑131 (MQ=255)
              ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                                                                                                                    >  2:748463/1‑131 (MQ=255)
               cTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAg                                                                                                                                <  1:285229/141‑1 (MQ=255)
               cTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAg                                                                                                                                <  1:603150/141‑1 (MQ=255)
               cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                                                                   >  2:870350/1‑130 (MQ=255)
               cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                                                                   >  2:318013/1‑130 (MQ=255)
               cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataagtaa                                                                                                                                   >  2:742876/1‑130 (MQ=255)
               cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataagtaa                                                                                                                                   >  2:742861/1‑130 (MQ=255)
                tATCTCTTAGCCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                                                  <  2:178232/141‑13 (MQ=255)
                tATCTCTCAGGCTTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCGCCAGTCATTATTGGCTGTGCTGTGGCTTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                                                  <  2:1064077/141‑13 (MQ=255)
                tATCTCTCAGGCGTCAATT‑‑GACGAAGGGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                                                  <  2:510179/141‑13 (MQ=255)
                tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGCTGCGTATTATACTGATGGTCTTAGTAAACGCaacaaataatcaat                                                                                                                                  >  2:882173/1‑129 (MQ=255)
                 aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTATACGCaacaaataattaata                                                                                                                                 <  1:607034/141‑14 (MQ=255)
                 aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                                                 >  1:279443/1‑128 (MQ=255)
                 aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                                                 >  1:361686/1‑128 (MQ=255)
                 aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                                                 >  2:474869/1‑128 (MQ=255)
                 aTCTCTCACGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATGATTTTCGTTCACATCATAGCACCAGTCATCAGTGGGTGTGCTGTTGCGTATTATACTTAGTGGCTTAGTAAACGCaacaaataattaata                                                                                                                                 >  2:815114/1‑128 (MQ=255)
                  tctctcAGGCTTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaatag                                                                                                                                <  2:681465/141‑15 (MQ=255)
                  tctctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTCTTGGCTAAGTAGACGCAAT‑‑ACAC‑‑AATAGGta                                                                                                                             >  2:897083/1‑140 (MQ=255)
                   ctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCAAC‑‑ACAT‑‑AAttaatagg                                                                                                                           >  2:340398/1‑134 (MQ=255)
                   ctctcAGGCGTCAATATAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCCAT‑‑ACAA‑‑AATAGGTGa                                                                                                                            >  1:820861/1‑141 (MQ=255)
                            gTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGTTTATCTTCGTTCAAATCCTATCTTCCGTGATCAGTGTCTGTGCTATCGCGTTTTTTACTTATTGTCTTACTAACCGCAAT‑‑ACAA‑‑AATAGGTGACATATGGcc                                                                                                                 <  2:458261/141‑1 (MQ=255)
                             tCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCAATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATCGGTGACATATAGCCg                                                                                                               >  1:1002530/1‑141 (MQ=255)
                                  ttAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAC‑‑ACTAGGTGACATATAGCCGCACCa                                                                                                          >  2:72035/1‑141 (MQ=255)
                                             gaggagGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAACAAAccc                                                                                               >  2:1071356/1‑141 (MQ=255)
                                              aggaggTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTATTATTGGTTTTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCCATAAAAATcccc                                                                                              >  1:526248/1‑141 (MQ=255)
                                               ggaggTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCt                                                                                             <  2:857951/141‑1 (MQ=255)
                                               ggaggTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCCGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATCGGTGACATATAGCCGCACCAATAAAAATCCCCt                                                                                             <  2:729976/141‑1 (MQ=255)
                                                gaggTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTc                                                                                            <  1:1101970/141‑1 (MQ=255)
                                                 aggTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCa                                                                                           >  2:903794/1‑141 (MQ=255)
                                                  ggTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCAc                                                                                          <  1:1035769/141‑1 (MQ=255)
                                                   cgataataAGGAAATGCTGGTCTGCGTTCACATCAATGCGCCCCTCATCAGCGGCTGCGCCTTTTCGTTTTCTTCTTACCGGCAAACTATACGCAAT‑‑AGAA‑‑AATATGTGACATATAGCTTCACCAATAAAAATCCCCTCAc                                                                                          <  2:857658/136‑1 (MQ=255)
                                                     gtATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTTGCTAAGTAGACTCAAT‑‑ACAG‑‑AATAGGTGACATATAGACGAACCAATAAAAATCCCCTGAcgaa                                                                                       >  1:194561/1‑138 (MQ=255)
                                                          taAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGc                                                                                    >  2:632019/1‑141 (MQ=255)
                                                          taAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTAACTACCGc                                                                                    >  2:621381/1‑141 (MQ=255)
                                                          taAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGAAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGc                                                                                    >  1:636235/1‑141 (MQ=255)
                                                           aaGGTAATGCTTGTTTGCGTCCACATCCTAGCACCAGTCATCAGTGGCTTTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACACATAGCCGCACCAATAACAATCCCCTCACTACCGCa                                                                                   <  2:213763/141‑1 (MQ=255)
                                                            aGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTTTCTCAGTAGACTCTAT‑‑ACGA‑‑ATTAGTGGCCATATATCCGCCCCAATAAAAATCCCCTCACTACCGCaa                                                                                  >  2:231315/1‑141 (MQ=255)
                                                            aGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCaa                                                                                  >  2:331773/1‑141 (MQ=255)
                                                            aGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCaa                                                                                  >  1:861554/1‑141 (MQ=255)
                                                            aGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCaa                                                                                  <  1:951174/141‑1 (MQ=255)
                                                                  tGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGGTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGAAACAATAAAAATCCCCTCACTACCGCAAATa                                                                               <  2:118213/141‑1 (MQ=255)
                                                                  tGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATa                                                                               <  1:719582/141‑1 (MQ=255)
                                                                   gATGCTGATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAg                                                                              <  1:390148/141‑1 (MQ=255)
                                                                        ttATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgg                                                                         >  2:988190/1‑141 (MQ=255)
                                                                        ttATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAAAAGTGAgg                                                                         >  1:1028616/1‑141 (MQ=255)
                                                                         tactttCGTTCACATCATCGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTGTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGCGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAggg                                                                        <  2:279408/138‑1 (MQ=255)
                                                                          aTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgggg                                                                       >  1:840814/1‑141 (MQ=255)
                                                                           tcttCGGTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                                                      <  1:897083/139‑1 (MQ=255)
                                                                               cGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTgg                                                                  >  2:778228/1‑141 (MQ=255)
                                                                               cGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTgg                                                                  >  1:444873/1‑141 (MQ=255)
                                                                               cGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACAGATAGCCGCACCAATAAAAATCCCCTCCCTACCGCAAAAAGTGAGGGGATTgg                                                                  >  1:21586/1‑141 (MQ=255)
                                                                                 ttCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATGGgtt                                                                >  1:597107/1‑140 (MQ=255)
                                                                                  tCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGgtgt                                                               >  1:622748/1‑141 (MQ=255)
                                                                                  tCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGgtgt                                                               >  2:183488/1‑141 (MQ=255)
                                                                                          tAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCCCACTACTGCAATAGTGTGGGGATTGGTGcatt                                                            >  1:216538/1‑140 (MQ=255)
                                                                                            gCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAAt                                                     <  1:342894/141‑1 (MQ=255)
                                                                                               ccAGTCATACATGGCTGTTCCATTTCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTGAAGACt                                                  >  2:1072782/1‑141 (MQ=255)
                                                                                                  gTCATCAGTGGCTGTGCCATTGCGTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAtt                                              >  1:389431/1‑141 (MQ=255)
                                                                                                     atcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAtttt                                            >  1:651664/1‑141 (MQ=255)
                                                                                                     atcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACACAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGCTTGGTGTATAAGTAAATACTTAtttt                                            >  1:1042337/1‑141 (MQ=255)
                                                                                                            gCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTc                                     <  2:840814/141‑1 (MQ=255)
                                                                                                               gtgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACATATTTTCGTTGTCTTa                                  >  1:185970/1‑141 (MQ=255)
                                                                                                                    aTTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATtat                             <  1:874452/141‑1 (MQ=255)
                                                                                                                      tGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATAc                           <  1:778228/141‑1 (MQ=255)
                                                                                                                        cGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTg                         <  1:988190/141‑1 (MQ=255)
                                                                                                                           ttttttCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTa                      <  2:861554/141‑1 (MQ=255)
                                                                                                                                   tATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCt              >  2:334082/1‑141 (MQ=255)
                                                                                                                                        gCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCtttttt         >  1:345864/1‑140 (MQ=255)
                                                                                                                                             gTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTaaaa    >  2:147867/1‑141 (MQ=255)
                                                                                                                                               aGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTAACGCAAATAGTGAGGGGATTGGTGGATAAGTAAATACTTATTTTCGTTGTCTTAATGATACTGCTAATTTTTCTTTTTGTAGAatg  >  1:352263/1‑140 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
TAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAAATA  >  CP000730/1929807‑1930085

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