Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,578,242 G→A 5.1% R699C (CGC→TGC)  PP_3155 ← ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,578,2420GA 5.1% 826.3 / 20.4 256R699C (CGC→TGC) PP_3155ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base G (117/125);  new base A (2/11);  total (120/136)
Fisher's exact test for biased strand distribution p-value = 2.28e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

Only 100 of 260 total aligned reads displayed.

ATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAGCGCTCCATCGTGTAG  >  NC_002947/3578095‑3578390
                                                                                                                                                   |                                                                                                                                                      
aTGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGcgcg                                                                                                                                                      <  1:3592129/150‑1 (MQ=255)
aTGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCGGGTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGACGGTATAGGAAACGCCGGCCTGCACCGTGcgcg                                                                                                                                                      >  1:1299927/1‑150 (MQ=255)
    actactCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTgc                                                                                                                                                  >  2:701281/1‑150 (MQ=255)
    actactCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTgc                                                                                                                                                  >  1:3372318/1‑150 (MQ=255)
    actactCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTgc                                                                                                                                                  >  1:2719733/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                                                                                                                                           >  2:2351117/1‑150 (MQ=255)
         tCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGCGCACGGTgcaggcc                                                                                                                                             <  1:1384386/150‑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)
          ccAACGCCGGGATAAAGCAGCAAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCggg                                                                                                                                          <  1:3839117/150‑1 (MQ=255)
           cAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCTgggg                                                                                                                                         >  1:1133837/1‑150 (MQ=255)
           cAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCgggg                                                                                                                                         >  2:5054743/1‑150 (MQ=255)
           cAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCgggg                                                                                                                                         >  1:4660702/1‑150 (MQ=255)
           cAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCgggg                                                                                                                                         >  1:1990075/1‑150 (MQ=255)
            aaCGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGt                                                                                                                                        <  1:957292/150‑1 (MQ=255)
                  gggATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACAt                                                                                                                                  >  1:1799554/1‑150 (MQ=255)
                     aTAAAGCAGCCAAACGACGGCATAGAGTCGTACTCTCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGc                                                                                                                               >  1:176208/1‑150 (MQ=255)
                        aaGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACCTTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGt                                                                                                                            <  2:3466164/150‑1 (MQ=255)
                        aaGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGt                                                                                                                            <  1:4074351/150‑1 (MQ=255)
                        aaGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGt                                                                                                                            >  1:1710525/1‑150 (MQ=255)
                         agcagcCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGtt                                                                                                                           >  1:525561/1‑150 (MQ=255)
                             gcCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAAc                                                                                                                       <  1:648191/150‑1 (MQ=255)
                                 aaCGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACgcgc                                                                                                                   <  2:4599512/150‑1 (MQ=255)
                                 aaCGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACgcgc                                                                                                                   <  2:273237/150‑1 (MQ=255)
                                         cATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCa                                                                                                           >  1:2868217/1‑150 (MQ=255)
                                                 tcagaacGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCtt                                                                                                   >  2:2684405/7‑150 (MQ=255)
                                                   tACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCttgt                                                                                                 >  1:4756307/1‑150 (MQ=255)
                                                   agaacGCTGTGCTGTAGAACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCttgt                                                                                                 >  1:245514/5‑150 (MQ=255)
                                                     cTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCttgttg                                                                                               >  2:5124681/1‑150 (MQ=255)
                                                     cTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCttgttg                                                                                               <  2:1528127/150‑1 (MQ=255)
                                                     cTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCttgttg                                                                                               >  1:3520111/1‑150 (MQ=255)
                                                        gCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAac                                                                                            >  1:2237538/1‑150 (MQ=255)
                                                           gTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACg                                                                                         >  2:2397518/1‑150 (MQ=255)
                                                            ttCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGt                                                                                        >  1:3507119/1‑150 (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)
                                                                   gcacagcGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCAc                                                                                 <  2:2986961/150‑1 (MQ=255)
                                                                      cagcGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAAGTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTGTTTCACGTc                                                                              >  2:4214393/1‑150 (MQ=255)
                                                                       agcGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTcc                                                                             >  1:1919539/1‑150 (MQ=255)
                                                                           ttCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGAt                                                                         <  2:4738493/150‑1 (MQ=255)
                                                                             cTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTc                                                                       <  1:208139/150‑1 (MQ=255)
                                                                              tGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCa                                                                      >  2:304430/1‑150 (MQ=255)
                                                                                     aCGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGctc                                                               <  2:3964287/150‑3 (MQ=255)
                                                                                       ggggCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACgt                                                             >  1:5407811/1‑150 (MQ=255)
                                                                                       ggggCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACgt                                                             >  2:1740713/1‑150 (MQ=255)
                                                                                         ggCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACgtgt                                                           >  2:5212556/1‑150 (MQ=255)
                                                                                           ccACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTc                                                         <  2:4892902/150‑1 (MQ=255)
                                                                                           ccACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTc                                                         <  2:2061339/150‑1 (MQ=255)
                                                                                           ccACTGCGCAGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTc                                                         <  1:4702231/150‑1 (MQ=255)
                                                                                            cACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCg                                                        >  1:3829559/1‑150 (MQ=255)
                                                                                              cTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGtt                                                      <  1:737559/150‑1 (MQ=255)
                                                                                               tGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTg                                                     >  1:3975331/1‑150 (MQ=255)
                                                                                                    ggCCCCGGCAATTCAGAAGGTATAGGACACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCtt                                                <  2:3127098/150‑1 (MQ=255)
                                                                                                     gCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTg                                               <  1:9020/150‑1 (MQ=255)
                                                                                                     gCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTg                                               <  1:1018584/150‑1 (MQ=255)
                                                                                                       cccGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTa                                             >  2:3133407/1‑150 (MQ=255)
                                                                                                       cccGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCACCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTa                                             <  2:4322994/150‑1 (MQ=255)
                                                                                                         cGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGa                                           >  1:4536695/1‑150 (MQ=255)
                                                                                                         cGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGa                                           >  1:3609019/1‑150 (MQ=255)
                                                                                                          ggCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGaa                                          <  2:2785325/150‑1 (MQ=255)
                                                                                                                 cAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAAc                                   <  2:2778950/150‑1 (MQ=255)
                                                                                                                  aGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACa                                  <  2:4952687/150‑1 (MQ=255)
                                                                                                                  aGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGCTGACCTTGTa                                             >  2:67077/1‑139 (MQ=255)
                                                                                                                  aGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCCTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACa                                  <  2:267290/150‑1 (MQ=255)
                                                                                                                    aaGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCTGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAgg                                <  2:1133837/150‑1 (MQ=255)
                                                                                                                     aGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGt                               <  1:1972380/150‑1 (MQ=255)
                                                                                                                        tataGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGa                            >  1:1531278/1‑150 (MQ=255)
                                                                                                                         ataGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGAc                           >  2:4790337/1‑150 (MQ=255)
                                                                                                                          taGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGAcc                          >  1:5142244/1‑150 (MQ=255)
                                                                                                                           aGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCa                         >  1:773437/1‑150 (MQ=255)
                                                                                                                            ggAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTGGAAGCTCAACAGGTCGACCAc                        >  2:1601495/1‑150 (MQ=255)
                                                                                                                            ggAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCAc                        >  2:3607085/1‑150 (MQ=255)
                                                                                                                             gAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACg                       >  2:3013347/1‑150 (MQ=255)
                                                                                                                              aaaCGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGc                      >  2:3216543/1‑150 (MQ=255)
                                                                                                                               aaCGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCt                     <  1:4343305/150‑1 (MQ=255)
                                                                                                                                  gCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTa                  >  2:831191/1‑150 (MQ=255)
                                                                                                                                  gCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTa                  >  2:3782894/1‑150 (MQ=255)
                                                                                                                                  gCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTa                  <  1:4167280/150‑1 (MQ=255)
                                                                                                                                  gCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTa                  <  2:1007222/150‑1 (MQ=255)
                                                                                                                                  gCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTa                  <  2:1710525/150‑1 (MQ=255)
                                                                                                                                  gCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTa                  <  2:2064643/150‑1 (MQ=255)
                                                                                                                                   ccGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAg                 <  1:2384967/150‑1 (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)
                                                                                                                                       ccTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAGCGCt             <  2:3325339/150‑1 (MQ=255)
                                                                                                                                        cTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAGCGCTc            >  1:4928283/1‑150 (MQ=255)
                                                                                                                                        cTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAGCGCTc            >  1:1554814/1‑150 (MQ=255)
                                                                                                                                          gCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAGCGCTCCa          <  2:2180022/150‑1 (MQ=255)
                                                                                                                                            aCCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAGCGCTCCATc        <  2:4512426/150‑1 (MQ=255)
                                                                                                                                             ccGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAGCGCTCCATCg       >  1:3383455/1‑150 (MQ=255)
                                                                                                                                               gTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAGCGCTCCATCgtg     >  1:2583014/1‑150 (MQ=255)
                                                                                                                                                  cgcgGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAGCGCTCCATCGTGTAg  >  2:411293/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
ATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAGCGCTCCATCGTGTAG  >  NC_002947/3578095‑3578390

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