Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,929,946 G→A 31.5% 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,9460GA31.5% 111.7 / 63.8 97intergenic (‑1029/‑139)USA300HOU_1794/USA300HOU_1796hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (33/33);  new base A (21/10);  total (54/43)
Fisher's exact test for biased strand distribution p-value = 1.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.56e-01

TAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATT  >  CP000730/1929807‑1930068
                                                                                                                                              |                                                                                                                                 
taataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGAT‑GCTTATTTTCGTTCACATAATAGCACCAGTTATCAGTGGCTGTGCCAATGCATTTTTTTCTTATTGGCTAAGTAg                                                                                                                                   <  1:28062/140‑1 (MQ=255)
 aataatAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAc                                                                                                                                 >  1:667840/1‑140 (MQ=255)
  ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCT‑TTA‑GTTCTTAAAGCAATTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACg                                                                                                                                <  1:154182/140‑1 (MQ=255)
  ataCTAGAGTAGCCTATCTCTCAGGCAAAAATT‑‑GAAGAAGAGAGGAGGTGAATTAA‑TTTGAA‑AAGCATTTTCAGTTGCATCAGATCACCAGTCATCAGTGGCTGTGGTGTTGCGTATGATACTTATTGGCTTAGTAAACg                                                                                                                                <  2:657606/140‑1 (MQ=255)
    agtagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGAT‑GCTTACTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGc                                                                                                                               <  1:116665/140‑1 (MQ=255)
     atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGACGAGGTGCATTAC‑TTTGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                                                                                             >  2:433154/2‑140 (MQ=255)
      tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaac                                                                                                                            >  2:353430/1‑139 (MQ=255)
         agacGCCTATCTCTCAGGCGACAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                                                                                                         <  1:179096/137‑5 (MQ=255)
         agTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGAT‑GCTGATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCATTTTTTTCTTATTGGCTAAGTAGACGCAATAc                                                                                                                          >  2:371126/1‑140 (MQ=255)
         agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                                                                                                         >  1:644172/1‑136 (MQ=255)
         agTCGCCTACCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTTATCAGTGGCTGTGTCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAc                                                                                                                          <  1:173036/140‑1 (MQ=255)
          gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                                                                                        >  1:631400/1‑135 (MQ=255)
          gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTTTGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                                                                                        >  1:438656/1‑135 (MQ=255)
          gTCGCCTATCTCTCAGACGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTTTGCT‑AATCATTTTCGATCACATCATAGCACCAGTCATCAGTGGCCGTGCTGTTGCGTATTATCCTTATTGGCTTAGTAAACGCaaaaaat                                                                                                                        >  2:188322/1‑135 (MQ=255)
          gTCGCCTATATCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTTTGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                                                                                        >  2:251746/1‑135 (MQ=255)
           tCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑aa                                                                                                                      >  1:284975/1‑140 (MQ=255)
           tCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCT‑TCTCAATTGCCGCCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                                                                                       <  1:484475/140‑7 (MQ=255)
           tCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTTTGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                                                                                       >  2:540239/1‑134 (MQ=255)
             gCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGTCTTTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACTCAATAC‑‑aa‑‑aa                                                                                                                  >  1:467819/1‑140 (MQ=255)
             gCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑aa‑‑aa                                                                                                                  <  2:177297/140‑1 (MQ=255)
             gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTTTGCT‑AATAATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat                                                                                                                     >  1:186500/1‑132 (MQ=255)
              cctctctttacggCGTAAATTGAGACCCAGAGAGGAGGTGTA‑TAA‑GGTGTT‑GGTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGAAGCAATAC‑‑AA‑‑AAt                                                                                                                 <  2:179563/129‑1 (MQ=255)
              ccTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGAT‑GCTTCTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCATATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AAt                                                                                                                 >  1:264990/1‑140 (MQ=255)
              ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                                    >  1:53523/1‑131 (MQ=255)
              ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                                    >  1:30862/1‑131 (MQ=255)
              ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                                    >  2:324772/1‑131 (MQ=255)
               cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                                                                                                   >  2:100951/1‑130 (MQ=255)
               cTAGCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCCGTCGATGTGCTGTTGCGTATTATACGTATTGGCTTAGTAAACGCaacaaataatta                                                                                                                   <  2:304254/140‑11 (MQ=255)
                tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTCC‑TTTGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGATGTGCTGTTGCGTATTATACTTATTAGCTTAGTAAACACaacaagtaattta                                                                                                                  >  1:649359/1‑129 (MQ=255)
                tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGATGTTGCGTATAATACTTATTGGCTTGGTAAACGCaacaaataattaa                                                                                                                  >  1:451416/1‑129 (MQ=255)
                tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCCGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                                                  >  1:418785/1‑129 (MQ=255)
                tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTTTGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                                                  >  1:307337/1‑129 (MQ=255)
                 aTCTCTCATGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                                 <  2:214678/140‑13 (MQ=255)
                 aTCTCTCAGTCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑CTTGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                                 <  2:532574/140‑13 (MQ=255)
                 aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattact                                                                                                                 >  2:624595/1‑128 (MQ=255)
                 aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                                 >  1:314945/1‑128 (MQ=255)
                   ctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaatag                                                                                                               <  2:660952/140‑15 (MQ=255)
                   ctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaatag                                                                                                               >  2:607649/1‑126 (MQ=255)
                        aGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ata                                                                                                      >  2:53056/1‑140 (MQ=255)
                             tccaTTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGcc                                                                                                <  2:654316/137‑1 (MQ=255)
                               aaTTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑Gc                                                                                             >  2:298503/1‑140 (MQ=255)
                                aattgacgcctgagagGAGGTGTA‑TAA‑GGTGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCa                                                                                            <  1:45457/129‑1 (MQ=255)
                                aTTTAGGCGCAGAGAGGAGGTGTA‑TAA‑GGTGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCa                                                                                            <  2:117695/140‑1 (MQ=255)
                                  ttAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGAT‑GCTTATTTTCGTTTACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCAcc                                                                                          >  2:355254/1‑140 (MQ=255)
                                        gCAGATACGAGTTGTA‑TAA‑GGTGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATCGGCTAAGTCGCCGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATaaa                                                                                    <  1:657444/140‑1 (MQ=255)
                                          agagagGAGGTGTA‑TAA‑GGTGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATaaaaa                                                                                  <  1:546191/140‑1 (MQ=255)
                                           gagagGAGGTGTA‑TAA‑GGTGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GAACCAATAAAAAt                                                                                 >  1:373733/1‑140 (MQ=255)
                                            agagGAGGTGTA‑TAA‑GGTGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATc                                                                                <  1:581282/140‑1 (MQ=255)
                                            agagGAGGTGTA‑TAA‑GGTGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATc                                                                                >  1:442154/1‑140 (MQ=255)
                                              gggaggGGTA‑TAA‑GGTGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCCAAAAAAATcac                                                                              >  2:541456/2‑138 (MQ=255)
                                              aggaggTGTA‑TAA‑GGTGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCCATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCAAACATAAAAATccc                                                                              >  2:228734/1‑140 (MQ=255)
                                               ggaggTGTA‑TAA‑GGTGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATcccc                                                                             <  1:420496/140‑1 (MQ=255)
                                                gaggTGTA‑TAA‑GGTGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCt                                                                            <  2:375683/140‑1 (MQ=255)
                                                  ggTGTA‑TAA‑GGTGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCCTTGCGTTTTTTTCTTATTGGCTAAGTAGAGGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCACTAAAAATCCCCCCa                                                                          >  2:293383/1‑140 (MQ=255)
                                                   gtgtA‑TAA‑GGTGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCAc                                                                         >  1:556422/1‑140 (MQ=255)
                                                    tgtA‑TAA‑GGTGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACt                                                                        <  1:208517/140‑1 (MQ=255)
                                                     gtA‑TAA‑GGTGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTa                                                                       >  1:634105/1‑140 (MQ=255)
                                                     gtA‑TAA‑GGTGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCAATTGCGTTTTTTTCTTATTGGCTAAGTAGAAGAAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCACCTCACTa                                                                       >  1:433186/1‑140 (MQ=255)
                                                     gtA‑TAA‑GGTGAT‑GCTTATTCTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGAAATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTa                                                                       >  1:354392/1‑140 (MQ=255)
                                                          aa‑GGTGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGc                                                                   <  1:359776/140‑1 (MQ=255)
                                                           a‑GGTGAT‑GCTTATTTTCTTTCACATCATAGCACCAGTCATCAGTGGCTGGGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑CCCCGGTCAC‑AGAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCa                                                                  <  2:607591/140‑1 (MQ=255)
                                                              gTGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCaaa                                                                >  2:255216/1‑140 (MQ=255)
                                                              gTGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTACGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCaaa                                                                >  1:17989/1‑140 (MQ=255)
                                                               tgtt‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTatttgcggt                                                               <  1:601680/137‑9 (MQ=255)
                                                                taT‑GCTTATTTTTGTTCACATCATAGCACCAGTCATCAGTGGCTTTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATa                                                              <  1:301929/139‑1 (MQ=255)
                                                                 cacggctATTTTCGTTCACATCATAGCGCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATa                                                              <  2:666234/134‑1 (MQ=255)
                                                                     cTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGa                                                          >  1:60442/1‑140 (MQ=255)
                                                                      ttATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGCC‑ATAT‑AGCC‑GCACCAATACACAACCCCTCACTACCGCAAGTaacgcg                                                         >  2:552250/1‑135 (MQ=255)
                                                                      ttATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATCGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATTAAAATCCCCTCACTACCACAAagagagag                                                         >  2:424809/1‑140 (MQ=255)
                                                                       tATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgg                                                        <  2:58478/140‑1 (MQ=255)
                                                                       tATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGGTTTGTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgg                                                        >  1:457144/1‑140 (MQ=255)
                                                                       tATTTCCGTTCCTATTATAGCCCCAGTCATCAGTGGCTGTGCCATTTCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgg                                                        <  2:355065/140‑1 (MQ=255)
                                                                         ttttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgggg                                                      <  1:130121/140‑1 (MQ=255)
                                                                         ttttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgggg                                                      <  1:290266/140‑1 (MQ=255)
                                                                               ttCACATCATAGCACCAATCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGgt                                                >  2:25765/1‑140 (MQ=255)
                                                                                     tcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCCAAAAGTGAGGGGATTGGTGTATaa                                          <  2:643649/140‑1 (MQ=255)
                                                                                     tcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGAATaa                                          >  2:665917/1‑140 (MQ=255)
                                                                                      catAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAg                                         <  1:146383/140‑1 (MQ=255)
                                                                                        tAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGA‑AAAAATCCCCACACTACTGC‑AATAGTGTGGGGATTGgtgcat                                            <  1:102347/140‑4 (MQ=255)
                                                                                           cACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCAACAATAAAAATCCCCTCACTACCGCAAATAATGAGGGGATTGGTGTATAAGTAAAt                                    >  1:196179/1‑140 (MQ=255)
                                                                                            aCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATa                                   >  1:185855/1‑140 (MQ=255)
                                                                                                gTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTa                               >  1:562388/1‑140 (MQ=255)
                                                                                                  catcaGTGGCTGCGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAT‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTAGCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAtt                             <  1:53056/140‑1 (MQ=255)
                                                                                                   atcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAttt                            <  1:240119/140‑1 (MQ=255)
                                                                                                   atcaGTGGCTGTGCAATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTATCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAttt                            <  1:343325/140‑1 (MQ=255)
                                                                                                    tcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTATGTAAACGTAATAC‑‑TA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATAC‑TAttttt                          <  2:467819/140‑2 (MQ=255)
                                                                                                     caGTGGCTGTGCCATTGCGTTTTTTTTTTGTTGGGTGGGAAGACGCAATAC‑‑AA‑‑AAGAGGAGAC‑GTAT‑AGCC‑GCACCAATGAAAATCCCCTCACTACCGCAAATAGTAGGGGGATTGGTGTATAAGGAAATACTTTTTTTc                          >  2:341920/1‑140 (MQ=255)
                                                                                                     caGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑CA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTc                          >  1:221683/1‑140 (MQ=255)
                                                                                                            tgtgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCt                   >  2:374244/1‑140 (MQ=255)
                                                                                                            tgtgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATACGGAGGGGGTAGGTGTATAAGTAAATACTTTTTTTTGTTTTCt                   >  2:370638/1‑140 (MQ=255)
                                                                                                              tgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAACACTTATTTTCGTTGTCTTa                 <  1:371126/140‑1 (MQ=255)
                                                                                                               gCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTaa                <  1:637807/140‑1 (MQ=255)
                                                                                                                 cATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGAGTATAAGTAAATACTTATTTTCGTTGTCTTAAtt              >  1:611717/1‑140 (MQ=255)
                                                                                                                 cATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAAAAAAAAACCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAAAACTTATTTTCGTTGTCTTAAtt              >  1:596654/1‑140 (MQ=255)
                                                                                                                  aTTGCGTTTTTTTCGTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTTCCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATta             <  1:180334/140‑1 (MQ=255)
                                                                                                                     gCGTTTTTTTCTCTTTGGCTAAGTAGACTCAATAC‑‑AA‑‑AATAGGGCAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTTCCGTAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTTTCTTAATTATAc          <  2:43113/140‑1 (MQ=255)
                                                                                                                             ttCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑ATTAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAATTTTTTGCTTTGTTTTATACTGCTAAtt  <  2:642826/140‑1 (MQ=255)
                                                                                                                             ttCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAAtt  <  1:25765/140‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                 
TAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAC‑‑AA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATT  >  CP000730/1929807‑1930068

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