Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 18,890 C→G 93.9% N182K (AAC→AAG USA300HOU_0014 → DHH subfamily phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073018,8900CG93.9% 238.4 / 5.3 77N182K (AAC→AAGUSA300HOU_0014DHH subfamily phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base G (45/27);  total (48/29)
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.01e-01

TATAATCAGTATCATTTCCAAGTACGTTATTCTGAGAATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGC  >  CP000730/18753‑19024
                                                                                                                                         |                                                                                                                                      
tataATCCGTATCATTTCCAAGTACGTTATTCTGAGAATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTAc                                                                                                                                     >  2:679492/1‑141 (MQ=255)
 ataATCAGCATCATTTTCAAGTACGTTATTCCGAGAATGATCACGTCCTCTATTTCTTTGATATAACTGAACAAGTACCAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACAATCTTTTGAGAGAAGTACg                                                                                                                                    <  2:189360/141‑1 (MQ=255)
  taATCAGTATCATTTCCAAGTACGTTATTCTGAGAATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAGAATTCTAAACCAATCATTACGAAATTATTTTTAGATAAGTACGa                                                                                                                                   >  2:443155/1‑141 (MQ=255)
  aaagcAGTATCATTTCCAAGTACGTTATTCTGAGAATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGa                                                                                                                                   <  2:13394/137‑1 (MQ=255)
  aaaTCAGTATCATTTCCAAGTACGTTATTCTGAGAATGATCACTGCCTCTATCTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGa                                                                                                                                   <  2:13340/140‑1 (MQ=255)
    aTCAGTATCATTTCTCAGTTCGGTATTCTGAGAATGATCACTGCATCTTTTTTCTTGATATAACTGAACAAGTACAAACAAATGAAGTTTATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATg                                                                                                                                 <  2:181102/141‑1 (MQ=255)
     tCAGTATCATTTCCAAGTACGTTATTCTGAGAATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATga                                                                                                                                <  1:790733/141‑1 (MQ=255)
     tCAGTATCATTTCCAAGTACGTTATTCTGAGAATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATga                                                                                                                                <  1:95698/141‑1 (MQ=255)
       aGTATCATTTCCAAGTACGTTATTCTGAGAATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGCCATTATTTTTAGATAAGTACGATgaga                                                                                                                              >  1:377450/1‑141 (MQ=255)
       aGTATCATTTCCAAGTACGTTATTCTGAGAATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATgaga                                                                                                                              >  1:445723/1‑141 (MQ=255)
              tttCCAAGTACGTTATTCTGAGGCTAAGCACGTCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCa                                                                                                                       >  1:206983/1‑141 (MQ=255)
              tttCCAAGTACGTTATTCTGAGAATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAGAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCa                                                                                                                       >  2:851551/1‑141 (MQ=255)
              tttCCAAGTACGTTATTCTGAGAATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTACACCCATAATCACGACATTATCTTTAGCTAACTACGATGAAATTACGCa                                                                                                                       >  2:590101/1‑141 (MQ=255)
              tttCCAAGTACGTTATTCTGAGAATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCa                                                                                                                       >  2:359730/1‑141 (MQ=255)
                  cAAGTCCGTTATAATGAGAATGATCGCTGCCTCTATTTCTTTAATATAACTGAACAAGTCCAAACAAATGAACTATATGAAAAGTCTAAACCAATCAGAGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAat                                                                                                                   <  2:279056/141‑1 (MQ=255)
                   aaGTACGTTATTCTGAGAATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAata                                                                                                                  >  1:407645/1‑141 (MQ=255)
                            aTTCTGAGAATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGACCTATACGACAATTCTAAACCAATCATTGCGACATTATTTTTACATAAGTACGATAAGATTACGCAAAATAGGAAGGATa                                                                                                         >  2:3849/1‑141 (MQ=255)
                            aTTCTGAGAATGATCACTGCCTCTATTTCTTTGATATAACGGAACAAGTACAGACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATa                                                                                                         >  2:651282/1‑141 (MQ=255)
                               cTGAGAATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGc                                                                                                      <  2:175175/141‑1 (MQ=255)
                                 gagaATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATAAGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGGAGCGATTACGCAAAAAATGAGTGAAACGCaa                                                                                                    >  1:153327/1‑140 (MQ=255)
                                     aTGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATATGCAACGtt                                                                                                >  1:436835/1‑141 (MQ=255)
                                         tCACTGCGTCTATTTCTTTGATATAACCGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACCATGAGAGGAGGCAAAATATGAATGATACGCAGCGTTCGGa                                                                                            <  1:122368/141‑1 (MQ=255)
                                         tCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGa                                                                                            <  1:726367/141‑1 (MQ=255)
                                         tCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGa                                                                                            >  2:335868/1‑141 (MQ=255)
                                         tCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGa                                                                                            >  1:24179/1‑141 (MQ=255)
                                          cACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGaa                                                                                           >  1:275334/1‑141 (MQ=255)
                                          cACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGaa                                                                                           >  2:210847/1‑141 (MQ=255)
                                              gCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCa                                                                                       <  2:409395/141‑1 (MQ=255)
                                                  ctATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGCTTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTc                                                                                   >  1:453703/1‑141 (MQ=255)
                                                  ctATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTc                                                                                   >  1:499799/1‑141 (MQ=255)
                                                    aTTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCaa                                                                                 >  2:386046/1‑141 (MQ=255)
                                                       tCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATgg                                                                              >  1:127937/1‑141 (MQ=255)
                                                        cTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGt                                                                             >  2:467546/1‑141 (MQ=255)
                                                            gATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAAcg                                                                         <  1:469563/141‑1 (MQ=255)
                                                             atatAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAAcgc                                                                        >  2:653382/1‑141 (MQ=255)
                                                               ctaACTGAACAAGGACAAACAAATGAACTATATGAAATTTCTAAACCAAGCATTGCGACCTTATTTTTGGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCgt                                                                      <  1:114695/140‑1 (MQ=255)
                                                                  aCTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAAACAATCATTGCTACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTc                                                                   >  2:611122/1‑141 (MQ=255)
                                                                         aaGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTc                                                            >  1:305099/1‑141 (MQ=255)
                                                                         aaGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGAAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTc                                                            >  1:415627/1‑141 (MQ=255)
                                                                            tACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGAt                                                         >  1:343934/1‑141 (MQ=255)
                                                                             acaaacaaATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCCAAATATGAATGATACGCAGCGTTCGGAAATAAACTCAATGGTAACGCGTGTCATTAGTCGATg                                                        >  1:337556/1‑141 (MQ=255)
                                                                             acaaacaaATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATg                                                        >  1:433317/1‑141 (MQ=255)
                                                                              caaacaaaTGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATgg                                                       <  1:306209/141‑1 (MQ=255)
                                                                              aaaaaaaaGGAACTATATGAAAATTCTAAACCAAGCATTGCGACATTAGTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGCAAATCAACTCAATGGTAACGCTTGTCATTAGTCGATgg                                                       <  1:483199/140‑1 (MQ=255)
                                                                               aaacaaaTGAATTATATGAACATTCTAAACCAATCATTGCGACATTATTTTTAGAGAAGTACGATGAGATTACGCAAAACATGAATGATACGCAGCGTTCGGAAATCAACTCAATGTTAACGCGTGTCATTAGTCGATggg                                                      <  1:554873/141‑1 (MQ=255)
                                                                                aacaaaTGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGc                                                     <  2:37164/141‑1 (MQ=255)
                                                                                 accaaTGAACTATATGAAAATTCGAAACCAATCATTGCGACCTTATTTTTAGATAAGTACCATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCa                                                    <  2:362477/138‑1 (MQ=255)
                                                                                   aaaTGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAAc                                                  >  1:358901/1‑141 (MQ=255)
                                                                                     aTGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATGAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCGACTg                                                >  1:501542/1‑141 (MQ=255)
                                                                                     aTGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGGTTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAGCTCAATGGTAACGCGTGTCATTAGTCGATGGGCGACTg                                                >  1:499449/1‑141 (MQ=255)
                                                                                     aTGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTg                                                <  2:499799/141‑1 (MQ=255)
                                                                                     aTGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTg                                                >  2:574071/1‑141 (MQ=255)
                                                                                     aTGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAACATGAGTGATACGAAGCGCTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAAcgg                                                >  1:114703/1‑139 (MQ=255)
                                                                                       gAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAg                                              <  1:82314/141‑1 (MQ=255)
                                                                                       gAACTATAGTAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGACAGGAACGCAGGGCTCGGAAATCATATACCTACAAACGCGTGTCATTAGTCGCTGCGCATTTGAg                                              <  2:610768/141‑1 (MQ=255)
                                                                                         aCTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCCAAATATGAATGATAGGCAGCGTGCGGAAATCAACTCAATGGTAACGCGGGTCATTAGTCGATGGGCAACTGAGta                                            >  1:162119/1‑141 (MQ=255)
                                                                                         aCTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAGTATGAATGATACGCGGGTAGCGATAATCAACTCAATGGTAACGAGTGTGATTAGTCGATGGGCAACTGAGta                                            >  1:162629/1‑141 (MQ=255)
                                                                                         aCTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGta                                            >  2:775805/1‑141 (MQ=255)
                                                                                          cTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGtat                                           <  2:445723/141‑1 (MQ=255)
                                                                                            atatGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGAGGGTCATTAATCGATGGGCAACTGAGTataa                                         >  1:398010/1‑141 (MQ=255)
                                                                                                gAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAatata                                     >  2:701914/1‑141 (MQ=255)
                                                                                                       cTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCa                              >  1:810241/1‑141 (MQ=255)
                                                                                                       cTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCa                              <  1:359730/141‑1 (MQ=255)
                                                                                                        tAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAAAGGTAAAGAGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCaa                             >  2:82676/1‑141 (MQ=255)
                                                                                                               aTCATTGCGACATTATTTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGGCGATGGGAAACGGAGCAGAATCTATTTTTCAACAGAGAc                      >  2:775170/1‑141 (MQ=255)
                                                                                                                tCATTGCGACATTATTTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACa                     >  1:186822/1‑141 (MQ=255)
                                                                                                                 taTTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAg                    <  2:519836/140‑1 (MQ=255)
                                                                                                                 cATTGCGGCATTTTTTTTAGATAAGTATGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAg                    <  2:88872/141‑1 (MQ=255)
                                                                                                                  aTTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGt                   >  1:740948/1‑141 (MQ=255)
                                                                                                                  aTTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGt                   <  2:821696/141‑1 (MQ=255)
                                                                                                                    tGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGGTc                 >  1:59379/1‑141 (MQ=255)
                                                                                                                     gCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTcc                <  1:386046/141‑1 (MQ=255)
                                                                                                                          attattTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAGTGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACGGTGCCGATCa           >  1:789627/1‑141 (MQ=255)
                                                                                                                               ttttAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCg      <  1:735782/141‑1 (MQ=255)
                                                                                                                                 ttAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAGCTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTa    >  1:216294/1‑141 (MQ=255)
                                                                                                                                 ttAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTa    >  1:456904/1‑141 (MQ=255)
                                                                                                                                   aGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTCCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGc  <  2:81113/141‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                      
TATAATCAGTATCATTTCCAAGTACGTTATTCTGAGAATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGC  >  CP000730/18753‑19024

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: