Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 2,177,618 A→C A126A (GCT→GCG U712_10895 ← Dihydrodipicolinate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,177,6180AC97.1% 220.8 / ‑4.3 70A126A (GCT→GCGU712_10895Dihydrodipicolinate reductase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (35/33);  total (36/34)
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

CTGTTTTCAGCGCCGTTCCGCTTGGTGCGTCAAGCTTCCGATCATGGTGTAGCTCAATAATCTCAACATCCTCAAAGTAGTTGGCAGCCATTTTCGAAAATTTCATCATCAGCACCGCACCAAGTGCAAAGTTTGGAGCGATGATGGCTCCGATTCCCTTCTCTTCTGTTAATGACGTAAGCTCTTTTAAATCAGCTTCTGAGAAACCGGTTGTTCCAACAACCGGACGAACGCCGTGCTCTAATGCAATTTTTGTATGTACTTTTCCG  >  CP006881/2177482‑2177750
                                                                                                                                        |                                                                                                                                    
cTGTTTTAAGCGCAGTTCCGCTTGGTGCGTCAAGCTTCTGGTCATGATGAAGCTCAATAATCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGgcgc                                                                                                                                    <  2:96216/139‑1 (MQ=255)
      taaGCGCAGTTCCGCTTGGTGCGTCAAGCTTCTGGTCATGATGAAGCTCAATAATCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCgatga                                                                                                                               <  8:438626/136‑1 (MQ=255)
             agTTCCGCTTGGTGCGTCAAGCTTCTGGTCATGATGAAGCTCAATAATCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTcc                                                                                                                        >  6:9732/2‑138 (MQ=255)
             agTTCCGCTTGGTGCGTCAAGCTTCTGGTCATGATGAAGCTCAATAATCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTcc                                                                                                                        >  2:149035/2‑138 (MQ=255)
             agTTCCGCTTGGTGCGTCAAGCTTCTGGTCATGATGAAGCTCAATAATCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCg                                                                                                                       <  3:358883/138‑1 (MQ=255)
             agTTCCGCTTGGTGCGTCAAGCTTCTGGTCATGATGAAGCTCAATAATCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCg                                                                                                                       <  3:82903/138‑1 (MQ=255)
               ttCCGCTTGGTGCGTCAAGCTTCTGGTCATGATGAAGCTCAATAATCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCTCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGAt                                                                                                                     >  5:347756/1‑139 (MQ=255)
                 ccGCTTGGTGCGTCAAGCTTCTGGTCATGATGAAGCTCAATAATCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGAtc                                                                                                                    <  3:69440/138‑2 (MQ=255)
                    cTTGGTGCGTCAAGCTTCTGGTCATGATGAAGCTCAATAATCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGCGCGCAAAATTTGGCGCGATGATGGCTCCGATCCctt                                                                                                               >  6:245617/1‑139 (MQ=255)
                         tGCGTCAAGCTTCTGGTCATGATGAAGCTCAATAATCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCtttttct                                                                                                           <  3:166552/139‑1 (MQ=255)
                         tGCGTCAAGCTTCTGGTCATGATGAAGCTCAATAATCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCtttttct                                                                                                           <  1:42407/139‑1 (MQ=255)
                         tGCGTCAAGCTTCTGGTCATGATGAAGCTCAATAATCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCtttttct                                                                                                           <  5:321494/139‑1 (MQ=255)
                           cGTCAAGCTTCTGGTCATGATGAAGCTCAATAATCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCtttttcttc                                                                                                         <  1:8004/139‑1 (MQ=255)
                           cGTCAAGCTTCTGGTCATGATGAAGCTCAATAATCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCtttttcttc                                                                                                         >  2:8004/1‑139 (MQ=255)
                           cGTCAAGCTTCTGGTCATGATGAAGCTCAATAATCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCtttttctt                                                                                                          <  3:267530/138‑1 (MQ=255)
                             tCAAGCTTCTGGTCATGATGAAGCTCAATAATCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGAtc                                                                                                                    <  7:6301/126‑2 (MQ=255)
                             tCAAGCTTCTGGTCATGATGAAGCTCAATAATCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGAtc                                                                                                                    >  8:6301/1‑125 (MQ=255)
                             tCAAGCTTCTGGTCATGATGAAGCTCAATAATCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTAATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATTGCTCCGAGACCTTTTTCTTCTg                                                                                                       >  4:267530/1‑139 (MQ=255)
                                     ctggtCATGATGAAGCTCAATAATCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTGTGTTAAAGat                                                                                               >  2:28701/5‑138 (MQ=255)
                                      tggtCATGATGAAGCTCAATAATCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGatg                                                                                              >  2:372904/4‑137 (MQ=255)
                                        gtCATGATGAAGCTCAATAATCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGat                                                                                               <  5:157955/135‑2 (MQ=255)
                                           atgatgAAGCTCACTAATCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGc                                                                                         <  6:13823/139‑1 (MQ=255)
                                                  aGCTCAATAATCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCgatga                                                                                                                               >  5:43617/1‑94 (MQ=255)
                                                  aGCTCAATAATCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCgatga                                                                                                                               <  6:43617/94‑1 (MQ=255)
                                                       aataatCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTa                                                                                                    >  4:60764/1‑116 (MQ=255)
                                                       aataatCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTa                                                                                                    >  4:54796/1‑116 (MQ=255)
                                                       aataatCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTa                                                                                                    <  3:54796/116‑1 (MQ=255)
                                                       aataatCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTa                                                                                                    <  3:60764/116‑1 (MQ=255)
                                                       aataatCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAAt                                                                               >  1:66193/1‑137 (MQ=255)
                                                       aataatCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAAt                                                                               <  2:40509/137‑1 (MQ=255)
                                                         taatCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCtttttct                                                                                                           <  4:344223/107‑1 (MQ=255)
                                                         taatCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCtttttct                                                                                                           >  3:344223/1‑107 (MQ=255)
                                                         taatCTCAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGc                                                                           <  3:286092/139‑1 (MQ=255)
                                                              tcAACATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCt                                                                       >  3:351243/1‑138 (MQ=255)
                                                                aaCATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCt                                                                       <  1:372904/136‑1 (MQ=255)
                                                                 aCATCCTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCTgag                                                                    >  1:340462/1‑138 (MQ=255)
                                                                     ccTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCTGAGAAAcc                                                               >  5:196194/1‑139 (MQ=255)
                                                                     ccTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATAGCACCGAGCTCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCTGAGAAAcc                                                               >  1:252267/1‑139 (MQ=255)
                                                                      cTCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCCCCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAATATGTGAGCTCTTTTATATCAGCTTCTGAGAAAcag                                                              >  7:277939/1‑137 (MQ=255)
                                                                       tCAAAATAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCTGAGATACCgg                                                             >  8:429855/1‑139 (MQ=255)
                                                                         aaaaTAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCTGAGAAACCGgtt                                                           <  2:80995/139‑1 (MQ=255)
                                                                         aaaaTAGTTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCTGAGAAACCGgtt                                                           <  6:171362/139‑1 (MQ=255)
                                                                            atagtTGGCAGCAATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCTGAGAAACCGgttgtt                                                        <  6:347756/138‑1 (MQ=255)
                                                                               gtTGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCTGAGAAACCGGTTGTTcc                                                      <  7:72192/138‑1 (MQ=255)
                                                                                ttGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCTGAGAAACCGGTTGTTcc                                                      <  7:144751/137‑1 (MQ=255)
                                                                                ttGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCTGAGAAACCGGTTGTTcc                                                      >  8:144751/1‑137 (MQ=255)
                                                                                ttGGCAGCCATTTTCGAAAATTTCATCATCAGCACCGCACCAAGTGCAAAGTTTGGAGCGATGATGGCTCCGATTCCctt                                                                                                               >  5:153727/1‑80 (MQ=255)
                                                                                ttGGCAGCCATTTTCGAAAATTTCATCATCAGCACCGCACCAAGTGCAAAGTTTGGAGCGATGATGGCTCCGATTCCctt                                                                                                               <  6:153727/80‑1 (MQ=255)
                                                                                 tGGCAGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTATGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCTGAGAAACCGgttgttacta                                                    >  8:143352/1‑134 (MQ=255)
                                                                                     aGCCATTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCTTAGAAACCGGTTGTTCCgacaact                                               >  5:312440/1‑138 (MQ=255)
                                                                                         aTTTTTGAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCTGAGAAACCGGTTGTTCCGACAACTGGAc                                           <  2:290875/139‑1 (MQ=255)
                                                                                             ttgAAAATTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCTGAGAAACCGGTTGTTCCGACAACTGGACgga                                        <  5:9732/136‑3 (MQ=255)
                                                                                                 aaaTTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCTGAGAAACCGGTTGTTCCGACAACTGGACGGACTCCg                                   <  2:340462/139‑1 (MQ=255)
                                                                                                 aaaTTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCTGAGAAACCGGTTGTTCCGACAACTGGACGGACTCCg                                   <  2:89560/139‑1 (MQ=255)
                                                                                                 aaaTTTCATCATCAGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCTGAGAAACCGGTTGTTCCGACAACTGGACGGACTCCg                                   <  8:349450/139‑1 (MQ=255)
                                                                                                         tcatcaGTATCGCACCGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCTGAGAAACCGGTTGTTCCGACAACTGGACGGACTCCGTGCTCTaa                           <  6:196194/139‑1 (MQ=255)
                                                                                                                     caccGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCTGAGAAACCGGTTGTTCCGACAACTGGACGGACTCCGTGCTCTAATGCAATTTTtgt               >  2:162856/1‑139 (MQ=255)
                                                                                                                     caccGAGCGCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATAAGCTTCTGAGAAACCGGTTGTTCCGACAACTGGACGGACTCCGTGCTCTAATGCAATTTTtg                >  3:335915/1‑138 (MQ=255)
                                                                                                                        cgagcgCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCTGAGAAACCGGTTGTTCCGACAACTGGACGGACTCCGTGCTCTAATGCAATTTTtgtatg            <  1:162856/134‑1 (MQ=255)
                                                                                                                         gagcgCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCTGAGAAACCGGTTGTTCCGACAACTGGACGGACTCCGTGCTCTAATGCAATTTTtgtatg            >  5:195886/5‑138 (MQ=255)
                                                                                                                          cgcgCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCTGAGATACCGGTTGTTCGGCCATCTGGACGGCCTCCGTGCTCTAATGCAATTTTtgtatgta          >  7:259039/4‑139 (MQ=255)
                                                                                                                          agcggaaaaTTTGGCGCGATGATGGCTCCGATCCCTGTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCGTCTGAGAAACCGGTTGTTCCGACACCTGGACGGAGTCCGTGCTCTAATGCAATTTTtgtatg            >  2:71260/6‑137 (MQ=255)
                                                                                                                          agcgCAAAATTTGGCGCGATGATGGCTCCGCTCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCTGAGAAACCGGTTGTTCCGACAACTGGACGGACTCCGTGCTCTAATGCAATTTTtgtatgta          >  8:375171/4‑139 (MQ=255)
                                                                                                                          agcgCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCTGAGAAACCGGTTGTTCCGACAACTGGAc                                           >  2:235218/4‑106 (MQ=255)
                                                                                                                          agcgCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCTGAGAAACCGGTTGTTCCGACAACTGGAc                                           <  1:235218/103‑1 (MQ=255)
                                                                                                                          agcgCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCTGAGAAACCGGTTGTTCCGACAACTGGACGGACTCCGTGCTCTAATGCAATTTTtgtatgta          >  1:192981/4‑139 (MQ=255)
                                                                                                                          agcgCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCTGAGAAACCGGTTGTTCCGACAACTGGACGGACTCCGTGCTCTAATGCAATTTTtgtatgta          >  1:216012/4‑139 (MQ=255)
                                                                                                                          agcgCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCTGAGAAACCGGTTGTTCCGACAACTGGACGGACTCCGTGCTCTAATGCAATTTTtgtatgta          >  3:126348/4‑139 (MQ=255)
                                                                                                                          agcgCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCTGAGAAACCGGTTGTTCCGACAACTGGACGGACTCCGTGCTCTAATGCAATTTTtgtatg            >  3:125217/4‑137 (MQ=255)
                                                                                                                          agcgCAAAATTTGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCTGAGAAACCGGTGGTGCCGACAACTGGACGGACTCCGTGCTCTACTGCAATTTTtgtatg            >  3:245601/4‑137 (MQ=255)
                                                                                                                                  atttGGCGCGATGATGGCTCCGATCCCTTTTTCTTCTGTTAAAGATGTGAGCTCTTTTAAATCAGCTTCTGAGAAACCGGTTGTTCCGACAACTGGACGGACTCCGTGCTCTAATGCAATTTTTGTATGTACTTTTCCg  <  3:290010/138‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                    
CTGTTTTCAGCGCCGTTCCGCTTGGTGCGTCAAGCTTCCGATCATGGTGTAGCTCAATAATCTCAACATCCTCAAAGTAGTTGGCAGCCATTTTCGAAAATTTCATCATCAGCACCGCACCAAGTGCAAAGTTTGGAGCGATGATGGCTCCGATTCCCTTCTCTTCTGTTAATGACGTAAGCTCTTTTAAATCAGCTTCTGAGAAACCGGTTGTTCCAACAACCGGACGAACGCCGTGCTCTAATGCAATTTTTGTATGTACTTTTCCG  >  CP006881/2177482‑2177750

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