Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,673,549 G→C 10.1% P107P (CCG→CCC PP_4137 → outer membrane siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,673,5490GC10.1% 197.2 / 13.0 69P107P (CCG→CCCPP_4137outer membrane siderophore receptor
Reads supporting (aligned to +/- strand):  ref base G (31/31);  new base C (4/3);  total (35/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 = 5.27e-01

ACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAGGCTCCC  >  NC_002947/4673413‑4673694
                                                                                                                                        |                                                                                                                                                 
aCATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGGCGCACGGt                                                                                                                                       >  4:92568/1‑149 (MQ=255)
    gCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCg                                                                                                                                   >  2:215075/1‑149 (MQ=255)
     cAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCg                                                                                                                                   <  2:309667/148‑1 (MQ=255)
     cAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGGCGCACGGttgag                                                                                                                                   <  4:168180/148‑4 (MQ=255)
      aGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATGGGCAGGTAACGATGGCATCAAGCTGTCTATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGgt                                                                                                                                 <  2:116395/149‑1 (MQ=255)
      aGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACg                                                                                                                                         <  7:281726/141‑1 (MQ=255)
      aGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACg                                                                                                                                         >  8:281726/1‑141 (MQ=255)
      aGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGgt                                                                                                                                 >  2:117326/1‑149 (MQ=255)
       gcgcGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTcc                                                                                                                                    <  1:87903/145‑1 (MQ=255)
       gcgcGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTcc                                                                                                                                    >  2:87903/1‑145 (MQ=255)
        cgcgTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGGCGCACGGttgaggccc                                                                                                                               <  3:92568/149‑8 (MQ=255)
         gcgTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTc                                                                                                                              <  1:273103/149‑1 (MQ=255)
                    ggCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGtt                                                                                                                   <  2:197503/149‑1 (MQ=255)
                      ccGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCg                                                                                                                 <  7:301003/149‑1 (MQ=255)
                        gAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGAc                                                                                                               >  1:106923/1‑149 (MQ=255)
                             gCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCt                                                                                                          <  3:329492/149‑1 (MQ=255)
                             gCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCt                                                                                                          >  7:253689/1‑149 (MQ=255)
                                aGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTc                                                                                                                                     <  5:146185/119‑1 (MQ=255)
                                aGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTc                                                                                                                                     >  6:146185/1‑119 (MQ=255)
                                 gCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCg                                                                                                      <  1:215075/149‑1 (MQ=255)
                                    aCCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGtt                                                                                                    >  2:62606/1‑148 (MQ=255)
                                     ccAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGa                                                                                                  >  7:325096/1‑149 (MQ=255)
                                         gACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTTACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCgg                                                                                              >  1:243755/1‑149 (MQ=255)
                                              gCTGGCCTACCAGCAAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGACAACACCATGTTCGACGGCCTGCCGTTGACCGGCCCgg                                                                                         >  6:181514/1‑149 (MQ=255)
                                                 ggCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGggcg                                                                                      <  1:117326/149‑1 (MQ=255)
                                                  gCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATAAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGggcgg                                                                                     >  1:95783/1‑149 (MQ=255)
                                                   ccTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGggcggc                                                                                    <  8:325096/149‑1 (MQ=255)
                                                    cTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCa                                                                                   >  7:114179/1‑149 (MQ=255)
                                                               gggTGTTGGCCAAATCGGATGGTAAGTATTGAATAAAGCTGTCTATACGCGGCTGGGGCCTCAACGGTGCGCCGGGGGCGAACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGaa                                                                        <  2:243755/149‑1 (MQ=255)
                                                               gggTGTCCGCCCAATCGGCAGGTAACGATGGAATCAAGCTGTCGATACGCGGCTCGGGCCGCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGaa                                                                        <  2:234224/149‑1 (MQ=255)
                                                                    ttCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGtt                                                                   >  7:224281/1‑149 (MQ=255)
                                                                    ttCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGtt                                                                   <  8:224281/149‑1 (MQ=255)
                                                                      cGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTc                                                                  <  4:311321/148‑1 (MQ=255)
                                                                       gCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCg                                                                 >  3:237130/1‑148 (MQ=255)
                                                                         ccAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAAc                                                              >  3:222743/1‑149 (MQ=255)
                                                                          cAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGCGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAAcc                                                             >  5:44579/1‑149 (MQ=255)
                                                                               ggCAGTTAACGATGGCAGCAAGCTGTCGATACGCGGCGGGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCtgt                                                        <  6:35984/149‑1 (MQ=255)
                                                                               ggCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCtgt                                                        >  6:253626/1‑149 (MQ=255)
                                                                                gCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCtgtg                                                       >  7:62039/1‑149 (MQ=255)
                                                                                  aGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGc                                                     <  2:106923/149‑1 (MQ=255)
                                                                                       aCGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGc                                                >  1:201460/1‑149 (MQ=255)
                                                                                       aCGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGc                                                <  5:253626/149‑1 (MQ=255)
                                                                                           tGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGt                                             <  2:95783/148‑1 (MQ=255)
                                                                                               aTCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCg                                         >  8:197422/1‑148 (MQ=255)
                                                                                                tCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGaa                                       >  1:250520/1‑149 (MQ=255)
                                                                                                  aaGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGt                                     >  3:14765/1‑149 (MQ=255)
                                                                                                           aTACGCGGCGCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACAATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCTCTGTGTCTGAGCCGTGCCGAAGTGCTGCGTgg                            <  5:181514/149‑1 (MQ=255)
                                                                                                            tACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGgc                           <  8:62039/149‑1 (MQ=255)
                                                                                                             aCGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGgcg                          >  4:17815/1‑149 (MQ=255)
                                                                                                              cgcgGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGgcgc                         <  2:201460/149‑1 (MQ=255)
                                                                                                              cgcgGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGgcgc                         <  2:250520/149‑1 (MQ=255)
                                                                                                              cgcgGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGgcg                          <  4:222743/148‑1 (MQ=255)
                                                                                                                      gggCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCt                 >  6:119443/1‑149 (MQ=255)
                                                                                                                        gCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGTACGCCTTACGAACTGt                                                                    >  6:103555/1‑96 (MQ=255)
                                                                                                                        gCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGTACGCCTTACGAACTGt                                                                    <  5:103555/96‑1 (MQ=255)
                                                                                                                         ccTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGaa                                       >  1:137052/1‑124 (MQ=255)
                                                                                                                           tCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGtgct                                  >  6:258693/1‑127 (MQ=255)
                                                                                                                           tCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGtgct                                  <  5:258693/127‑1 (MQ=255)
                                                                                                                              aCCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCaa         <  4:237130/149‑1 (MQ=255)
                                                                                                                               ccGTGCGCCGGGGGCGCCCGGTTCCGGTGTCTACACCATGTTCCACGGCCCGCCGTTGAGCGGCCCGGGCGGGACGCCGTACGAACTGTTCGAACCGCAGTGGCTGAGCCCAGCCGAACTGCAGCGTGGCGCCAATGGCTTTGACCAAg        >  1:190686/1‑149 (MQ=255)
                                                                                                                               ccGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAg        <  4:125348/149‑1 (MQ=255)
                                                                                                                               ccGTGCGCCCCCGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAg        >  2:54624/1‑149 (MQ=255)
                                                                                                                                 gTGCGCCCCCGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAg        <  4:133347/147‑1 (MQ=255)
                                                                                                                                 gTGCGCCCCCGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAg        >  3:133347/1‑147 (MQ=255)
                                                                                                                                 gTGCGCCCCCGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAg        >  3:128294/1‑147 (MQ=255)
                                                                                                                                    cgcCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAGGCTcc   <  5:119443/149‑1 (MQ=255)
                                                                                                                                     gcCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAGGccagg  >  5:65941/1‑145 (MQ=255)
                                                                                                                                      ccGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCt                                                          >  4:319372/1‑92 (MQ=255)
                                                                                                                                      ccGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCt                                                          <  3:319372/92‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                                 
ACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAGGCTCCC  >  NC_002947/4673413‑4673694

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: