Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,233,504 T→A 83.1% V470D (GTT→GAT)  USA300HOU_1157 → non‑specific serine/threonine protein kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,233,5040TA83.1% 207.6 / 33.5 89V470D (GTT→GAT) USA300HOU_1157non‑specific serine/threonine protein kinase
Reads supporting (aligned to +/- strand):  ref base T (5/10);  new base A (32/41);  total (37/52)
Fisher's exact test for biased strand distribution p-value = 5.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

TGGTGAACGTGTTGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCTTTAGGCATTAAGCAAGTTTATGTAG  >  CP000730/1233367‑1233644
                                                                                                                                         |                                                                                                                                            
tGGTGAACGTGTTGAACGTGGCGACAGTGTTGATGTTGATATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATAc                                                                                                                                           <  1:886313/141‑1 (MQ=255)
  gTGAACGTGTTGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGa                                                                                                                                         <  1:696061/141‑1 (MQ=255)
   tGAACGTGTTGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGAt                                                                                                                                        <  1:1034410/141‑1 (MQ=255)
   tGAACGTGTTGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGAt                                                                                                                                        <  2:876609/141‑1 (MQ=255)
    gAACGTGTTGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGAtt                                                                                                                                       <  1:99768/141‑1 (MQ=255)
      aCGTGTTGAACGTGGTGATAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGa                                                                                                                                     <  2:63086/141‑1 (MQ=255)
      aCGTGTTGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACgcaaga                                                                                                                                     >  2:895716/1‑136 (MQ=255)
      aCGTGTTGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGa                                                                                                                                     >  2:1040469/1‑141 (MQ=255)
         tgtTGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGaaaa                                                                                                                                  >  1:157385/1‑141 (MQ=255)
                cGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGtatat                                                                                                                           <  2:245601/141‑1 (MQ=255)
                      gACAGTGTTGATGTTGATATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATAATAAAAAAAAACAAGTATATaataat                                                                                                                     <  2:719749/141‑1 (MQ=255)
                       aCAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATc                                                                                                                    >  1:404212/1‑141 (MQ=255)
                          gtgtTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAg                                                                                                                 <  2:590856/141‑1 (MQ=255)
                           tgtTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGACGATACGATTGAAAAAGTATATAATAATCAAgc                                                                                                                <  2:338395/141‑1 (MQ=255)
                             ttGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAgcgc                                                                                                              >  1:202808/1‑141 (MQ=255)
                              tGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGCTACGATTGAAAAAGTATATAATAATCAAGCGcc                                                                                                             <  2:1050950/141‑1 (MQ=255)
                               gATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCa                                                                                                            >  2:653968/1‑141 (MQ=255)
                                  gttgttATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAAAAATCAAGCGCCaaaa                                                                                                         >  1:971539/1‑141 (MQ=255)
                                   ttgttATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATATTCATGCGCCaaatg                                                                                                        >  2:932282/1‑139 (MQ=255)
                                   ttgttATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCaaaaa                                                                                                        <  2:672480/141‑2 (MQ=255)
                                   gtgttATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAg                                                                                                        <  2:226340/140‑1 (MQ=255)
                                      ttATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGAt                                                                                                     >  1:125505/1‑141 (MQ=255)
                                      ttATATCAAAGGGCCATGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCTAAAGGAt                                                                                                     <  1:662375/141‑1 (MQ=255)
                                        atatCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTCTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTGAAATATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATAc                                                                                                   <  1:508645/141‑1 (MQ=255)
                                           tCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACAtt                                                                                                >  1:232471/1‑141 (MQ=255)
                                            cAAAGGGCCCTGAAAAGGTTAACATGCCAAATGTCTTTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGTAAAAGTATATAATAATCAACAGCGAAAAGGATACAtta                                                                                               >  2:164762/1‑140 (MQ=255)
                                            cAAAGGGCCCTGAAAAGGTTAAAATGCCACGTCTCACATGCGTACCTCCAGAGGGCGGAAATCGGAAATTAAAAGCGGTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTg                                                                                               >  2:1019098/1‑141 (MQ=255)
                                            cAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTACAGATGAGACGATTGAAATAGTATACAATAATCAAGCGCCAAAAGGGTACATTg                                                                                               <  1:672032/141‑1 (MQ=255)
                                            cAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTg                                                                                               <  1:790764/141‑1 (MQ=255)
                                            cAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACAAAAAAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTg                                                                                               <  2:132520/141‑1 (MQ=255)
                                            cAAAGGGCCCTGAAAAGGTTAAAATCCCCCCTGCCCTTTGGTTAACTAAGGCGTAATCCCTGCAGAAACTAAAATCGCCAGGTCTTAAAGATGACACGATTGAAAAAGTATATACTAATCAAGCGCCAAAAGGACCCAttc                                                                                               >  2:1019218/1‑140 (MQ=255)
                                            cAAAGGGCCCTGAAAAGGATAAAAAGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGTTCTTAAAGATAAAAAAAATAAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTg                                                                                               <  2:458990/141‑1 (MQ=255)
                                            cAAAGGGACCTGAAAACGATAAAATGACAAATGACATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGAAAAAAAAAAAAAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTg                                                                                               <  2:1035583/141‑1 (MQ=255)
                                             aaaGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGc                                                                                              >  1:1030544/1‑141 (MQ=255)
                                             aaaGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGc                                                                                              >  1:1030581/1‑141 (MQ=255)
                                              aaGGGCCTTGAAAAGGTTAAAATGCCAAATGTCATTTGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCa                                                                                             <  1:74935/141‑1 (MQ=255)
                                                gggCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGCTAGGTCTTAAAGATGAAACGATCGAAAAAGCATATAATAATCAAGCGCCAAAAGAATACATAGCaaa                                                                                           >  2:28125/1‑141 (MQ=255)
                                                    ccTGAAAAGGTTAAAAGGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCaa                                                                                       >  1:30868/1‑141 (MQ=255)
                                                     cTGAAAAGGTTAAACTGCCAAATGTCATGGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTgcaaagcaaa                                                                                      >  1:691394/1‑141 (MQ=255)
                                                     cTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCaaa                                                                                      >  2:988408/1‑141 (MQ=255)
                                                     cTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCaaa                                                                                      >  2:545095/1‑141 (MQ=255)
                                                     cTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAAACaaa                                                                                      >  2:20656/1‑141 (MQ=255)
                                                      tGAAAAGGTTAAAATTCTAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGCTACATTGCAAATCaaaa                                                                                     >  2:757277/1‑140 (MQ=255)
                                                      tGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAg                                                                                     <  2:280252/141‑1 (MQ=255)
                                                      tGAAAAGGTTAAAAGGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAg                                                                                     <  2:88689/141‑1 (MQ=255)
                                                       gAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAgt                                                                                    >  1:985998/1‑141 (MQ=255)
                                                       gAAAAGGTTAAAATGCCAAATGTCAGTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAgt                                                                                    <  2:603084/141‑1 (MQ=255)
                                                            ggTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAAc                                                                               <  1:541186/141‑1 (MQ=255)
                                                            ggTCAAAATGCCAAATGACATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATAAAAAAAATAAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAAc                                                                               <  2:746451/141‑1 (MQ=255)
                                                              ttAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCTTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCg                                                                             <  2:540931/141‑1 (MQ=255)
                                                               tAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGc                                                                            <  1:517286/141‑1 (MQ=255)
                                                                      ccAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCGGAAATGAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACGGCAAAAACt                                                                     >  2:529710/1‑141 (MQ=255)
                                                                      ccAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTACAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACt                                                                     >  2:852791/1‑141 (MQ=255)
                                                                      ccAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACt                                                                     >  1:193700/1‑141 (MQ=255)
                                                                      ccAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACt                                                                     <  1:626373/141‑1 (MQ=255)
                                                                      ccAAATGTCATTGGATTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACt                                                                     >  1:556079/1‑141 (MQ=255)
                                                                        aaaaGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAATTACTGa                                                                   <  2:666584/141‑1 (MQ=255)
                                                                           tgataTTGGTTTACCTAAGGAGCAAGCCGTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAAGCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAAt                                                                <  1:603530/137‑1 (MQ=255)
                                                                           tGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAAt                                                                <  1:271068/141‑1 (MQ=255)
                                                                                 tGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTAt                                                          <  1:953178/141‑1 (MQ=255)
                                                                                  ggTTTACCTAAGGAGGAAGCCTTGCAGAAAATAAAATCGTTAGGTCTTAAAGAAAAAAAAAAAAAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTAtt                                                         <  2:1033349/141‑1 (MQ=255)
                                                                                  ggTTGACCTAATGATGAAGCCTTGCAGAAATTAAAATCGTTAGGTATCAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTAtt                                                         <  1:54659/141‑1 (MQ=255)
                                                                                   gTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTc                                                        <  1:318642/141‑1 (MQ=255)
                                                                                    tttACCTAAGGAGGAAGCCTGGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCa                                                       <  2:533904/141‑1 (MQ=255)
                                                                                         cTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGCTTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTACTGAtt                                                  >  2:692025/1‑141 (MQ=255)
                                                                                           aaGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCt                                                <  1:22853/141‑1 (MQ=255)
                                                                                            aggaggAAGCCTTGCAGAAGTTAAAATCGTTAGGTCATAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTa                                               >  1:561967/1‑141 (MQ=255)
                                                                                             ggaggaAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCTAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTaa                                              <  2:404212/141‑1 (MQ=255)
                                                                                              gaggaAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAat                                             >  1:651275/1‑141 (MQ=255)
                                                                                              gaggaAGCCTTGCAGAAATTAAAATCGTTAAGTATTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAATGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAat                                             <  2:669992/141‑1 (MQ=255)
                                                                                               aggaAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCACCCGCCTTAATGTAAGAGAGGAATGCGAAGTGAACTAATGAATACTGAAATCGCTATTCATGTTTCTAata                                            <  2:748122/141‑1 (MQ=255)
                                                                                               aggaAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAata                                            <  2:913644/141‑1 (MQ=255)
                                                                                                ggaAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAatat                                           <  2:2898/141‑1 (MQ=255)
                                                                                                    gCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTaaa                                       >  1:101338/1‑141 (MQ=255)
                                                                                                    gCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTaaa                                       >  1:658161/1‑141 (MQ=255)
                                                                                                          cAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTatat                                 <  1:413368/141‑1 (MQ=255)
                                                                                                          cAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTatat                                 <  1:93868/141‑1 (MQ=255)
                                                                                                                  aggaTCGTTAGGTCTTGAAGATGATACGATTGAAAGAGTATATAGTAGTCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCttt                         <  2:290587/138‑1 (MQ=255)
                                                                                                                  aaaaTCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCttt                         <  1:862334/141‑1 (MQ=255)
                                                                                                                  aaaaTCGTAAGGCTTTAAAGCTGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCttt                         <  2:971539/141‑1 (MQ=255)
                                                                                                                    aaTCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCTTTAg                       >  1:204037/1‑141 (MQ=255)
                                                                                                                    aaTCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCTTTAg                       <  2:1035129/141‑1 (MQ=255)
                                                                                                                     aTCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAGTATTAAACTATATGAATCTTTAgg                      >  2:663308/1‑141 (MQ=255)
                                                                                                                      tCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAGTCTTTAGGc                     >  2:369111/1‑141 (MQ=255)
                                                                                                                        gTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATAGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCTTGAGGCAt                   >  1:96795/1‑141 (MQ=255)
                                                                                                                          tAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCTTTAGGCATTa                 >  2:1012778/1‑141 (MQ=255)
                                                                                                                          tAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCTTTAGGCATTa                 <  1:658546/141‑1 (MQ=255)
                                                                                                                                  aaaGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCTTTAGGCATTAAGCAAGtt         <  1:624060/141‑1 (MQ=255)
                                                                                                                                  aaaGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCTTTAGGCATTAAGCAAGtt         <  1:797896/141‑1 (MQ=255)
                                                                                                                                     gatgatACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCTTTAGGCATTAAGCAAGTTTAt      >  2:519623/1‑141 (MQ=255)
                                                                                                                                     gatgatACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACAGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCTTTAGGCATTAAGCAAGTTTAt      >  1:715401/1‑141 (MQ=255)
                                                                                                                                         ttACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATATTTAGGCATTAAGCAAGTTTATGTag  >  2:121427/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                                            
TGGTGAACGTGTTGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCTTTAGGCATTAAGCAAGTTTATGTAG  >  CP000730/1233367‑1233644

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: