Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 1,111,634 G→A 6.8% S38L (TCG→TTG)  U712_05455 ← putative spore germination protein gerPF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068811,111,6340GA 6.8% 201.7 / 3.5 72S38L (TCG→TTG) U712_05455putative spore germination protein gerPF
Reads supporting (aligned to +/- strand):  ref base G (36/31);  new base A (2/3);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

TCTTTGCAGAGCGGCTGCTTTTATCCGTTTCCAACCATATCCTGATCGTTGACATCGGGGTCTATATAATTTGTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTTACGACTCCGCCCGATATGGAATTGATAGCTATAGGCCCGACAATTGCTGGCATAAACGACACCCTTTATCGAAGTCT  >  CP006881/1111508‑1111770
                                                                                                                              |                                                                                                                                        
tctTTGCAGAGCGGCTGCTTTTATCCGTTTCCAACCATATCCTGATCGTTGACATCGGGGTCTATATAATTTGTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAgcgc                                                                                                                                <  8:242490/137‑1 (MQ=255)
     gCAGAGCGGCTGCTTTTATCCGTTTCCAACCATATCCTGATCGTTGACATCGGGGTCTATATAATTTGTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGAtttt                                                                                                                         >  1:47132/1‑139 (MQ=255)
       agagCGGCTGCTTTTATCCGTTTCCAACCATATCCTGATCGTTGACATCGGGGTCTATATAATTTGTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGCACCGACCGCAGATTTTg                                                                                                                        >  7:190949/1‑138 (MQ=255)
          gCGGCTGCTTTTATCCGTTTCCAACCATATCCTGATCGTTGACATCGGGGTCTATATAATTTGTCGCCTT‑ACTGCATTATTTAATAGAACGAAATCTCCCGTGTTTCCTGCACCCGAACCGACCGCAGATTTTGAagag                                                                                                                   >  7:47289/1‑139 (MQ=255)
                     tATCCGTTTCCAACCATATCCTGATCGTTGACATCGGGGTCTATATAATTTGTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGGTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCg                                                                                                            >  5:199317/1‑136 (MQ=255)
                      aTCCGTTTCCAACCATATCCTGATCGTTGACATCGGGGTCTATTTAATTTGTCGCGTTGACTGCATTATTTAATAGTAGGAAATCTCCCGTGTTTCCTTCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGgag                                                                                                         <  5:274967/138‑1 (MQ=255)
                       tCCGTTTCCAACCATATCCTGATCGTTGACATCGGGTTCTATATAATTTGTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGgaga                                                                                                        <  1:320497/138‑1 (MQ=255)
                       tCCGTTTCCAACCATATCCTGATCGTTGACATCGGGGTCTATATAATTTGTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGgaga                                                                                                        <  3:265295/138‑1 (MQ=255)
                          gTTTCCAACCATATCCTGATCGTTGACATCGGGGTCTATATAATTTGTCGCGTTGACTGCATTATTTAAGAGAAGGTAATCTCCCGTGTTTCCTGCTCCCGACCCGAGCGCAGATTTTGAAGCGCTTTTCGGAGAAAgg                                                                                                    >  5:141873/1‑139 (MQ=255)
                           tttCCAACCATATCCTGATCGTTGACATCGGGGTCTATATAATTTGTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCACATTTTGAAGAGCTTTTCGGAGAAAGGt                                                                                                   >  5:285692/1‑139 (MQ=255)
                           tttCCAACCATATCCTGATCGTTGACATCGGGGTCTATATAATTTGTCGCGTTGACGGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGCGCGCCGATTTTGAAGAGCTTTTCGGAGAAAgg                                                                                                    >  3:17999/1‑138 (MQ=255)
                                     tatCCTGATCGTTGACATCGGGGTCTATATAATTTGTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCCCGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGaaa                                                                                          >  8:217706/1‑137 (MQ=255)
                                       tCCTGATCGTTGACATCGGGGTCTATATAATTTGTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCTAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCAc                                                                                       >  1:201142/1‑139 (MQ=255)
                                         cTGATCGTTGACATCGGGGTCTATATAATTTGTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCa                                                                                     <  1:201623/139‑1 (MQ=255)
                                         cTGATCGTTGACATCGGGGTCTATATAATTTGTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCa                                                                                     <  1:102296/139‑1 (MQ=255)
                                            aTCGTTGACATCGGGGTCTATATAATTTGTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGaa                                                                                                                      >  2:113242/1‑103 (MQ=255)
                                            aTCGTTGACATCGGGGTCTATATAATTTGTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGaa                                                                                                                      <  1:113242/103‑1 (MQ=255)
                                            aTCGTTGACATCGGGGTCTATATAATTTGTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGgaga                                                                                                        >  6:172208/1‑117 (MQ=255)
                                            aTCGTTGACATCGGGGTCTATATAATTTGTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGgaga                                                                                                        <  5:172208/117‑1 (MQ=255)
                                            aTCGTTGACATCGGGGTCTATATAATTTGTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCa                                                                                     <  8:168127/136‑1 (MQ=255)
                                            aTCGTTGACATCGGGGTCTATATAATTTGTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCa                                                                                     >  7:168127/1‑136 (MQ=255)
                                            aTCGTTGACATCGGGGTCTATATAATTTGTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAg                                                                                  >  7:71017/1‑139 (MQ=255)
                                                ttGACATCGGGGTCTATATAATTTGTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCa                                                                                     <  1:111222/132‑1 (MQ=255)
                                                ttGACATCGGGGTCTATATAATTTGTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCa                                                                                     >  2:111222/1‑132 (MQ=255)
                                                 tGACATCGGGGTCTATATAATTTGTCGCGTTGCCTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTTAc                                                                             >  2:38400/1‑139 (MQ=255)
                                                 tGACATCGGGGTCTATATAATTTGTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTAGGAGAAAGGTAAAAGGAATCACCAAAGTTTAc                                                                             >  8:11362/1‑139 (MQ=255)
                                                 tGACATCGGGGTCTATATAATTTGTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCAAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCa                                                                                     >  8:352071/1‑131 (MQ=255)
                                                 tGACATCGGGGTCTATATAATTTGTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCAAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCa                                                                                     >  8:356857/1‑131 (MQ=255)
                                                 tGACATCGGGGTCTATATAATTTGTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCAAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCa                                                                                     <  7:352071/131‑1 (MQ=255)
                                                 tGACATCGGGGTCTATATAATTTGTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCAAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCa                                                                                     <  7:356857/131‑1 (MQ=255)
                                                          ggTCTATATAATTTGTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTTACGACTCCGcc                                                                    <  7:34212/139‑1 (MQ=255)
                                                                tataATTTGTCGCGTT‑ACTGCATTATTTAATAGAAGGAAAGCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTAACGACTCCGCCCGata                                                               <  8:47289/137‑1 (MQ=255)
                                                                     tttGTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGAGTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTTACGACTCCGCCCGATATGGAAt                                                         >  5:205023/1‑139 (MQ=255)
                                                                       tGTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTTACGACTCCGCCCGATATGGAAtt                                                        >  5:184111/1‑138 (MQ=255)
                                                                        gTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGgaga                                                                                                        >  6:195085/1‑89 (MQ=255)
                                                                        gTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGgaga                                                                                                        <  5:195085/89‑1 (MQ=255)
                                                                        gTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCa                                                                                     <  1:84837/108‑1 (MQ=255)
                                                                        gTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCa                                                                                     >  2:84837/1‑108 (MQ=255)
                                                                        gTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTTACGACTCCGCCCGATATGGAAt                                                         >  4:284681/1‑136 (MQ=255)
                                                                        gTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTTACGACTCCGCCCGATATGGAAt                                                         <  3:284681/136‑1 (MQ=255)
                                                                        gTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTTACGACTCCGCCCGATATGGAATTGa                                                      >  4:200154/1‑139 (MQ=255)
                                                                         tCGCGTTTACTGCATTATTTAATAGACGCAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGTAGAAAGGTAAAAGGAATCACCAAAGTTTACGACTCCGCCCGATATGGAATTGAt                                                     <  6:205023/139‑1 (MQ=255)
                                                                         tCGCGTTGACTGCATTATTTATAAGAAGGAAATCTCCCGTTTTTCCTGCTCCCGAACCTAGGGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTTACGACTCCGCCCGATATGGAATTGAt                                                     <  2:201142/139‑1 (MQ=255)
                                                                         tCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTTACGACTCCGCCCGATATGGAATTGAt                                                     <  6:184111/139‑1 (MQ=255)
                                                                         tCGCGTTGACTGCATTATTAAATAGAAGGAAATCTCCCGTGTTACCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGTTAAAAGGAATCACCATAGTTTACGACTCCGCCCGAAATGGAATTGAt                                                     <  3:125838/139‑1 (MQ=255)
                                                                         tCGCGGTGACGGCATTATTTAATAGAACGAAATCTCCCGTGTTTCCTGCTCCCGACCGGAGCGCAGATTTTGAAGAGGTTTTCGTAGAAAGGTAAAAGGAATCACCAAAGTTTACGACTCCGCCCGATATGGAATTGAt                                                     <  4:317375/139‑1 (MQ=255)
                                                                         gcgctTTGACTTCATTATTTAATAGAAGCAAATCTCCCGTGTTTCCTGCGCCCGAACCGCGCGCAGATTTTGAAGAGCTTTTCGGACAAAGGAAAAAGCAATCACCAAAGTTTACGACTCCGCCCGATATGGAATTGAt                                                     <  1:174181/138‑1 (MQ=255)
                                                                          cgcgTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTTACGACTCCGCCCGATATGGAATTGa                                                      >  6:261443/1‑137 (MQ=255)
                                                                             ggtGACTGAATTATTTAATAGAAGGAAAGCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTTACGACTCCGCCCGATATGGAATTGATAGCt                                                 <  5:252730/137‑1 (MQ=255)
                                                                                             aaTAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGCGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTTACTACTCCGCGCGATATGGAATTGATAGCTATTGGCCCGACAATTg                                 >  2:28162/1‑139 (MQ=255)
                                                                                               tAGATGGACATCTCCTGTGGTTCCGGCTCCCAAACGGAGGGCAGTTTTTGATCCTGTTTTCGGAGAACGGTAAAAGGAATACCCAATTTTAACGACTCCGCCCGATATGGAAGTGATAGCAATAGGCCGGACATTTGCt                               <  8:216173/139‑1 (MQ=255)
                                                                                               tAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTTACGACTcc                                                                       <  2:296642/99‑1 (MQ=255)
                                                                                               tAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTTACGACTcc                                                                       >  1:296642/1‑99 (MQ=255)
                                                                                               tAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTAGGAGAAAGGTAAAAGGAATCACCAAAGTTTACGACTCCGCCCGATATGGAATTGATAGCTATAGGCCCGACAATTGCt                               >  8:213890/1‑139 (MQ=255)
                                                                                                  aaGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCa                                                                                     >  2:43232/1‑82 (MQ=255)
                                                                                                  aaGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCa                                                                                     <  1:43232/82‑1 (MQ=255)
                                                                                                  aaGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTTACGACTCCGCCCGATATg                                                             <  7:68602/106‑1 (MQ=255)
                                                                                                  aaGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTTACGACTCCGCCCGATATg                                                             >  8:68602/1‑106 (MQ=255)
                                                                                                    ggAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTTACGACTCCGCCCGATATGGAATTGATAGCTATAGGCCCGACAATTGCTGGCa                           >  5:137212/1‑138 (MQ=255)
                                                                                                    ggAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTTACGACTCCGCCCGATATGGAATTGATAGCTATAGGCCCGACAATTGCTGGCa                           >  7:9933/1‑138 (MQ=255)
                                                                                                       aaTCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTTACGACTCCGcc                                                                    <  2:4709/94‑1 (MQ=255)
                                                                                                       aaTCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTTACGACTCCGcc                                                                    >  1:4709/1‑94 (MQ=255)
                                                                                                            cccGTGTTTCCTGCTCCCGACCCGAGCGCAGATTTTGAAGAGCTATTCGGAGAAAGGTAAAAGGAATCACCAAAGTAAACGACTCCGCCCGATATGGAATTGATAGCTATAGGCCCGACAATTGCTGGCATAAACGaca                  >  1:37596/1‑139 (MQ=255)
                                                                                                            cccGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTTACGCCTCCGCCCGATATGGAATTGATAGCTATAGGCCCGACAATTGCTGGCATAAACGaca                  >  7:327229/1‑139 (MQ=255)
                                                                                                            cccGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTTACGACTCCGCCCGATATGGAATTGATAGCTATAGGCCCGACAATTGCTGGCATAAACGaca                  <  8:9933/139‑1 (MQ=255)
                                                                                                            cccGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTTACGACTCCGCCCGATATGGAATTGATAGCTATAGGCCCGACAATTGCTGGCATAAACGaca                  >  5:312476/1‑139 (MQ=255)
                                                                                                              cGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTTACGACTCCGCCCGATATGGAATTGATAGCTATAGGCCCGACAATTGCTGGCATAAACGACAcc                <  6:137212/139‑1 (MQ=255)
                                                                                                                     ccTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTTACGACTCCGCCCGATATGGAATTGATAGCTATAGGCCCGACAATTGCTGGCATAAACGACACCCTTTATc         >  3:64692/1‑139 (MQ=255)
                                                                                                                       tGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTTACGACTCCGCCCGATATGGAATTGATAGCTATAGGCCCGACAATTGCTGGCATAAACGACACCCTGTATCg        >  5:343357/1‑138 (MQ=255)
                                                                                                                          tCCCGAACCGAGCGCAGATTTTGACGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTTACGACTCCGCCCGATATGGAATTGATAGCTATAGGCCCGACAATTGCTGGCATAAACGACACCCTTTATCGAAGt    <  3:200154/139‑1 (MQ=255)
                                                                                                                            ccGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTTACGACTCCGCCCGATATGGAATTGATAGCTATAGGCCCGACAATTGCTGGCATAAACGACACCCTTTATCGAAGtct  <  6:312476/139‑1 (MQ=255)
                                                                                                                            ccGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTTACGACTCCGCCCGATATGGAATTGATAGCTATAGGCCCGACAATTGCTGGCATAAACGACACCCTTTATCGAAGtct  <  5:261443/139‑1 (MQ=255)
                                                                                                                              |                                                                                                                                        
TCTTTGCAGAGCGGCTGCTTTTATCCGTTTCCAACCATATCCTGATCGTTGACATCGGGGTCTATATAATTTGTCGCGTTGACTGCATTATTTAATAGAAGGAAATCTCCCGTGTTTCCTGCTCCCGAACCGAGCGCAGATTTTGAAGAGCTTTTCGGAGAAAGGTAAAAGGAATCACCAAAGTTTACGACTCCGCCCGATATGGAATTGATAGCTATAGGCCCGACAATTGCTGGCATAAACGACACCCTTTATCGAAGTCT  >  CP006881/1111508‑1111770

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: ‑