Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,458,481 T→C 6.5% R32R (AGA→AGG)  USA300HOU_2319 ← aldose 1‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,458,4810TC 6.5% 204.4 / 2.9 73R32R (AGA→AGG) USA300HOU_2319aldose 1‑epimerase
Reads supporting (aligned to +/- strand):  ref base T (33/34);  new base C (2/3);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATTTAGCATTATCAATTCTACCGCTATATCTTCCTATAGAAGCACCAAATTTAAAAGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTGTCTTTGGTGCTCAACTTCAACTATCATTTTTACTTCTCCCTTCTAACCACAAGTGTTCAAGCTCTG  >  CP000730/2458343‑2458615
                                                                                                                                          |                                                                                                                                      
aaaTTTAGCATTATCAATTCTACCGCTATATCTTCCTATAGAAGCACCAAATTTAAAAGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTActct                                                                                                                                      >  1:692429/1‑141 (MQ=255)
    ttAGCATTATCAATTCTACCGCTATATCTTCCTATAGAAGCACCAAATTTAAAAGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTTTCATGATATTTCCAAGACACTACCCTTGCt                                                                                                                                  >  2:277867/1‑141 (MQ=255)
    ttAGCATTATCAATTCTACCGCTATATCTTCCTATAGAAGCACCAAATTTAAAAGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACCCTTGcc                                                                                                                                  >  2:680907/1‑140 (MQ=255)
      aGCATTAGTAATTCTACCGCTATATCTTCCTATAGAAGCACTAAATTAAAAAGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTcc                                                                                                                                <  2:645188/141‑1 (MQ=255)
        cATTATCAATTCTACCGCTATATCTTCCTATAGAAGCACCAAATTTAAAAGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCAt                                                                                                                              >  1:631466/1‑141 (MQ=255)
        cATTATCAATTCTACCGCTATATCTTCCTATAGAAGCACCAAATTTAAAAGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCAt                                                                                                                              >  1:71710/1‑141 (MQ=255)
        cATTATCAATTCCACCGCTATATCTTCCTATAGAAGCACCAAATTTAAAAGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCAt                                                                                                                              <  1:130184/141‑1 (MQ=255)
           tATCAATTCTACCGCTATATCTTCCTATAGAAGCACCAAATTTAAAAGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAAt                                                                                                                           <  2:136718/141‑1 (MQ=255)
           tATCAATTCCACCCCTATATCTTCCTATAGAAGCACCAAATTTAAAAGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAAt                                                                                                                           <  2:585057/141‑1 (MQ=255)
               aaTTCTACTGCTATATCTTCCTATAAAAGCACCAAATTTAAAAGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGt                                                                                                                       <  2:124340/141‑1 (MQ=255)
                 ttCTACCGCTATATCTTCCTATAGAAGCACCAAATTTAAAAGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTaa                                                                                                                     <  1:116401/141‑1 (MQ=255)
                   cTACCGCTATATCTTCCTATAGAAGCACCAAATTTAAAAGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTAACATGATATTTCCAAGACACTACTCTTGTTCCATAATTCGTaaaa                                                                                                                   >  2:416734/1‑141 (MQ=255)
                        gCTATATCTTCCTATAGAAGCACCAAATTTAAAAGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAaataat                                                                                                              >  1:158471/1‑141 (MQ=255)
                         cTATATCTTCCTATAGAAGCACCAAATTTAAAAGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAAtt                                                                                                             <  1:765625/141‑1 (MQ=255)
                         cTATATCTTCCTATAGAAGCACCAAATTTAAAAGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAAtt                                                                                                             <  1:260035/141‑1 (MQ=255)
                         cTATATCTTCCTATAGAAGCACCAAATTTAAAAGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCCCCATAATTCGTAAAAATAAtt                                                                                                             >  1:74454/1‑141 (MQ=255)
                               cTTCCTATAGAAGCACCAAATTTAAAAGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGtt                                                                                                       >  2:683574/1‑141 (MQ=255)
                                       aGAAGCACCAAATTTAAAAGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCAtt                                                                                               <  1:400901/141‑1 (MQ=255)
                                          aGCACCAAATTTAAAAGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATc                                                                                            >  1:455337/1‑141 (MQ=255)
                                               cAAATTTAAAAGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAAttt                                                                                       >  1:658707/1‑141 (MQ=255)
                                                aaaTTTAAAAGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTCAAAATAATTTTAGTTTCATCATTATCAAtttt                                                                                      >  2:290949/1‑141 (MQ=255)
                                                aaaTTTAAAAGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACAATACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAAtttt                                                                                      >  1:184204/1‑141 (MQ=255)
                                                   tttAAAAGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGAt                                                                                   >  2:280245/1‑141 (MQ=255)
                                                    ttAAAAGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCCTGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAACAATTTTAGTTTCAAAATTATAAATTTTGAtt                                                                                  >  1:383953/1‑141 (MQ=255)
                                                    ttAAAAGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATATTAGGTTCATCCTTATCAATTTTGAtt                                                                                  >  2:165644/1‑141 (MQ=255)
                                                     tAAAAGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTa                                                                                 <  1:754938/141‑1 (MQ=255)
                                                     tAAAAGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTa                                                                                 <  1:722031/141‑1 (MQ=255)
                                                      aacaGGCTTGCTATGATAAAATTCATGCGGTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTaa                                                                                <  2:283933/138‑1 (MQ=255)
                                                      aaaaGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTaa                                                                                <  2:478362/141‑1 (MQ=255)
                                                      aaaaGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTaa                                                                                <  1:766851/141‑1 (MQ=255)
                                                      aaaaGGATTACTACGATAAAATTCATCCGCTTCAACAACATTCCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTaa                                                                                <  2:42721/141‑1 (MQ=255)
                                                        aaGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTaaaa                                                                              >  2:498962/1‑140 (MQ=255)
                                                        aaGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAAt                                                                              >  1:511909/1‑141 (MQ=255)
                                                        aaGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAAt                                                                              >  1:691559/1‑141 (MQ=255)
                                                         aGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTATTTTCATGATTATCAATTTTGATTCAATc                                                                             >  1:79617/1‑141 (MQ=255)
                                                               aCTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACAcc                                                                       <  1:527025/141‑1 (MQ=255)
                                                                    gATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTGTc                                                                  <  2:621690/141‑1 (MQ=255)
                                                                          aTTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTGTCTTTGGt                                                            >  2:440527/1‑141 (MQ=255)
                                                                          aTTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTGTCTTTGGt                                                            >  2:549396/1‑141 (MQ=255)
                                                                          aTTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTGTCTTTGGt                                                            >  2:798694/1‑141 (MQ=255)
                                                                          aTTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTGTCTTTGGt                                                            >  1:375762/1‑141 (MQ=255)
                                                                             cATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTGTCTTTGGTGCt                                                         <  1:683574/141‑1 (MQ=255)
                                                                              ctCCGCTTCAACAACATTTCCAAGAACAATAGTATTATCATGATATTTCCACGACTCTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATTGATACCAATCTTGATTAAATCTACACCTTGTCTTTGGTGCtc                                                        <  2:397343/140‑1 (MQ=255)
                                                                              aTCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTGTCTTTGGTGCtc                                                        <  1:339083/141‑1 (MQ=255)
                                                                                ccGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTGTCTTTGGTGCtcaa                                                      >  2:735228/1‑141 (MQ=255)
                                                                                ccGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTGTCTTTGGTGCtcaa                                                      >  1:44142/1‑141 (MQ=255)
                                                                                ccGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGAGTAAATCTACACCTTGTCTTTGGTGCtcaa                                                      >  2:11585/1‑141 (MQ=255)
                                                                                 cggttCCACAGCATTTGCAAGAACAAGATTATTATTGTGATATTTTCAAGACACTACTCTTGCTCCATAATTCGTTAATATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTGTCTTTGGTGCtcaac                                                     <  2:126477/138‑1 (MQ=255)
                                                                                      caacaaCATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACCCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTGTCTTTGGTGCtcaacttcaa                                                <  1:277867/141‑1 (MQ=255)
                                                                                      caacaaCATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACCCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTGTCTTTGGTGCtcaacttcaa                                                <  1:680907/141‑1 (MQ=255)
                                                                                       ccaagcaTTGCCAAGATTAATATAATTATCATTATATTTCCAAGACACTACTCGTGCTCCATAATTCGTAAAAATAATTGTATTTTCATCATTATCAAGTTTGATTAAATCTACACCTTGTCTTTGGTGCtcaacttcaac                                               <  1:566714/136‑1 (MQ=255)
                                                                                       aacaacATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTGTCTTTGGTGCtcaacttcaac                                               <  1:669318/141‑1 (MQ=255)
                                                                                       aacaacATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTCTGATTAAATCTACACCTTGTCTTTGGTGCtcaacttcaac                                               <  2:158471/141‑1 (MQ=255)
                                                                                        acaacaTTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTGTCTTTGGTGCtcaacttcaact                                              >  1:8846/1‑141 (MQ=255)
                                                                                        acaacaTTTCCAAGAACAAAATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTGTCTTTGGTGCtcaacttcaact                                              >  1:33573/1‑141 (MQ=255)
                                                                                               ttCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTGTCTTTGGTGCTCAAGTTCAACTATCAttt                                       <  1:443309/141‑1 (MQ=255)
                                                                                               ttCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTGTCTTTGGTGCTCAACTTCAACTATTAttt                                       <  1:552213/141‑1 (MQ=255)
                                                                                                cccAAGAACAATATTATTATCATGCTATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTGTCTTTGGTGCTCAACTTCAACTATCAtttt                                      <  1:409549/140‑1 (MQ=255)
                                                                                                cccAAGAAACATATTATAATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTGTCTTTGGTGCTCAACTTCAACTATCAtttt                                      <  1:409630/140‑1 (MQ=255)
                                                                                                   aaGAACAATATTCATATCATGATATTTCCACGACACTACTCTTGCTCCATAATTCGTAAAAATAATCTTAGTTTCATTACTATCAATTTAGATTAAACCTACACCTTGTCCTCGGTGCTCAACTTCAACTATCATTTTTAc                                   <  2:399963/141‑1 (MQ=255)
                                                                                                   aaGAAAAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTGTCTTTGGTGCTCAACTTCAACTATCATTTTTAc                                   >  1:439888/1‑141 (MQ=255)
                                                                                                      aaCAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTGTCTTTGGTGCTCAACTTCAACTATCATTTTTACTtc                                <  1:594359/141‑1 (MQ=255)
                                                                                                      aaCAATATTATTATCATGATATGTCCAAGAAACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTGTCTTTGGTGCTCAACTTCAACTATCATTTTTACTtc                                <  1:559851/141‑1 (MQ=255)
                                                                                                        cAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTACAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTGTCTTTGGTGCGCAACTTCAACTATCATTTTTACTtctc                              >  2:556714/1‑141 (MQ=255)
                                                                                                           tattattatCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTGTCTTTGGTGCTCAACTTCAACTATCATTTTTACTTCTCCCt                           <  2:692429/141‑1 (MQ=255)
                                                                                                              tattatCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAGATCTACACCTTGTCTTTGGTGCTCAACTTCAACTATCCTTTTTACTTCTCtcttct                        >  2:469921/1‑141 (MQ=255)
                                                                                                              tattatCATGATATTTCCAAGACACTACACTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTGTCTTTGGTGCTCAACTTCAACTATCATTTTTACTTCTCCCTTCt                        <  2:509035/141‑1 (MQ=255)
                                                                                                                   tCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTGTCTTTGGTGCTCAACTTCAACTATCATTTTTACTTCTCCCTTCTAACca                   <  2:455337/141‑1 (MQ=255)
                                                                                                                     aTGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTGTCTTTGGTGCTCAACTTCAACTATCATTTTTACTTCTCCCTTCTAACcaca                 >  1:178401/1‑141 (MQ=255)
                                                                                                                        atatTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTTTCTTTGGTGCCCAACTTCAACTATCATTTTTACTTCTCCCTTCTAACCCCAAgg              >  2:82908/1‑140 (MQ=255)
                                                                                                                          atTTCTAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTGTCTTTGGTGCTCAACTTCAACTATCATTTTTACTTCTCCCTTCTAACCACAAgtgt            <  1:312966/141‑1 (MQ=255)
                                                                                                                               cAAGACCCTACCCTTGCTCCATAATTCGTAAAAATAATTTTAGGTTCATCATTATCAATTTTGATTAAATCTACACCTTGTCTTTGTTGCTCAACTTCAACTATCATTTTTACTTCTCCCTTCTAGCCACAAGTGTTCAAg       <  1:525821/141‑1 (MQ=255)
                                                                                                                                   acacTACTCTTGCTCCATAATTCGTAAAAATAAATTTCGTTTCCTCATTATCAATTTGGATTAAATCTACACCTTGTCTTTGGTGCTAAACTTTAACTATCATGGTTCCTTTTCCCTTGTAACCGCAAGTGTTCAAGttct   >  2:191029/1‑141 (MQ=255)
                                                                                                                                    cacTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTGTCTTTGGTGCTCAACTTCAACTATCATTTTTACTTCTCCCTTCTAACCACAAGTGTTCAAGCTctg  >  1:723958/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                      
AAATTTAGCATTATCAATTCTACCGCTATATCTTCCTATAGAAGCACCAAATTTAAAAGGATTACTATGATAAAATTCATCCGCTTCAACAACATTTCCAAGAACAATATTATTATCATGATATTTCCAAGACACTACTCTTGCTCCATAATTCGTAAAAATAATTTTAGTTTCATCATTATCAATTTTGATTAAATCTACACCTTGTCTTTGGTGCTCAACTTCAACTATCATTTTTACTTCTCCCTTCTAACCACAAGTGTTCAAGCTCTG  >  CP000730/2458343‑2458615

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: ‑