Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 324,136 A→T 5.5% V35V (GTT→GTA)  PP_0267 ← ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947324,1360AT 5.5% 232.7 / 2.5 73V35V (GTT→GTA) PP_0267ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base A (35/34);  new base T (2/2);  total (37/36)
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.61e-01

AGCGTGTTGCTGCCGTAATCCAGCACATTGGAGGTTTTCACCGTATCGATGGCCTGCGGCACATAACGCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAGGGGTGAAAGACGGCAGGGTACAAGCATGGGCACATCCGGGACGGGCAGTTGAAAATAAAGTGCGAATG  >  NC_002947/323990‑324281
                                                                                                                                                  |                                                                                                                                                 
aGCGTGTTGCTGCCGTAATCCAGCACATTGGAGGTTTTCACCGTATCGATGGCCTGCGGCACATAACGCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAAc                                                                                                                                                 >  3:49772/1‑149 (MQ=255)
 gCGTGTTGCTGCCGTAATCCAGCACATTGGAGGTTTTCACCGTATCGATGGCCTGCGGCACATAACGCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAAcc                                                                                                                                                <  3:345611/149‑1 (MQ=255)
          tgcCGTAATCCAGCACATTGGAGGTTTTCACCGTATCGATGGCCTGCGGCACATAACGCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGCCTcc                                                                                                                                       >  3:264318/1‑149 (MQ=255)
           gcCGTAATCCAGCACATTGGAGGTTTTCACCGTATCGATGGCCTGCGGCACATAACGCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGTCGGTTACCgccgacag                                                                                                                                        >  6:154456/1‑145 (MQ=255)
           gcCGTAATCCAGCACATTGGAGGTTTTCACCGTATCGATGGCCTGCGGCACATAACGCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGTCGGTTACCgccgacag                                                                                                                                        <  5:154456/147‑3 (MQ=255)
               tAATCCAGCACATTGGAGGTTTTCACCGTATCGATGGCCTGCGGCACATAACGCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAg                                                                                                                                  >  4:143305/1‑149 (MQ=255)
                  tCCAGCACATTGGAGGTTTTCACCGTATCGATGGCCTGCGGCACATAACGCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAgtgt                                                                                                                               >  3:227619/1‑149 (MQ=255)
                    cAGCACATTGGAGGTTTTCACCGTATCGATGGCCTGCGGCACATAACGCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTgtgg                                                                                                                             >  8:74896/1‑149 (MQ=255)
                    cAGCACATTGGAGGTTTTCACCGTATCGATGGCCTGCGGCACATAACGCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTgtgg                                                                                                                             >  4:306988/1‑149 (MQ=255)
                    cAGCACATTGGAGGTTTTCACCGTATCGATGGCCTGCGGCACATAACGCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTgtgg                                                                                                                             >  3:298283/1‑149 (MQ=255)
                     aGCACATTGGAGGTTTTCACCGTATCGATGGCCTGCGGCACATAACGCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTgtggt                                                                                                                            <  7:74896/149‑1 (MQ=255)
                            tGGAGGTTTTCACCGTATCGATGGCCTGCGGCACATAACGCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAgtgt                                                                                                                               >  4:241304/1‑139 (MQ=255)
                            tGGAGGTTTTCACCGTATCGATGGCCTGCGGCACATAACGCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAgtgt                                                                                                                               <  3:241304/139‑1 (MQ=255)
                             ggAGGTTTTCACCGTATCGATGGCCTGCGGCACATAACGCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTggcggcg                                                                                                                     <  4:49772/148‑1 (MQ=255)
                             ggAGGTTTTCACCGTATCGATGGCCTGCGGCACATAACGCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTggcggcg                                                                                                                     <  8:121000/148‑1 (MQ=255)
                               aGGTTTTCACCGTATCGATGGCCTGCGGCACATAACGCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCa                                                                                                                  <  8:41341/149‑1 (MQ=255)
                                  ttttCACCGTATCGATGGCCTGCGGCACATAACGCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTg                                                                                                               <  4:298283/149‑1 (MQ=255)
                                    ttCACCGTATCGATGGCCTGCGGCACATAACGCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGc                                                                                                             >  1:162276/1‑149 (MQ=255)
                                     tCACCGTATCGATGGCCTGCGGCACATAACGCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGc                                                                                                             <  3:143305/148‑1 (MQ=255)
                                        ccGTATCGATGGCCTGCGGCACATAACGCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAgg                                                                                                         <  2:162276/149‑1 (MQ=255)
                                              cGATGGCCTGCGGCACATAACGCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCgg                                                                                                                                     <  3:291427/115‑1 (MQ=255)
                                              cGATGGCCTGCGGCACATAACGCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCgg                                                                                                                                     >  4:291427/1‑115 (MQ=255)
                                              cGATGGCCTGCGGCACATAACGCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGgcg                                                                                                   >  4:263750/1‑149 (MQ=255)
                                                aTGGCCTGCGGCACATAACGCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCt                                                                                                 >  6:206816/1‑149 (MQ=255)
                                                     cTGCGGCACATACCGCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAAca                                                                                            <  3:226992/149‑1 (MQ=255)
                                                        cGGCACATAACGCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAAcagca                                                                                         <  3:234666/149‑1 (MQ=255)
                                                          gCACATAACGCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACagcagc                                                                                       <  7:49370/149‑1 (MQ=255)
                                                           cacaTAACGCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACagcagc                                                                                       >  4:85894/1‑148 (MQ=255)
                                                               tAACGCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAg                                                                                  <  5:206816/149‑1 (MQ=255)
                                                                   gCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGcc                                                                               >  2:126631/1‑148 (MQ=255)
                                                                    cACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGa                                                                             <  7:206706/149‑1 (MQ=255)
                                                                      cGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGc                                                                           <  8:200403/149‑1 (MQ=255)
                                                                         gCGTCCAGGTGCGGGTTGCGCTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGATAACCGACGACTCCGGCAGTGTGCTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGACGCGACAGGCCGAGCGt                                                                         >  8:63701/1‑148 (MQ=255)
                                                                             cgtggTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAggg                                                                    <  1:126631/146‑1 (MQ=255)
                                                                              gAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAggg                                                                    <  3:85894/148‑1 (MQ=255)
                                                                                    cGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAGGGGTGa                                                                >  5:338618/1‑146 (MQ=255)
                                                                                     gggTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAGGGGTGaaa                                                              <  6:363683/147‑1 (MQ=255)
                                                                                     gggTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAGGGGTGaaa                                                              >  5:363683/1‑147 (MQ=255)
                                                                                      ggTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAGGGGTGAAAGAc                                                           >  6:106758/1‑149 (MQ=255)
                                                                                           cGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAGGGGTGAAAGACGGCa                                                       <  3:263750/148‑1 (MQ=255)
                                                                                            ggTGCTCACCTCCTTCTCCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAGGGGTGAAAGACGGCAgg                                                     <  5:193493/149‑1 (MQ=255)
                                                                                            ggTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAGGGGTGAAAGACGGCAgg                                                     <  3:306988/149‑1 (MQ=255)
                                                                                                cTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGCGCGTCAGGCGTGAAAGACGGCAGGGTa                                                  >  7:281767/1‑148 (MQ=255)
                                                                                                 tCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCAGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAGGGGTGAAAGACGGCAGGGTAc                                                 >  2:52136/1‑148 (MQ=255)
                                                                                                  cACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAGGGGTGAAAGACGGCAGGGTACaa                                               >  6:214329/1‑149 (MQ=255)
                                                                                                   aCCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAgg                                                                                 >  4:282917/1‑114 (MQ=255)
                                                                                                   aCCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAgg                                                                                 <  3:282917/114‑1 (MQ=255)
                                                                                                       cctTCACCCGCGGGTTGGCCTTCTCGCGTTCGCTGTCAGCGGTAACCGACTACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGAAGCGACAGGCCTTGCGTCAGGGGTGAAAGACGGCAGGGTACAAGCATg                                          <  1:52136/149‑1 (MQ=255)
                                                                                                        ctTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGcc                                                                               >  6:12297/1‑111 (MQ=255)
                                                                                                        ctTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGcc                                                                               <  5:12297/111‑1 (MQ=255)
                                                                                                        ctTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAGGGGTGAAAGACGGCAGGGTACAAGCATgg                                         <  6:338618/149‑1 (MQ=255)
                                                                                                         ttCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTTGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAGGGGTGAAAGACGGCAGGGTACAAGCATggg                                        >  5:212387/1‑149 (MQ=255)
                                                                                                             cccGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAGGGGTGAAAGACGGCAGGGTACAAGCATGGGcaca                                    >  1:254612/1‑149 (MQ=255)
                                                                                                              ccGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGACAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAGGCGTGAAAGACGTCAGGGTACAAGCATGGGCACAt                                   >  5:182227/1‑149 (MQ=255)
                                                                                                                    gTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAGGGGTGAAAGACGGCAGGGTACAAGCATGGGCACATCCGGGa                             >  2:209290/1‑149 (MQ=255)
                                                                                                                     ttGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAGGGGTGAAAGACGGCAGGGTACAAGCATGGGCACATCCGGGAc                            <  5:214329/149‑1 (MQ=255)
                                                                                                                     ttGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAGGGGTGAAAGACGGCAGGGTACAAGCATGGGCACATCCGGGAc                            >  4:50524/1‑149 (MQ=255)
                                                                                                                        tCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAGGGGTGAAAGACGGCAGGGTACTAGCATGGGCACATCCGGGACgg                          <  2:254612/148‑1 (MQ=255)
                                                                                                                         ccTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAGGGGTGAAAGACGGCAGGGTACaa                                               >  5:177426/1‑126 (MQ=255)
                                                                                                                         ccTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAGGGGTGAAAGACGGCAGGGTACaa                                               <  6:177426/126‑1 (MQ=255)
                                                                                                                         ccTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAGGGGTGAAAGACGGCAGGGTACAAGCATGGGCACATCCGGGACggg                         >  7:36994/1‑148 (MQ=255)
                                                                                                                          cTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAGGGGTGAAAGACGGCAGGGTACAAGCATGGGCACATCCGGGACGGGCa                       <  8:36994/149‑1 (MQ=255)
                                                                                                                           ttCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAGGGGTGAAAGACGGCAGGGTACAAGCATGGGCACATCCGGGACGGGCAg                      >  3:146308/1‑149 (MQ=255)
                                                                                                                                gcgTTCGCTGTCGGCGGTAACCGACGACTCCGGCATTGTTGTGGCGGGGGCACTGGCCAGGCTTGCGCTGAAGAGCAGCGACAGGCCGAGCGTCAGGGGTGATAGACGGCAGGGTACATGCATGGGCACATCCGGGACGGGCAGTTGa                  <  7:63701/148‑1 (MQ=255)
                                                                                                                                 cgTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGaa                                                                                              <  3:218008/71‑1 (MQ=255)
                                                                                                                                 cgTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGaa                                                                                              >  4:218008/1‑71 (MQ=255)
                                                                                                                                     cggagtcgtcGGTTACCGCCGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAGGGGTGAAAGACGGCAGGGTACAAGCATGGGCACATCCGGGACGGGCAGTTGAAAATa             >  4:122398/5‑148 (MQ=255)
                                                                                                                                       cTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAGGGGTGAAAGACGGCAGGGTACAAGCATGGGCACATCCGGGACGGGCAGTTGAAAATAAAg          <  5:106758/149‑1 (MQ=255)
                                                                                                                                        tGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAGGGGTGAAAGACGGCAGGGTACAAGCATGGGCACATCCGGGACGGGCAGTTGAAAATAAAGt         >  5:16312/1‑149 (MQ=255)
                                                                                                                                        tGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAGGGGTGAAAGACGGCAGGGTACAAGCATGGGCACATCCGGGACGGGCAGTTGAAAATAAAGt         >  2:319283/1‑149 (MQ=255)
                                                                                                                                            gtcggTTACCGCCGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAGGGGTGAAAGACGGCAGGGTACAAGCATGGGCACATCCGGGACGGGCAGTTGAAAATAAAGTGCGa     <  3:122398/147‑1 (MQ=255)
                                                                                                                                               ggTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAGGGGTGAAAGACGGCAGGGTACAAGCATGGGCACATCCGGGACGGGCAGTTGAAAATAAAGTGCGAATg  >  6:343434/1‑149 (MQ=255)
                                                                                                                                               ggTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAGGGGTGAAAGACGGCAGGGTACAAGCATGGGCACATCCGGGACGGGCAGTTGAAAATAAAGTGCGAATg  <  5:343434/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AGCGTGTTGCTGCCGTAATCCAGCACATTGGAGGTTTTCACCGTATCGATGGCCTGCGGCACATAACGCACGGGCGTCGAGGTGCGGGTTGCGGTGCTCACCTCCTTCACCCGCGGGTTGTCCTTCTCGCGTTCGCTGTCGGCGGTAACCGACGACTCCGGCAGTGTGGTGGCGGCGGCACTGGCCAGGCTGGCGCTGAACAGCAGCGACAGGCCGAGCGTCAGGGGTGAAAGACGGCAGGGTACAAGCATGGGCACATCCGGGACGGGCAGTTGAAAATAAAGTGCGAATG  >  NC_002947/323990‑324281

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