Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,573,247 T→C 6.9% T700A (ACG→GCG)  PP_3155 ← ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,573,2470TC 6.9% 208.0 / 3.9 72T700A (ACG→GCG) PP_3155ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base T (36/31);  new base C (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 = 9.67e-01

CATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAGCGCTCCATCGT  >  NC_002947/3573102‑3573394
                                                                                                                                                 |                                                                                                                                                   
cATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTgcg                                                                                                                                                  >  8:60250/1‑149 (MQ=255)
 aTGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTgcgc                                                                                                                                                 >  1:189932/1‑149 (MQ=255)
 aTGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTgcgc                                                                                                                                                 <  7:60250/149‑1 (MQ=255)
     actactCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGCGCACGGTg                                                                                                                                             >  2:54391/1‑149 (MQ=255)
     actactCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGCACACGGTg                                                                                                                                             >  4:276152/1‑149 (MQ=255)
       tactCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGt                                                                                                                                                     >  5:102538/1‑139 (MQ=255)
       tactCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGt                                                                                                                                                     <  6:102538/139‑1 (MQ=255)
       tactCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTgcg                                                                                                                                           >  2:344451/1‑149 (MQ=255)
         ctCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGCGCACGGTgcagg                                                                                                                                         <  7:353130/149‑4 (MQ=255)
         ctCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGCGCACGGTgcagg                                                                                                                                         <  1:54391/149‑4 (MQ=255)
         ctCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGCACACGGTgcagg                                                                                                                                         <  3:276152/149‑4 (MQ=255)
          tCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCg                                                                                                                                        >  5:186422/1‑149 (MQ=255)
            cAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCggg                                                                                                                                      <  4:266659/149‑1 (MQ=255)
            cAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCggg                                                                                                                                      >  3:266659/1‑149 (MQ=255)
            cAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACAGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCg                                                                                                                                        <  2:290830/147‑1 (MQ=255)
             aaCGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCgggg                                                                                                                                     >  8:86262/1‑149 (MQ=255)
                          agcagcCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGt                                                                                                                        <  7:242055/149‑1 (MQ=255)
                           gcagcCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGtt                                                                                                                       <  6:143970/149‑1 (MQ=255)
                                cAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAAcg                                                                                                                  >  3:137148/1‑149 (MQ=255)
                                         gCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGt                                                                                                         <  8:296863/149‑1 (MQ=255)
                                            tAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCa                                                                                                       >  8:2950/1‑148 (MQ=255)
                                            tAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCa                                                                                                       >  2:59602/1‑148 (MQ=255)
                                                 tCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATAcccc                                                                                                 <  6:186422/149‑1 (MQ=255)
                                                    agaacGCTGTGCTGTAGAACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCttg                                                                                              >  5:403498/5‑149 (MQ=255)
                                                         gCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGaa                                                                                         <  1:344451/149‑1 (MQ=255)
                                                           tGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAaca                                                                                       >  5:417997/1‑149 (MQ=255)
                                                                tACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCt                                                                                                <  4:182675/135‑1 (MQ=255)
                                                                tACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCt                                                                                                >  3:182675/1‑135 (MQ=255)
                                                                tACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGttt                                                                                  >  1:240836/1‑149 (MQ=255)
                                                                tACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGttt                                                                                  >  2:371900/1‑149 (MQ=255)
                                                                  cagcacagcGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAt                                                                                                      <  7:346426/127‑1 (MQ=255)
                                                                  cagcacagcGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAt                                                                                                      >  8:346426/1‑127 (MQ=255)
                                                                  cagcacagcGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGttttt                                                                                <  2:240836/149‑1 (MQ=255)
                                                                  cagcacagcGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGttt                                                                                  <  5:339012/147‑1 (MQ=255)
                                                                     cacagcGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGttt                                                                                  >  2:49859/1‑144 (MQ=255)
                                                                     cacagcGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGttt                                                                                  <  1:49859/144‑1 (MQ=255)
                                                                     cacagcGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCAc                                                                             <  7:263520/149‑1 (MQ=255)
                                                                        agcGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTc                                                                          >  1:379462/1‑149 (MQ=255)
                                                                         gcGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTcc                                                                         <  7:2950/149‑1 (MQ=255)
                                                                          cGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCa                                                                        <  7:86262/149‑1 (MQ=255)
                                                                          cGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCa                                                                        <  6:403498/149‑1 (MQ=255)
                                                                            ttCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGa                                                                      <  1:371900/149‑1 (MQ=255)
                                                                             tCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGAt                                                                     >  2:128453/1‑149 (MQ=255)
                                                                               tGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGAtt                                                                    >  8:235107/1‑148 (MQ=255)
                                                                                      aCGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGc                                                             <  1:59602/148‑1 (MQ=255)
                                                                                           gCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACgtgt                                                       <  2:189932/149‑1 (MQ=255)
                                                                                             cACTGCGCGGCCCCGGCAATTCAGAAGGAATAGGAAACGCCGGCCTGCACCGTGCGCGCTGCGCCGGGGTATACATAGATGTTGAACGCGCCCTCGACATACACCTTGATCAACACGTTATAAACCTCCAGATTCAACCCCACgtctac                                                     >  3:333860/1‑145 (MQ=255)
                                                                                              aCTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCttgttg                                                                                           <  3:120651/110‑1 (MQ=255)
                                                                                              aCTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCttgttg                                                                                           >  4:120651/1‑110 (MQ=255)
                                                                                                      gCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCtt                                            >  3:45822/1‑149 (MQ=255)
                                                                                                             cAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACCCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAg                                     <  1:128453/149‑1 (MQ=255)
                                                                                                                ttCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCt                                   >  7:301476/1‑148 (MQ=255)
                                                                                                                  cAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCaa                                >  1:25235/1‑149 (MQ=255)
                                                                                                                   aGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGaa                                      <  2:175826/142‑1 (MQ=255)
                                                                                                                   aGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGaa                                      >  1:175826/1‑142 (MQ=255)
                                                                                                                    gAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCaa                                >  6:69126/1‑147 (MQ=255)
                                                                                                                     aaGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAg                             >  1:381579/1‑149 (MQ=255)
                                                                                                                       ggTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCa                                                                        >  2:244077/1‑104 (MQ=255)
                                                                                                                       ggTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCa                                                                        <  1:244077/104‑1 (MQ=255)
                                                                                                                          ataGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGa                        >  3:408014/1‑149 (MQ=255)
                                                                                                                              gAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTc                          <  3:163651/143‑1 (MQ=255)
                                                                                                                              gAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTc                          >  4:163651/1‑143 (MQ=255)
                                                                                                                              gAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCAc                    >  4:162721/1‑149 (MQ=255)
                                                                                                                              gAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCAc                    <  7:235107/149‑1 (MQ=255)
                                                                                                                              gAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGATCAACAGGTCGACAAc                    >  6:221101/1‑149 (MQ=255)
                                                                                                                              gAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATACCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACCTATTTCACGTCCAGATTCAACCGCACCTGATCCTTCAACTTGTAGAAGCTCAACAGGTCGACCAc                    >  8:265585/1‑149 (MQ=255)
                                                                                                                                   gCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGt               <  3:162721/149‑1 (MQ=255)
                                                                                                                                     cGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTg                                           >  7:4694/1‑119 (MQ=255)
                                                                                                                                     cGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTg                                           <  8:4694/119‑1 (MQ=255)
                                                                                                                                      ggCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAAGACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAgc            <  7:265585/149‑1 (MQ=255)
                                                                                                                                        ccTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAgcgc          <  4:45822/149‑1 (MQ=255)
                                                                                                                                           gCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAGCGCTc        >  8:122468/1‑148 (MQ=255)
                                                                                                                                              ccGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAGCGCTCCATc    <  6:417997/149‑1 (MQ=255)
                                                                                                                                                gTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGAGCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAGCGCTCCATCgt  >  7:72698/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAGCGCTCCATCGT  >  NC_002947/3573102‑3573394

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