Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,578,239 T→C 5.1% 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,578,2390TC 5.1% 834.8 / 19.8 257T700A (ACG→GCG) PP_3155ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base T (117/126);  new base C (2/11);  total (120/137)
Fisher's exact test for biased strand distribution p-value = 2.32e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01

Only 100 of 258 total aligned reads displayed.

ATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAGCGCTCCATCGTGT  >  NC_002947/3578095‑3578388
                                                                                                                                                |                                                                                                                                                       
aTGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCGGGTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGACGGTATAGGAAACGCCGGCCTGCACCGTGcgcg                                                                                                                                                    >  1:1299927/1‑150 (MQ=255)
    actactCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTgc                                                                                                                                                >  2:701281/1‑150 (MQ=255)
      tactCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTgc‑‑gc                                                                                                                                            >  1:4952687/1‑150 (MQ=255)
      tactCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACAGTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGCGCACGGTgcag                                                                                                                                              <  2:1402901/150‑3 (MQ=255)
      tacgCCAACGCCGGGATAACGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGCGCACGGTgcag                                                                                                                                              <  2:1403090/150‑3 (MQ=255)
        ctCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGCGCACGGTgcaggc                                                                                                                                            <  2:3292665/150‑5 (MQ=255)
         tCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCgg                                                                                                                                         >  1:3327786/1‑150 (MQ=255)
         tCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCgg                                                                                                                                         >  2:3324623/1‑150 (MQ=255)
         tCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGCGCACGGTgcaggcc                                                                                                                                           <  1:1384386/150‑6 (MQ=255)
         accAACGCCGGGATAAAGCAGCCACACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGCGCACGGTgcaggcc                                                                                                                                           <  2:2746647/149‑6 (MQ=255)
          ccACCGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACCTGACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCggg                                                                                                                                        <  2:5294274/150‑1 (MQ=255)
          ccAACGCCGGGATAACGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGGTCTACCGCACAGCGTTCGGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCggg                                                                                                                                        <  1:3683717/150‑1 (MQ=255)
          ccAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGCGCACGGTgcaggccg                                                                                                                                          <  2:538399/150‑7 (MQ=255)
          ccAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCCGCACCGCGCACGGTgcaggccg                                                                                                                                          <  2:4340838/150‑7 (MQ=255)
          ccAACGCCGGGATAAAGCAGCAAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCggg                                                                                                                                        <  1:3839117/150‑1 (MQ=255)
          aaaaCGACGTGATAAATCAGCCAAACGACGGCATAGAGTCGTACTCGATGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAATGTATAGGAAACGCCGGCCTGCACAGTGCGCTGTGC‑‑GCCggg                                                                                                                                        <  2:3913535/148‑1 (MQ=255)
           cAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCTgggg                                                                                                                                       >  1:1133837/1‑150 (MQ=255)
           cAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCgggg                                                                                                                                       >  1:1990075/1‑150 (MQ=255)
           cAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCgggg                                                                                                                                       >  2:5054743/1‑150 (MQ=255)
            aaCGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGt                                                                                                                                      <  1:957292/150‑1 (MQ=255)
                 cGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGTGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACa                                                                                                                                 <  2:1731230/150‑1 (MQ=255)
                  gggATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACAt                                                                                                                                >  1:1799554/1‑150 (MQ=255)
                     aTAAAGCAGCCAAACGACGGCATAGAGTCGTACTCTCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGc                                                                                                                             >  1:176208/1‑150 (MQ=255)
                      tAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCt                                                                                                                            <  2:4103212/150‑1 (MQ=255)
                        aaGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGt                                                                                                                          >  1:1710525/1‑150 (MQ=255)
                        aaGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGt                                                                                                                          <  1:4074208/150‑1 (MQ=255)
                         agcagcCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCgg                                                                                                                                         <  1:347574/134‑1 (MQ=255)
                         agcagcCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGtt                                                                                                                         >  1:525561/1‑150 (MQ=255)
                          gcagcCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTg                                                                                                                        >  1:4615763/1‑150 (MQ=255)
                             gcCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAAc                                                                                                                     <  1:648191/150‑1 (MQ=255)
                                 aaCGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACgcgc                                                                                                                 <  2:273237/150‑1 (MQ=255)
                                 aaCGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACgcgc                                                                                                                 <  2:4599512/150‑1 (MQ=255)
                                         cATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCa                                                                                                         >  1:2868217/1‑150 (MQ=255)
                                              tgTCGTACTCGCTGTTCTACATCACAGCGTTCTGACATTACGGGGCCTCTGCGCGGCCCCGGATATTCAGAAGGTATAGGATACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCTTagct                                                                                                    <  2:2915629/149‑4 (MQ=255)
                                              agTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATAccc                                                                                                    <  2:3195964/150‑1 (MQ=255)
                                                 tcagaacGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCtt                                                                                                 >  2:2684405/7‑150 (MQ=255)
                                                     cTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCttgttg                                                                                             <  2:1528127/150‑1 (MQ=255)
                                                     cTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCttgttg                                                                                             >  1:3520111/1‑150 (MQ=255)
                                                     cTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCttgttg                                                                                             >  2:5124681/1‑150 (MQ=255)
                                                     cTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCttgttg                                                                                             >  1:4148294/1‑150 (MQ=255)
                                                       cGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGaa                                                                                           <  1:3607947/150‑1 (MQ=255)
                                                        gCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAac                                                                                          >  1:2237538/1‑150 (MQ=255)
                                                           gTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACg                                                                                       >  2:2397518/1‑150 (MQ=255)
                                                            ttCTACAGCACAGCGTTCTGACAGTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGt                                                                                      <  2:2322149/150‑1 (MQ=255)
                                                              cTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCTGCCTGCACCGTGCGCGGTGC‑‑GCCGTGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGttt                                                                                    <  2:1868707/150‑1 (MQ=255)
                                                               tACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATCCATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGtttt                                                                                   >  2:3490959/1‑150 (MQ=255)
                                                               tACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGtttt                                                                                   <  2:4394743/150‑1 (MQ=255)
                                                               tACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGtttt                                                                                   >  2:4160368/1‑150 (MQ=255)
                                                                acagcacagcGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGttttt                                                                                  >  1:2664049/1‑150 (MQ=255)
                                                                 cagcacagcGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTc                                                                                 <  1:3738127/150‑1 (MQ=255)
                                                                   gcacagtGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCAc                                                                               <  2:3462153/150‑1 (MQ=255)
                                                                      cagcGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAAGTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTGTTTCACGTc                                                                            >  2:4214393/1‑150 (MQ=255)
                                                                       agcGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTcc                                                                           >  1:1919539/1‑150 (MQ=255)
                                                                        gcGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCa                                                                          <  2:2667486/150‑1 (MQ=255)
                                                                           ttCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGAt                                                                       <  2:4738590/150‑1 (MQ=255)
                                                                            tCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTcaa                                                                          >  1:1552981/1‑144 (MQ=255)
                                                                              tGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCa                                                                    >  2:304430/1‑150 (MQ=255)
                                                                                       ggggCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACgt                                                           >  1:4331278/1‑150 (MQ=255)
                                                                                       ggggCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACgt                                                           >  1:5407811/1‑150 (MQ=255)
                                                                                         ggCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACgtgt                                                         >  2:5212556/1‑150 (MQ=255)
                                                                                           ccACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTc                                                       <  2:2061339/150‑1 (MQ=255)
                                                                                           ccACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTc                                                       <  2:4892902/150‑1 (MQ=255)
                                                                                           ccACTGCGCAGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTc                                                       <  1:4702267/150‑1 (MQ=255)
                                                                                              cTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGtt                                                    <  1:737559/150‑1 (MQ=255)
                                                                                                    ggCCCCGGCAATTCAGAAGGTATAGGAACCGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCCTACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCtt                                              <  2:3126890/150‑1 (MQ=255)
                                                                                                     gCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTg                                             <  1:1018584/150‑1 (MQ=255)
                                                                                                     gCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTg                                             <  1:9020/150‑1 (MQ=255)
                                                                                                       cccGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTa                                           <  1:2051785/150‑1 (MQ=255)
                                                                                                       cccGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTa                                           >  2:3133407/1‑150 (MQ=255)
                                                                                                         cGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGa                                         >  1:4536695/1‑150 (MQ=255)
                                                                                                          ggCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGaa                                        <  2:2785325/150‑1 (MQ=255)
                                                                                                           gCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAg                                       <  1:4697889/150‑1 (MQ=255)
                                                                                                                tCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCaa                                  <  1:1942147/150‑1 (MQ=255)
                                                                                                                  aGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACa                                <  2:4952687/150‑1 (MQ=255)
                                                                                                                  aGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGCTGACCTTGTa                                           >  2:67077/1‑139 (MQ=255)
                                                                                                                    aaGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCTGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAgg                              <  2:1133837/150‑1 (MQ=255)
                                                                                                                        tataGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGa                          >  1:1531278/1‑150 (MQ=255)
                                                                                                                         ataGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGAc                         >  2:4790337/1‑150 (MQ=255)
                                                                                                                          taGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGAcc                        <  2:3520111/150‑1 (MQ=255)
                                                                                                                           aGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCa                       >  1:773437/1‑150 (MQ=255)
                                                                                                                            ggAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTGGAAGCTCAACAGGTCGACCAc                      >  2:1601495/1‑150 (MQ=255)
                                                                                                                             gAAACGCCGGCCTGCACCGTGCGCGTTGC‑‑GCCGGTGTATACATTGCTGTTGAACGCGGCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGTTTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACg                     <  1:3580314/150‑1 (MQ=255)
                                                                                                                             gAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACg                     >  2:3013347/1‑150 (MQ=255)
                                                                                                                                  gCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTa                >  2:831191/1‑150 (MQ=255)
                                                                                                                                  gCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTa                >  2:3782894/1‑150 (MQ=255)
                                                                                                                                  gCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTa                <  2:2064643/150‑1 (MQ=255)
                                                                                                                                  gCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTa                <  2:1710525/150‑1 (MQ=255)
                                                                                                                                  gCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTa                <  2:1007222/150‑1 (MQ=255)
                                                                                                                                   ccGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAg               <  1:2384967/150‑1 (MQ=255)
                                                                                                                                    cccggcgCACCGCGCACGGTGCAGGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTa                >  1:3888159/7‑150 (MQ=255)
                                                                                                                                    cccggcgCACCGCGCACGGTGCAGGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTa                >  1:3228487/7‑150 (MQ=255)
                                                                                                                                    cGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAgc              <  2:4660702/150‑1 (MQ=255)
                                                                                                                                     ccggcgCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAgcg             >  1:5608110/6‑150 (MQ=255)
                                                                                                                                        cTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAGCGCTc          >  1:4928178/1‑150 (MQ=255)
                                                                                                                                          gCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAGCGCTCCa        <  2:2180022/150‑1 (MQ=255)
                                                                                                                                          gCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAGCGCTCCa        <  1:5054743/150‑1 (MQ=255)
                                                                                                                                            aCCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAGCGCTCCATc      <  2:4512426/150‑1 (MQ=255)
                                                                                                                                            aCCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAGCGCTCCATc      <  1:2397518/150‑1 (MQ=255)
                                                                                                                                             ccGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAGCGCTCCATCg     >  1:3383455/1‑150 (MQ=255)
                                                                                                                                               gTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAGCGCTCCATCgtg   >  1:2583014/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                       
ATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAGCGCTCCATCGTGT  >  NC_002947/3578095‑3578388

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