Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,929,946 G→A 43.4% 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,9460GA43.4% 36.9 / 100.4 92intergenic (‑1029/‑139)USA300HOU_1794/USA300HOU_1796hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (22/30);  new base A (22/17);  total (44/48)
Fisher's exact test for biased strand distribution p-value = 2.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.89e-01

ATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAAATAT  >  CP000730/1929806‑1930086
                                                                                                                                                |                                                                                                                                                   
ataataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATGTTCGTTCGCAGCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTCAGTAg                                                                                                                                                     >  2:505018/1‑141 (MQ=255)
 taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCAGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAc                                                                                                                                                   >  2:191068/1‑141 (MQ=255)
  aataatAGAGTTGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACg                                                                                                                                                  >  1:307755/1‑141 (MQ=255)
  aataatAGAGTGGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACg                                                                                                                                                  <  2:239390/141‑1 (MQ=255)
   ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                                 <  1:810716/141‑1 (MQ=255)
   ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCGAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATTCTTATTGGCTTAGTAAACGc                                                                                                                                                 <  1:542625/141‑1 (MQ=255)
    taatagAGTCGCCTATCTCTCAGGCCTCAATT‑‑GAAGAAGAGAGGTGGTGCATTACTTTGTTAATCATTTTCCTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                                                                                                                <  2:379340/138‑1 (MQ=255)
     aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                                                                                                               >  2:650488/3‑141 (MQ=255)
     aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                                                                                                               <  1:111478/139‑1 (MQ=255)
     aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCCTAGTAAACGCaa                                                                                                                                               >  2:39152/3‑141 (MQ=255)
      atatagTCCCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaac                                                                                                                                              <  1:107199/137‑2 (MQ=255)
      atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaac                                                                                                                                              >  2:396951/2‑140 (MQ=255)
        agagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATa                                                                                                                                             >  2:836573/1‑141 (MQ=255)
         gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                                                                                                                           >  2:642553/1‑137 (MQ=255)
         gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                                                                                                                           <  1:376056/141‑5 (MQ=255)
         gagTCGCCGATGTGTCAGGCTGCAATCGAGGCGCAGAGGGGAGGTGTATAAGGGGATGCTTATTTTCGTTCACATCATAGCTCCAGTCATTAGTGGCTGTGCCAGTGCGTTTTTGTCTGATGGGCTAAGTAGACGCAATAc                                                                                                                                            <  2:171833/141‑1 (MQ=255)
          agTCGCCTATCTGTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCAGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                                                                                                          >  1:592453/1‑136 (MQ=255)
          agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCAAATCATAGCACCAGTCATCAGTGGCTGTGCGGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaact                                                                                                                                          >  2:90573/1‑136 (MQ=255)
          agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                                                                                                          >  2:198806/1‑136 (MQ=255)
           gTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGGTGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑aa                                                                                                                                        >  2:325926/1‑141 (MQ=255)
           gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCAACAa‑‑aa‑‑a                                                                                                                                     >  1:13125/1‑141 (MQ=255)
           gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                                                                                                         <  2:347035/141‑7 (MQ=255)
            tCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAGTAC‑‑aa‑‑a                                                                                                                                     >  1:26459/1‑141 (MQ=255)
            tCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑aa‑‑a                                                                                                                                     >  2:314793/1‑141 (MQ=255)
             cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat                                                                                                                                       >  1:580713/1‑133 (MQ=255)
              gCCTGTCTCTCAGGCGTCAATT‑‑GGCGAGGAGAGGAGGTGCATTACTTTGCTAATCCTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                                                      <  2:133413/141‑10 (MQ=255)
              gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                                                      <  2:568671/141‑10 (MQ=255)
              gCCTATCTCTCAGGCGTAAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                                                      <  2:641506/141‑10 (MQ=255)
               ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                                                                                                                     >  1:311927/1‑131 (MQ=255)
               ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                                                                                                                     <  1:814002/141‑11 (MQ=255)
               ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCCTAGTAAACGCaccaaataatta                                                                                                                                     >  1:191735/1‑130 (MQ=255)
                cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAGTCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                                                                    <  2:571334/141‑12 (MQ=255)
                cTAGCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                                                                    <  1:775842/141‑12 (MQ=255)
                atATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTACTAAACGCaccaaataattaa                                                                                                                                    >  2:514129/2‑129 (MQ=255)
                 tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                                                   <  2:59972/141‑13 (MQ=255)
                 tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaaa                                                                                                                                   >  2:284760/1‑129 (MQ=255)
                 tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGATGTGCTGTTGCTTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                                                   >  1:781601/1‑129 (MQ=255)
                 tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                                                   <  2:42152/141‑13 (MQ=255)
                 tATCTCTCAGGCGTCAA‑‑TAGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                                                   <  1:735220/141‑13 (MQ=255)
                  ttctctcAGGCGTCAATTTAGACGCAGAGCGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGt                                                                                                                               <  1:230311/140‑1 (MQ=255)
                  aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTCATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                                                  >  2:131690/1‑128 (MQ=255)
                  aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                                                  >  2:791499/1‑128 (MQ=255)
                  aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                                                  >  2:198995/1‑128 (MQ=255)
                  aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTAGTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                                                  >  2:751092/1‑128 (MQ=255)
                  aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                                                  >  1:417073/1‑128 (MQ=255)
                    ctctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAAGTGATGCTTATCTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCCATCC‑‑AA‑‑AATACGTgg                                                                                                                             >  2:139871/1‑140 (MQ=255)
                        cAGGCTTGAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACata                                                                                                                        <  2:663794/141‑1 (MQ=255)
                        cAGGCGTTGATTTAGACGAAGAATTTAGTTTTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACata                                                                                                                        <  2:625527/141‑1 (MQ=255)
                          ggCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTCTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGAGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACAtata                                                                                                                     >  1:788125/1‑141 (MQ=255)
                              tCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCg                                                                                                                <  2:122296/141‑1 (MQ=255)
                                   ttAGCCGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCa                                                                                                           <  1:331499/141‑1 (MQ=255)
                                   ttAGCCGCAGAGAGGAGGTGTACAAGGTGCTGCTTATTTTCGAAAGAGTCATCGTTTATGGCCCCGATCGCTGTGCCATTGCGTTTTGTGCTAGTTTTCCAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCa                                                                                                           <  2:804414/141‑1 (MQ=255)
                                    tggACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCaa                                                                                                          <  1:423127/139‑1 (MQ=255)
                                         gCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATaaaa                                                                                                     <  2:515932/141‑1 (MQ=255)
                                         gCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATaaaa                                                                                                     <  2:576682/141‑1 (MQ=255)
                                                ggaggTGTATAAGGTGATGCTGATTTTCGTTCGCATCATAGAACCAGTAATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCt                                                                                              <  2:610102/141‑1 (MQ=255)
                                                  aggTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTTTTTTTTGTTTAGTAGACGCAATAC‑‑AA‑‑AATGGGTGACATATACCCGCGCAAGTACAAATCCCCTcc                                                                                            >  2:787829/1‑140 (MQ=255)
                                                   tgtgtATAACGTTATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCAc                                                                                           <  1:52103/140‑1 (MQ=255)
                                                   ggTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCAc                                                                                           <  1:347823/141‑1 (MQ=255)
                                                       tataAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTAcc                                                                                       <  2:83127/141‑1 (MQ=255)
                                                           aaGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTCCCGCa                                                                                    >  2:22551/1‑141 (MQ=255)
                                                                    aTGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGt                                                                              >  2:522931/1‑141 (MQ=255)
                                                                        ttATTTTGGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgg                                                                          <  2:144275/141‑1 (MQ=255)
                                                                        ttATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgg                                                                          <  2:686996/141‑1 (MQ=255)
                                                                         tagtttCGTTCACATCATAGCACCAGTCATCAGTGGTTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAggg                                                                         <  2:136150/138‑1 (MQ=255)
                                                                         tATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAggg                                                                         >  2:229317/1‑141 (MQ=255)
                                                                          aTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCCATAAAAATCCCCTCACTACCGGAAATAGTGAgggg                                                                        <  2:720683/141‑1 (MQ=255)
                                                                           ttttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATGGTGAggggg                                                                       >  1:669930/1‑140 (MQ=255)
                                                                           ttttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                                                       <  1:314793/141‑1 (MQ=255)
                                                                           ttttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                                                       <  2:292295/141‑1 (MQ=255)
                                                                               cGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGGAGACGAAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCACCTCACTACCGCAAATAGTTAGGGGATTgg                                                                   >  2:91341/1‑141 (MQ=255)
                                                                                      atcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATaa                                                            <  2:713905/141‑1 (MQ=255)
                                                                                       tcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTATTTTTCTTATTGGCTTAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCCCACCAATAAAAATCCCCTCACTACCCACAATAGCGAGGGGATTGGTTTACaaa                                                           >  2:146387/1‑140 (MQ=255)
                                                                                         atAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGcat                                                              >  1:26559/1‑141 (MQ=255)
                                                                                            gCACCAGTCGGCACTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAAt                                                      <  1:603952/141‑1 (MQ=255)
                                                                                             cACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCGCCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATa                                                     >  1:653156/1‑141 (MQ=255)
                                                                                             cACCAGTCATCAGTGCCTGTGCCATTGCGTTTTTTTCATATTGCCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGAGATTGGTGTATAAGTAAATa                                                     <  1:589259/141‑1 (MQ=255)
                                                                                             cACCAATCATCACGGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTAAAGCAATCGATTATGTTATTGGTGTATAAGTAAATa                                                     <  1:405245/141‑1 (MQ=255)
                                                                                                cAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACtt                                                  <  2:658926/141‑1 (MQ=255)
                                                                                                 aGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTa                                                 >  1:79512/1‑141 (MQ=255)
                                                                                                 aGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGGGAGGGGAGTGGTGTATAAGTAAATACttt                                                 >  1:548900/1‑140 (MQ=255)
                                                                                                    catcaGTTGGTGTTCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAttt                                              <  2:564504/141‑1 (MQ=255)
                                                                                                     atcaGTGGCTGTGCCAGTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATGGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTTTAAGTAAATACTTAtttt                                             >  1:393119/1‑141 (MQ=255)
                                                                                                          tGGCTGTGCCATTGCGTTTTTTTTTTTTTTTTTAAGTAGACGCAATAC‑‑AA‑‑AGTAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTTCCGCAAATAGTGAGGGGGTTGGTGTATAAGTACATACTTTTTTTCGttt                                        >  2:572003/1‑140 (MQ=255)
                                                                                                              tgtgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTg                                        >  2:260277/1‑137 (MQ=255)
                                                                                                              tgtgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTg                                        <  1:260277/137‑1 (MQ=255)
                                                                                                              tgtgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCtt                                    <  1:385095/141‑1 (MQ=255)
                                                                                                               gtgCCATTGCGTTTCTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTa                                   >  1:546464/1‑141 (MQ=255)
                                                                                                                       gCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACt                           >  2:595940/1‑141 (MQ=255)
                                                                                                                            tttttCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTaa                      <  1:229317/141‑1 (MQ=255)
                                                                                                                            tttttCTTACTGGCTAAGTACACGCAATAC‑‑AA‑‑CATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGTAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTaa                      <  1:139871/141‑1 (MQ=255)
                                                                                                                                    aTTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCtt              >  1:513588/1‑141 (MQ=255)
                                                                                                                                          tAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTa        <  2:546464/141‑1 (MQ=255)
                                                                                                                                                gACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAAatat  <  2:79512/141‑1 (MQ=255)
                                                                                                                                                gACGCAATAA‑‑CA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAAatat  <  2:334782/141‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
ATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAAATAT  >  CP000730/1929806‑1930086

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: