Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,834,352 C→G 6.3% V339V (GTC→GTG vgrG‑II → type VI secretion system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,834,3520CG 6.3% 204.7 / 2.5 63V339V (GTC→GTGvgrG‑IItype VI secretion system protein
Reads supporting (aligned to +/- strand):  ref base C (27/32);  new base G (2/2);  total (29/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 = 9.91e-01

TTCACTGCCACCCGCCTCACGGCCCTGCGCCATGATGCACGCATTGCCGAAGTGGAAGGCGACGATGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAGGCACCCGCTACAAACAGGAAGCATTGCTGGTTCCCGGCAGGACCGAATGGCGCCCGGCCCCCCTGCCCAAGCCAC  >  NC_002947/3834209‑3834491
                                                                                                                                               |                                                                                                                                           
ttCACTGCCACCCGCCTCACGGCCCTGCGCCATGATGCACGCATTGCCGAAGTGGAAGGCGACGATGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGt                                                                                                                                        <  7:307741/149‑1 (MQ=255)
ttCACTGCCACCCGCCTCACGGCCCTGCGCCATGATGCACGCATTGCCGAAGTGGAAGGCGACGATGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGt                                                                                                                                        <  7:123288/149‑1 (MQ=255)
   aCTGCCACCCGCCTCACGGCCCTGCGCCATGATGCACGCATTGCCGAAGTGGAAGGCGACGATGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGc                                                                                                                                      <  2:130531/148‑1 (MQ=255)
   aCTGCCACCCGCCTCACGGCCCTGCGCCATGATGCACGCATTGCCGAAGTGGAAGGCGACGATGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGc                                                                                                                                      <  1:399200/148‑1 (MQ=255)
     tGCCACCCGCCTCACGGCCCTGCGCCATGATGCACGCATTGCCGAAGTGGAAGGCGACGATGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCAt                                                                                                                                    >  2:308661/1‑148 (MQ=255)
         aCCCGCCTCACGGCCCTGCGCCATGATGCACGCATTGCCGAAGTGGAAGGCGACGATGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACg                                                                                                                                >  6:9158/1‑148 (MQ=255)
          cccGCCTCACGGCCCTGCGCCATGATGCACGCATTGCCGAAGTGGAAGGCGACGATGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGaa                                                                                                                              <  8:100583/149‑1 (MQ=255)
          cccGCCTCACGGCCCTGCGCCATGATGCACGCATTGCCGAAGTGGAAGGCGACGATGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGaa                                                                                                                              <  4:159376/149‑1 (MQ=255)
           ccGCCTCACGGCCCTGCGCCATGATGCACGCATTGCCGAAGTGGAAGGCGACGATGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGagg                                                                                                                             >  5:357623/1‑147 (MQ=255)
           ccGCCTCACGGCCCTGCGCCATGATGCACGCATTGCCGAAGTGGAAGGCGACGATGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGCTGACCGTGcgcca                                                                                                                                  <  3:257817/144‑5 (MQ=255)
           ccGCCTCACGGCCCTGCGCCATGATGCACGCATTGCCGAAGTGGAAGGCGACGATGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGCTGACCGTGcgcca                                                                                                                                  >  4:257817/1‑140 (MQ=255)
                  aCGGCCCTGCGCCATGATGCACGCATTGCCGAAGTGGAAGGCGACGATGTGCGCCTGCACCCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCAcc                                                                                                                        >  3:323425/1‑147 (MQ=255)
                  aCGGCCCTGCGCCATGATGCACGCATTGCCGAAGTGGAAGGCGACGATGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCTACGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCAcc                                                                                                                        <  5:285752/147‑1 (MQ=255)
                  aCGGCCCTGCGCCATGATGCACGCATTGCCGAAGTGGAAGGCGACGATGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCAcc                                                                                                                        <  1:226142/147‑1 (MQ=255)
                  aCGGCCCTGCGCCATGATGCACGCATTGCCGAAGTGGAAGGCGACGATGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCAcc                                                                                                                        >  3:54894/1‑147 (MQ=255)
                   cGGCCCTGCGCCATGATGCACGCATTGCCGAAGTGGAAGGCGACGATGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAg                                                                                                                     >  3:388523/1‑149 (MQ=255)
                    ggCCCTGCGCCATGATGCACGCATTGCCGAAGTGGAAGGCGACGATGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGt                                                                                                                    >  2:206601/1‑149 (MQ=255)
                       ccTGCGCCATGATGCACGCATTGCCGAAGTGGAAGGCGACGATGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCa                                                                                                                 >  3:57446/1‑149 (MQ=255)
                        cTGCGCCATGATGCACGCATTGCCGAAGTGGAAGGCGACGATGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCAc                                                                                                                <  1:206601/149‑1 (MQ=255)
                              cATGATGCACGCATTGCCGAAGTGGAAGGCGACGATGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCt                                                                                                          >  1:269231/1‑149 (MQ=255)
                                tgatgCACGCATTGCCGAAGTGGAAGGCGACGATGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCa                                                                                                                          >  1:387402/1‑131 (MQ=255)
                                tgatgCACGCATTGCCGAAGTGGAAGGCGACGATGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCa                                                                                                                          <  2:387402/131‑1 (MQ=255)
                                  atgCACGCATTGCCGAAGTGGAAGGCGACGATGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCa                                                                                                       <  1:127615/148‑1 (MQ=255)
                                  atgCACGCATTGCCGAAGTGGAAGGCGACGATGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAg                                                                                                      <  7:73586/149‑1 (MQ=255)
                                              ccGAAGTGGATGGCGACGATGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGgccgcc                                                                                          >  1:240982/1‑149 (MQ=255)
                                              ccGAAGTGGAAGGCGACGATGTGCGCCTGCAACCGGGGCAGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGgccgcc                                                                                          >  4:88457/1‑149 (MQ=255)
                                                       aaGGCGACGATGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACAGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCgcag                                                                                 >  1:130288/1‑149 (MQ=255)
                                                        aGGCGACGATGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCgcagc                                                                                <  2:269231/149‑1 (MQ=255)
                                                        aGGCGACGATGTGCGCCTGAAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAATACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCGTCAGGAAGAGGCCGCCGACTAgcag                                                                                 >  4:403895/1‑148 (MQ=255)
                                                           cgacgaTGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAg                                                                             <  3:88457/149‑1 (MQ=255)
                                                           cgacgaTGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAg                                                                             <  1:378710/149‑1 (MQ=255)
                                                           cgacgaTGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAg                                                                             <  5:9158/149‑1 (MQ=255)
                                                            gacgaTGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCAc                                                                                                                         >  6:40589/1‑104 (MQ=255)
                                                            gacgaTGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCAc                                                                                                                         <  5:40589/104‑1 (MQ=255)
                                                              cgaTGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAGGCa                                                                          >  7:297918/1‑149 (MQ=255)
                                                                      gcCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAg                                                                             <  8:306058/138‑1 (MQ=255)
                                                                      gcCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAg                                                                             >  7:306058/1‑138 (MQ=255)
                                                                              ccGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAGGCACCCGCTACAAACAgg                                                           >  2:361003/1‑148 (MQ=255)
                                                                                  ggCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAGGCACCCGCTACAAACAGGAAGCa                                                      <  1:308661/149‑1 (MQ=255)
                                                                                  ggCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAGGCACCCGCTACAAACAGGAAGCa                                                      <  4:57446/149‑1 (MQ=255)
                                                                                    cTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAACAAGGCACCCGCTAAAAACAGGAAGCAtt                                                    >  1:140411/1‑149 (MQ=255)
                                                                                      gAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAGGCACCCGCTACAAAc                                                              >  5:215091/1‑137 (MQ=255)
                                                                                      gAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAGGCACCCGCTACAAAc                                                              <  6:215091/137‑1 (MQ=255)
                                                                                       aGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAGGCACCCGCTACAAACAGGAAGCATTGCt                                                 <  6:357623/149‑1 (MQ=255)
                                                                                             aCCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAGGCACCCGCTACAAACAGGAAg                                                        >  6:170851/1‑136 (MQ=255)
                                                                                             aCCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAGGCACCCGCTACAAACAGGAAg                                                        <  5:170851/136‑1 (MQ=255)
                                                                                               ccTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAGGCACGCGCTACAAACAGGAAGCATTGCTGGTTCCCg                                         >  4:198477/1‑149 (MQ=255)
                                                                                               ccTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAGGCACCCGCTACAAACAGGAAGCATTGCTGGTTCCCg                                         <  1:361003/149‑1 (MQ=255)
                                                                                               ccTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAGGCACCCGCTACAAACAGGAAGCATTGCTGGTTCCCg                                         >  8:54706/1‑149 (MQ=255)
                                                                                               ccTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAACAAGGCACCCGCTACAAACAGGAAGCATTGCTGGTTCCCg                                         <  2:140411/149‑1 (MQ=255)
                                                                                                   aTCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCGAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAGGCACCCGCTACAAACAGGAAGCATTGCTGGTTCCCGGCAg                                     <  4:323425/149‑1 (MQ=255)
                                                                                                   aTCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAGGCACCCGCTACAAACAGGAAGCATTGCTGGTTCCCGGCAg                                     >  7:50591/1‑149 (MQ=255)
                                                                                                      gAGCACCCGCGTGAAGACCTCACCGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAGGCACCCGCTACAAACAGGAAGCATTGCTGGTTCCCGGCAGGa                                   <  8:376748/148‑1 (MQ=255)
                                                                                                            ccGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAGGCACCCGCTACAAACAGGAAGCATTGCTGGTTCCCGGCAGGACCGAATg                            >  8:355305/1‑149 (MQ=255)
                                                                                                                  gAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAGGCAccc                                                                       <  1:68806/100‑1 (MQ=255)
                                                                                                                  gAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAGGCAccc                                                                       >  2:68806/1‑100 (MQ=255)
                                                                                                                   aaGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAGGCACCCGCTACAAACAGGAAGCATTGCTGGTTCCCGGCAGGACCGAATGGCGccc                      <  4:54894/148‑1 (MQ=255)
                                                                                                                    aGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAGGCACCCGCTACAAACAGGAAGCATTGCTGGTTCCCGGCAGGACCGAATGGCGCCCgg                    <  4:388523/149‑1 (MQ=255)
                                                                                                                          cAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAGGCACCCGCTACAAACAGGAAGCATTGCTGGTTCCCGGCAGGACCGAATGGCGCCCGGcccccc              <  6:80846/149‑1 (MQ=255)
                                                                                                                                 cACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAGGCACCCGCTACAAACAGGAAGCATTGCTGGTTCCCGGCAGGACCGAATGGCGCCCGGCCCCCCTGccc         >  2:220335/1‑147 (MQ=255)
                                                                                                                                    tgatgcACGCTGACCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAGGCACCCGCTACAAACAGGAAGCATTGCTGGTTCCCGGCAGGACCGAATGGCGCCCGGCCCCCCTGc           <  1:147797/138‑1 (MQ=255)
                                                                                                                                    tgatgcACGCTGACCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAGGCACCCGCTACAAACAGGAAGCATTGCTGGTTCCCGGCAGGACCGAATGGCGCCCGGCCCCCCTGc           >  2:147797/5‑142 (MQ=255)
                                                                                                                                      gcgcACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAGGCACCCGCTACAAACAGGAAGCATTGCTGGTTCCCGGCAGGACCGAATGGCGCCCGGCCCCCCTGCCCAAGCCAc  <  2:240982/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                           
TTCACTGCCACCCGCCTCACGGCCCTGCGCCATGATGCACGCATTGCCGAAGTGGAAGGCGACGATGTGCGCCTGCAACCGGGGCTGAGCTTCACCCTCATCGAGCACCCGCGTGAAGACCTCAACGTGCACTGGCGCACGGTCAGCGTGCATCACGAAGGCACCCAGTTCACCAGCCTGCAGGAAGAGGCCGCCGACTCGCAGCAAGGCACCCGCTACAAACAGGAAGCATTGCTGGTTCCCGGCAGGACCGAATGGCGCCCGGCCCCCCTGCCCAAGCCAC  >  NC_002947/3834209‑3834491

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