Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 940,348 T→C 5.5% R413R (CGT→CGC USA300HOU_0903 → NhaC family sodium:proton (Na+:H+) antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730940,3480TC 5.5% 259.8 / 4.6 90R413R (CGT→CGCUSA300HOU_0903NhaC family sodium:proton (Na+:H+) antiporter
Reads supporting (aligned to +/- strand):  ref base T (44/41);  new base C (3/2);  total (47/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

CAATTGGACTAGATACAATGGCATTAATCGCATTGATTGGAACAGCGAGTGCATTAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGACTCATTTATTAGAGACATTTGAGATGTCAATAGATCAC  >  CP000730/940209‑940486
                                                                                                                                           |                                                                                                                                          
cAATTGGACTAGATACAATGGCATTAATCGCATTGATTGGAACAGCGAGTGCATTAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTg                                                                                                                                           <  1:469742/141‑1 (MQ=255)
      gACTAGATACAATGGCATTAATCGCATTGATTGGAACAGCGAGTGCATTAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACAt                                                                                                                                     >  2:617237/1‑141 (MQ=255)
       aCTAGATACAATGGCATTAATCGCATTGATTGGAACAGCGAGTGCATTAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACAtg                                                                                                                                    >  1:444752/1‑141 (MQ=255)
        cTAGATACAATGGCATTAATCGCATTGATTGGAACAGCGAGTGCATTAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACAtgt                                                                                                                                   >  2:30268/1‑141 (MQ=255)
         tAGATACAATGGCATTAATCGCATTGATTGGAACAGCGAGTGCATTAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACAtgtg                                                                                                                                  >  2:22649/1‑141 (MQ=255)
          aGATACAATGGCATTAATCGCATTGATTGGAACAGCGAGTGCATTAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATgtgt                                                                                                                                 <  2:599951/141‑1 (MQ=255)
           gATACAATGGCATTAACCGTCCAGATGGGAACAGCGACTTCATTAGATTACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTa                                                                                                                                <  1:584961/141‑1 (MQ=255)
             tACAATGGCATTAATCGCATTGATTGGAACAGCGAGTGCATTAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTAcc                                                                                                                              >  1:588250/1‑141 (MQ=255)
              aatattgtaattaTCGCATTGATTGGAACAGCGAGTGCATTAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCa                                                                                                                             <  1:600126/129‑1 (MQ=255)
               cAATGGCATTAATCGCATTGATTGGAACAGCGAGTGCATTAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCaa                                                                                                                            >  1:831400/1‑141 (MQ=255)
                  tGGCATTAATCGCATTGATTGGAACAGCGAGTGCATTAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAAct                                                                                                                         <  1:807087/141‑1 (MQ=255)
                   ggCATTAATCGCATTGATTGGAACAGCGAGTGCATTAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAActt                                                                                                                        <  2:822429/141‑1 (MQ=255)
                    gCATTAATCGCATTGATTGGAACAGCGAGTGCATTAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAActtc                                                                                                                       >  2:738473/1‑141 (MQ=255)
                    gCATTAATCGCATTGATTGGAACAGCGAGTGCATTAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAActtc                                                                                                                       >  2:462188/1‑141 (MQ=255)
                      aTTAATCGCATTGATTGGAACAGCGAGTGCATTAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAActtctt                                                                                                                     <  2:1045124/141‑1 (MQ=255)
                       ttAATCGCATTGATTGGAACAGCGAGTGCATTAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTg                                                                                                                    <  2:671911/141‑1 (MQ=255)
                         aaTCGCATTGATTGGAACAGCGAGTGCATTAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGtt                                                                                                                  >  1:272016/1‑141 (MQ=255)
                           tCGCATTGATTGGAACAGCGAGTGCATTAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGtttt                                                                                                                <  1:738473/141‑1 (MQ=255)
                                 tgattgGAACAGCGAGTGCATTAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTataata                                                                                                          >  1:570150/1‑141 (MQ=255)
                                 tgagtgGAACAGCGAGTGCATTAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATCAAAACGTAACTTCTTGTTTTataata                                                                                                          <  1:526506/141‑1 (MQ=255)
                                  gattgGAACAGCGAGTGCATTAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAatat                                                                                                         <  1:692589/141‑1 (MQ=255)
                                    ttgGAACAGCGAGTGCATTAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTc                                                                                                       >  1:353772/1‑141 (MQ=255)
                                      ggAACAGCGAGTGCATTAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCt                                                                                                     <  2:673378/141‑1 (MQ=255)
                                      ggAACAGCGAGTGCATTAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCt                                                                                                     >  1:385877/1‑141 (MQ=255)
                                         aCAGCGAGTGCATTAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTa                                                                                                  >  2:939699/1‑141 (MQ=255)
                                         aCAGCGAGTGCATTAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTa                                                                                                  >  2:3766/1‑141 (MQ=255)
                                          cAGCGAGTGCATTAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTaa                                                                                                 >  2:856761/1‑141 (MQ=255)
                                          cAGCGAGTGCATTAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTaa                                                                                                 >  1:1044817/1‑141 (MQ=255)
                                           aGCGAGTGCATTAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAAt                                                                                                >  1:534060/1‑141 (MQ=255)
                                                gcgCATTAGGTGACTCATGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATAAGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGAttt                                                                                           >  1:721370/3‑141 (MQ=255)
                                                   cATTAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCg                                                                                        <  2:926605/141‑1 (MQ=255)
                                                     ttAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGt                                                                                      <  2:671169/141‑1 (MQ=255)
                                                        ggTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACt                                                                                   >  1:182262/1‑141 (MQ=255)
                                                          tGACTCAGGTTCGCCTGCAAGTGATTCAACACTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTAt                                                                                 <  2:15421/141‑1 (MQ=255)
                                                             cTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATtgc                                                                              >  2:968749/1‑141 (MQ=255)
                                                             cTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATtgc                                                                              <  1:340429/141‑1 (MQ=255)
                                                              tCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATtgct                                                                             >  1:843579/1‑141 (MQ=255)
                                                              tCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACGTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATtgct                                                                             <  1:855154/141‑1 (MQ=255)
                                                              tCAGGTTCGCCTGAATGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATtgct                                                                             <  2:984453/141‑1 (MQ=255)
                                                                  gTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTa                                                                         >  2:931132/1‑141 (MQ=255)
                                                                    tCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATg                                                                       <  2:160144/141‑1 (MQ=255)
                                                                      gCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGt                                                                     <  2:230929/141‑1 (MQ=255)
                                                                      gCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGt                                                                     <  2:339804/141‑1 (MQ=255)
                                                                        cTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTAc                                                                   >  2:808157/1‑141 (MQ=255)
                                                                           cAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGCGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACtat                                                                >  2:165235/1‑141 (MQ=255)
                                                                            aaGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACtata                                                               <  1:833658/141‑1 (MQ=255)
                                                                              gTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAg                                                             <  1:715948/141‑1 (MQ=255)
                                                                               tGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGc                                                            <  2:746501/141‑1 (MQ=255)
                                                                                gATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCa                                                           <  1:3766/141‑1 (MQ=255)
                                                                                    cAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATaaa                                                       >  1:138464/1‑141 (MQ=255)
                                                                                    cAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATaaa                                                       >  2:818185/1‑141 (MQ=255)
                                                                                          tAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTg                                                 >  2:196146/1‑141 (MQ=255)
                                                                                          tAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTg                                                 <  1:1034344/141‑1 (MQ=255)
                                                                                            ggACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGaa                                               >  1:1016590/1‑141 (MQ=255)
                                                                                              aCCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGaaaa                                             >  1:765805/1‑141 (MQ=255)
                                                                                               ccAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAat                                            >  2:1007251/1‑141 (MQ=255)
                                                                                                  aCTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATg                                         >  2:982317/1‑141 (MQ=255)
                                                                                                   cTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTTGGTCTTATTGATGCTATGGCATTATAGGCATCTTGAGGTGAAAATATGa                                        >  1:98838/1‑141 (MQ=255)
                                                                                                   cTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGa                                        <  2:881861/141‑1 (MQ=255)
                                                                                                    tGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGAc                                       <  1:574755/141‑1 (MQ=255)
                                                                                                    tGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATAATTTGCGGTACTACTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGAc                                       <  2:779187/141‑1 (MQ=255)
                                                                                                     gCGGTATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGACt                                      <  2:735427/141‑1 (MQ=255)
                                                                                                     gCGGGATTAAATGTTGATGGCCAACATGATCATATACGCGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGACt                                      >  1:823543/1‑141 (MQ=255)
                                                                                                     gCGGGATGAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGATACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGACt                                      <  2:380847/141‑1 (MQ=255)
                                                                                                      cGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGACTc                                     <  2:593279/141‑1 (MQ=255)
                                                                                                       gggATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGACTCa                                    >  2:26726/1‑141 (MQ=255)
                                                                                                       gggATTAAATGTTGATGGCCAACATGATCATATACGCGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGACTCa                                    <  2:241537/141‑1 (MQ=255)
                                                                                                       gggATTAAATGTTGATGGCCAACATGATCATATACGCGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGACTCa                                    <  2:241536/141‑1 (MQ=255)
                                                                                                        ggATTAAATGTTGATGGCCAACATGATCATATACGCGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGACTCAt                                   >  1:857111/1‑141 (MQ=255)
                                                                                                         gATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGACTCAtt                                  >  1:867129/1‑141 (MQ=255)
                                                                                                          aTTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGACTCAttt                                 <  1:730975/141‑1 (MQ=255)
                                                                                                            tAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGACTCATttat                               <  1:295144/141‑1 (MQ=255)
                                                                                                             aaaTGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTCTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGACTCATttctt                              >  2:1001997/1‑141 (MQ=255)
                                                                                                               aTGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGACTCATTTATTag                            <  1:65244/141‑1 (MQ=255)
                                                                                                               aTGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGACTCATTTATTag                            <  1:255117/141‑1 (MQ=255)
                                                                                                                tGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGACTCATTTATTaga                           <  1:495151/141‑1 (MQ=255)
                                                                                                                 gTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGACTCATTTATTagag                          >  1:93952/1‑141 (MQ=255)
                                                                                                                   tGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGACTCATTTATTAGAGAc                        <  2:71456/141‑1 (MQ=255)
                                                                                                                   tGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGACTCATTTATTAGAGAc                        <  1:975520/141‑1 (MQ=255)
                                                                                                                    gATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGACTCATTTATTAGAGACa                       >  2:140942/1‑141 (MQ=255)
                                                                                                                       ggCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGACTCATTTATTAGAGACAttt                    >  2:128224/1‑141 (MQ=255)
                                                                                                                       ggCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGACTCATTTATTAGAGACAttt                    >  2:876396/1‑141 (MQ=255)
                                                                                                                         ccAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGACTCATTTATTAGAGACATTTga                  <  2:515530/141‑1 (MQ=255)
                                                                                                                          cAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGACTCAGTTATTAGAGACATTTgag                 >  1:39562/1‑141 (MQ=255)
                                                                                                                           aaCATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGACTCATTTATTAGAGACATTTgaga                >  2:560867/1‑141 (MQ=255)
                                                                                                                              aTGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGACTCATTTATTAGAGACATTTGAGATGt             >  2:927054/1‑141 (MQ=255)
                                                                                                                                   cATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGACTCATTTATCAGAGACATTTGAGATGTCAATa        >  1:443949/1‑141 (MQ=255)
                                                                                                                                    atatACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGACTCATTTATTAGAGACATTTGAGATGTCAATAg       >  2:343312/1‑141 (MQ=255)
                                                                                                                                     tataCGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGACTCATTTATTAGAGACATTTGAGATGTCAATAGa      <  2:609226/141‑1 (MQ=255)
                                                                                                                                         cGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGACTCATTTATTAGAGACATTTGAGATGTCAATAGATCAc  >  2:255556/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                          
CAATTGGACTAGATACAATGGCATTAATCGCATTGATTGGAACAGCGAGTGCATTAGGTGACTCAGGTTCGCCTGCAAGTGATTCAACATTAGGACCAACTGCGGGATTAAATGTTGATGGCCAACATGATCATATACGTGATACATGTGTACCAAACTTCTTGTTTTATAATATTCCTTTAATGATTTTCGGTACTATTGCTGCTATGGTACTATAAGCATAAAGAGGTGAAAATATGACTCATTTATTAGAGACATTTGAGATGTCAATAGATCAC  >  CP000730/940209‑940486

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