Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 2,198,953 C→T K10K (AAG→AAA U712_11015 ← Chemotaxis protein methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,198,9530CT98.9% 295.1 / NA 94K10K (AAG→AAAU712_11015Chemotaxis protein methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (59/34);  minor base A (1/0);  total (60/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

CTTGATCCTTTTCTAATGCCGCCGCAAAGTCCTTAAAGCTTTGAAACCCCTTTTTCTCATAAAGCGACGTCAGCCTTCTCTTCATTTGTGCTTCTTTATATAATGTCAAATCGACTCCGGTTAATTGTTTCCACTTCGTTGTAAATACGCTGTATGTATCCATCTCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAAATTGTATTGTTTCAATAAAAAAATCAGCTTTGCAGCT  >  CP006881/2198820‑2199078
                                                                                                                                     |                                                                                                                             
cTTGATCCTTTTCCAATGCCGCAGCAAAGTCCTTAAAGCTTTGGAACCCCTTTTTCTCATAAAGTGACGTTAGCCTTCTCTTCATTTGTGGTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGt                                                                                                                          >  3:40109/1‑139 (MQ=255)
cTTGATCCTTTTCCAATGCCGCAGCAAAGTCCTTAAAGCTTTGGAACCCCTTTTTCTCATAAAGTGACGTTAGCCTTCTCTTCATTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGt                                                                                                                          >  4:13094/1‑139 (MQ=255)
cTTGATCCTTTTCCAATGCCGCAGCAAAGTCCTTAAAGCTTTGGAACCCCTTTTTCTCATAAAGTGACGTTAGCCTTCTCTTCATTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGt                                                                                                                          >  1:99688/1‑139 (MQ=255)
cTTGATCCTTTTCCAATGCCGCAGCAAAGTCCTTAAAGCTTTGGAACCCCTTTTTCTCATAAAGTGACGTTAGCCTTCTCTTCATTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGt                                                                                                                          >  7:54449/1‑139 (MQ=255)
cTTGATCCTTTTCCAATGCCGCAGCAAAGTCCTTAAAGCTTTGGAACCCCTTTTTCTCATAAAGTGACGTTAGCCTTCTCTTCATTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGt                                                                                                                          >  8:356225/1‑139 (MQ=255)
cTTGATCCTTTTCCAATGCCGCAGCAAAGTCCTTAAAGCTTTGGAACCCCTTTTTCTCATAAAGTGACGTTAGCCTTCTCTTCATTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATGTCGt                                                                                                                          >  5:89655/1‑139 (MQ=255)
cTTGATCCTTTTCCAATGCCGCAGCAAAGTCATTAAAGCTTTGGAACCCCTTTTTCTCATAAAGTCACGTTAGCCTTCTCTTCATTTGTGCTTCTTTATATAGGGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGt                                                                                                                          <  1:184476/139‑1 (MQ=255)
   gATCCTTTTCCAATGCCGCAGCAAAGTCCTTAAAGCTTTGGAACCCCTTTTTCTCATAAAGTGACGTTAGCCTTCTCTTCATTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTg                                                                                                                        >  3:386215/1‑138 (MQ=255)
   gATCCTTTACCAATGCCGCAGCAAAGTCCTTAAAGCTTTGGAACCCATTTTTCTCATAAAGTGACGTTAGCCTTCTCTTCATTTGTGCTTCTTTAAAAAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTg                                                                                                                        <  1:130246/138‑1 (MQ=255)
        ttttCCAATGCCGCAGCAAAGTCCTTAAAGCATTCGAACCACTTTTTCTCATAAAGAGACGTTACCCTTATCTTCATTAGTGCTTCTTTTATTAGTGTTAAATCGACTCCGGTTAAAACTTTCCAATTCGTAGTaaaa                                                                                                                   >  7:142967/1‑137 (MQ=255)
         tttCCAATGCCGCAGCAAAGTCCTTAAAGCTTTGGAACCCCTATTTCTCATAAAGTGACGTTAGCCTACACTTCATTTGTGCTTCTTTATATAGTGATAAAAAGACTCCGGTTAATAGTTTCCATTTCGATGTAAATAc                                                                                                                 >  5:136668/1‑139 (MQ=255)
           tccaaTGCCGCAGCAATGTCATTAACGCTTTGTAACCCCTTTTTGTCATAATGTGAGGTTCGCCTTCTCTTCATTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACg                                                                                                                <  8:54449/135‑1 (MQ=255)
           tccaaTGCCGCAGCAAAGTCCTTAAAGCTTTGGAACCCCTTTTTCTCATAAAGTGACTTTAGCCTTCTCTTCATTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACg                                                                                                                <  3:114875/135‑1 (MQ=255)
           tccaaTGCCGCAGCAAAGTCCTTAAAGCTTTGGAACCCCTTTTTCTCATAAAGTGACGTTAGCCTTCTCTTCATTTGTGCTTCTTTATATCGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACg                                                                                                                <  5:188295/135‑1 (MQ=255)
           tccaaTGCCGCAGCAAAGTCCTTAAAGCTTTGGAACCCCTTTTTCTCATAAAGTGACGTTAGCCTTCTCTTCATTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACg                                                                                                                <  4:435642/135‑1 (MQ=255)
            gaattGCGGCAGCAATGTCGTTAATGCTTTGGATCCCCTTTTTTTCATAATGTGACGTTTGCGTTCTCTTCATTTGTGCTTCTTTATATTGTGTTAATTCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGc                                                                                                               <  6:136668/134‑1 (MQ=255)
            cctaTGCCGCAGCAATGTCCTTAAAGCTTTTGAACCGCTTTTCCAGATAATTTTAAGTTAGCGTTGTCTTTATTTGTGGTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATTCGc                                                                                                               <  4:325497/135‑1 (MQ=255)
                tGCCGCAGCAACGTCCTTAAAGCTTTGGAACCCCTTTTTCTCATAAAGTGACGTTAGCCTTCTCTTCATTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCtgtat                                                                                                          <  6:412512/139‑1 (MQ=255)
                 gccgaagcAAAGTCCTTAAAGCTTTGGAACCCCTTTTTCTCATAAAGTGACGTTAGCCTTCTCTTCATTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCtgtat                                                                                                          <  8:216992/132‑1 (MQ=255)
                  ccgcagCAAAGTCCTTAAAGCTTTGGAACCCCTTTTTCTCATAAAGTGACGTTAGCCTTCTCTTCATTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAATTACGCtgtatg                                                                                                         >  8:413744/1‑138 (MQ=255)
                  ccgcaccAAAGTCCTTAAAGTTTTGGAACCCCTTTTTCTAATAAAGTGAGGTTAGCCTTCTCTTCTTTTGTGATTCTTTATATAGTTTAAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCtgtatg                                                                                                         <  3:39347/132‑1 (MQ=255)
                    gcagcAAAGTCCTTAAATCTTTGTAAGCCGTTTATCTCATAAAGTGACGTAAGCCGTCTATTCATTGGTGCTTCTTTATATAGTTTTAACTCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTgtatgtat                                                                                                      <  8:321512/136‑1 (MQ=255)
                    gcagcAAAGTCCTTAAAGCTTTGGAACCCCTTTTTCTCATAAAGTGACGTTAGCCTTCTCTTCATTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAATTACGCTgtatgtat                                                                                                      >  4:79918/4‑139 (MQ=255)
                     cagaaaaGTCCTTAAAGCTTTGGAACCCCTTTTTCTCATAAAGTGACGTTAGCCTTCTCTTAATTTGTGCTTATTTATATAGTTTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATc                                                                                                     <  6:253794/135‑1 (MQ=255)
                        cAATGTCTTAAAAGTTTTGGAACCCTTTTTTCTAATAATGTTCGTTTAGCGTTCTATTCTTTTTTGCTTATTTATTTTGTGTAAACTCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCAt                                                                                                  <  4:40109/139‑1 (MQ=255)
                        aaatttCCTTAAAGGTTTGGAACCCCTTTTTCTCATAACGTGACGTTAGCCTTCTCTTCATTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCAt                                                                                                  <  4:198385/134‑1 (MQ=255)
                         aaaGTCCTTAAAGCTTTGGAACCCCTTTTTCTCATAAAGTGACGTTAGCCTTCTCTTCATTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCAtc                                                                                                 <  5:208442/139‑1 (MQ=255)
                            gTCCTTAAAGCTTTGGAACCCCTTTTTCTCATAAAGTGACGTTAGCCTTCTCTTCATTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCAtctc                                                                                               >  3:419173/1‑138 (MQ=255)
                             tCCTTAAAGCTTTGGAACCCCTTTTTCTCATAAAGTGACGTTAGCCTTCTCTTCATTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCAtctcta                                                                                             <  3:13094/139‑2 (MQ=255)
                                       tttGGAACCCCTTTTTCTCATAAAGTGACGTTAGCCTTCTCTTCATTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGt                                                                                 >  5:98219/1‑139 (MQ=255)
                                         tggaaCCCCTTTTTCTCATAAAGTGACGTTAGCCTTCTCTTCATTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGtct                                                                               >  4:84557/4‑139 (MQ=255)
                                                 cTTTTTCTCATAAAGTGACGTTAGCCTTCTCTTCATTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCt                                                                       >  5:348675/1‑139 (MQ=255)
                                                    tttCTCATAAAGTGACGTTAGCCTTCTATTCTTTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCttt                                                                     <  3:163253/138‑1 (MQ=255)
                                                    tttCTAATAAAGTGACGTAAGCATTATCTTCATTTTTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCttt                                                                     <  2:99688/138‑1 (MQ=255)
                                                         cATAAAGTGACGTTAGCCTTCTCTTCATTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTcg                                                               >  1:335978/1‑138 (MQ=255)
                                                          aTAAAGTGACGTTAGCCTTCTCTTCATTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCtgtatgta                                                                                                       >  7:82927/1‑100 (MQ=255)
                                                          aTAAAGTGACGTTAGCCTTCTCTTCATTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCtgtatgta                                                                                                       <  8:82927/100‑1 (MQ=255)
                                                          aTAAAGTGACGTTAGCCTTCTCTTCATTTGTGCTTCTTTATATAGGGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTc                                                                <  7:92376/137‑1 (MQ=255)
                                                                  aCGTTAGCCTTCTCTTCATTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACg                                                                                                                <  7:59642/83‑1 (MQ=255)
                                                                  aCGTTAGCCTTCTCTTCATTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACg                                                                                                                >  8:59642/1‑83 (MQ=255)
                                                                  aCGTTAGCCTTCTCTTCATTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCt                                                                           >  8:354736/1‑118 (MQ=255)
                                                                  aCGTTAGCCTTCTCTTCATTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCt                                                                           <  7:354736/118‑1 (MQ=255)
                                                                  aCGTTAGCCTTCTCTTCATTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGc                                                      >  4:195976/1‑139 (MQ=255)
                                                                  aCGTTAGCCTTCTCTTCATTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGc                                                      >  6:87011/1‑139 (MQ=255)
                                                                       aGCCTTCTCTTCATTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGtt                                                 >  2:347332/1‑139 (MQ=255)
                                                                             ctctTCATTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTAt                                             >  8:397323/1‑137 (MQ=255)
                                                                                   aTTTGTGCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGaaa                                      <  6:348675/138‑1 (MQ=255)
                                                                                        tgCTTCTTTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAAattgtat                                >  3:49274/1‑139 (MQ=255)
                                                                                          cttcttTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAAattgtattg                              >  8:236898/1‑139 (MQ=255)
                                                                                          cttcttTATATAGTGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAAattgtattg                              >  6:208587/1‑139 (MQ=255)
                                                                                                      gtGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAAATTGTATTGTTTCAATaaaa                   >  5:259418/2‑138 (MQ=255)
                                                                                                       tGTTAACTCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATc                                                                                                     <  2:390839/57‑1 (MQ=255)
                                                                                                       tGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATc                                                                                                     >  1:390839/1‑57 (MQ=255)
                                                                                                       tGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTg                                                                  >  5:26046/1‑90 (MQ=255)
                                                                                                       tGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTg                                                                  <  6:26046/90‑1 (MQ=255)
                                                                                                       tGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTTGTATAGCTTGTTTTATGAGGAAATTGTATTGTTTCAATaaaaa                  >  5:17894/1‑138 (MQ=255)
                                                                                                       tGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGt                                                           <  4:422658/97‑1 (MQ=255)
                                                                                                       tGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGt                                                           >  3:422658/1‑97 (MQ=255)
                                                                                                       tGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATATCTTGTTTTATCAGGAAATTGTTTTGTTTCAATaaaaa                  >  5:424497/1‑138 (MQ=255)
                                                                                                       tGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGGTTGTTTTATCTGGATATTGTATTGTTTCAATaaaaa                  >  6:59876/1‑138 (MQ=255)
                                                                                                       tGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAgg                                         >  7:220728/1‑115 (MQ=255)
                                                                                                       tGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAgg                                         <  8:220728/115‑1 (MQ=255)
                                                                                                       tGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAATTTGTATTGTTTCAATaaaaa                  >  3:80069/1‑138 (MQ=255)
                                                                                                       tGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAATTTGTATTGTTTCAATaaaaa                  >  8:114192/1‑138 (MQ=255)
                                                                                                       tGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAAATTGTATTGTTTCAATaaaaa                  >  8:312777/1‑138 (MQ=255)
                                                                                                       tGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAAATTGTATTGTTTCAATaaaaa                  >  4:155533/1‑138 (MQ=255)
                                                                                                       tGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAAATTGTATTGTTTCAATaaaaa                  >  3:424560/1‑138 (MQ=255)
                                                                                                       tGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAAATTGTATTGTTTCAATaaaaa                  >  1:325620/1‑138 (MQ=255)
                                                                                                       tGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAAATTGTATTGTTTCAATaaaaa                  >  1:408996/1‑138 (MQ=255)
                                                                                                       tGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAAATTGTATTGTTTCAATaaaaa                  >  4:331008/1‑138 (MQ=255)
                                                                                                       tGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAAATTGTATTGTTTCAATaaaaa                  >  7:377017/1‑138 (MQ=255)
                                                                                                       tGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAAATTGTATTGTTTCAATaaaaa                  >  3:71179/1‑138 (MQ=255)
                                                                                                       tGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAAATTGTATTGTTTCAATaaaaa                  >  2:62448/1‑138 (MQ=255)
                                                                                                       tGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAAATTGTATTGTTTCAATaaaaa                  >  1:83634/1‑138 (MQ=255)
                                                                                                       tGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAAATTGTATTGTTTCAATaaaaa                  >  5:138477/1‑138 (MQ=255)
                                                                                                       tGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAAATTGTATTGTTTCAATaaaaa                  >  7:131636/1‑138 (MQ=255)
                                                                                                       tGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAAATTGTATTGTTTCAATaaaaa                  >  5:389816/1‑138 (MQ=255)
                                                                                                       tGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAAATTGTATTGTTTCAATaaaaa                  >  1:109896/1‑138 (MQ=255)
                                                                                                       tGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAAATTGTATTGTTTCAATaaaaa                  >  3:319919/1‑138 (MQ=255)
                                                                                                       tGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAAATTGTATTGTTTCAATaaaaa                  >  6:187284/1‑138 (MQ=255)
                                                                                                       tGTTAAATCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATAAGGAAATTGTATTGTTTCAATaaaaa                  >  2:109273/1‑138 (MQ=255)
                                                                                                         ttaaaTCGACTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAAATTGTATTGTTTCAATaaaaaaa                >  3:202169/3‑138 (MQ=255)
                                                                                                                 aCTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTc                                                                >  2:381119/1‑82 (MQ=255)
                                                                                                                 aCTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTc                                                                <  1:381119/82‑1 (MQ=255)
                                                                                                                 aCTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTa                                              >  1:54123/1‑100 (MQ=255)
                                                                                                                 aCTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTa                                              <  2:54123/100‑1 (MQ=255)
                                                                                                                 aCTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAAATTGTATTGTTTCAATAAAAAAATc              >  3:371848/1‑132 (MQ=255)
                                                                                                                 aCTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAAATTGTATTGTTTCAATAAAAAAATc              <  4:371848/132‑1 (MQ=255)
                                                                                                                 aCTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAAATTGTATTGTTTCAATAAAAAAATCAGCtttt       >  2:304753/1‑138 (MQ=255)
                                                                                                                 aCTCCGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCAAAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAAATTGTATTGTTTCAATAAAAAAATCAGCttt        >  4:317645/1‑138 (MQ=255)
                                                                                                                     cGGTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAAATTGTATTGTTTCAATAAAAAAATCAGCTTTTCAGc   <  8:377017/139‑1 (MQ=255)
                                                                                                                      ggTTAATTGTTTCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAAATTGTATTGTTTCAATAAAAAAATCAGCTTTTCAGCt  <  3:84557/139‑1 (MQ=255)
                                                                                                                                 tCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAAATTGTATTGTTTCAATAAAAAAATCAGCtttt       >  7:143722/1‑122 (MQ=255)
                                                                                                                                 tCCATTTCGTTGTAAATACGCTGTATGTATCCA‑‑TCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAAATTGTATTGTTTCAATAAAAAAATCAGCtttt       <  8:143722/123‑2 (MQ=255)
                                                                                                                                     |                                                                                                                             
CTTGATCCTTTTCTAATGCCGCCGCAAAGTCCTTAAAGCTTTGAAACCCCTTTTTCTCATAAAGCGACGTCAGCCTTCTCTTCATTTGTGCTTCTTTATATAATGTCAAATCGACTCCGGTTAATTGTTTCCACTTCGTTGTAAATACGCTGTATGTATCCATCTCTCTACTCTCCCAGTCTCGCTGTCTTTTTGTCCTAGTATAGCTTGTTTTATCAGGAAATTGTATTGTTTCAATAAAAAAATCAGCTTTGCAGCT  >  CP006881/2198820‑2199078

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