Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,397,108 T→A 7.2% V17D (GTC→GAC) ‡ tolA → colicin S4/filamentous phage transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,397,1080TA 7.2% 439.2 / 16.1 135V17D (GTC→GAC) ‡tolAcolicin S4/filamentous phage transport protein
Reads supporting (aligned to +/- strand):  ref base T (69/56);  new base A (5/5);  total (74/61)
Fisher's exact test for biased strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 137 total aligned reads displayed.

GCGACAAGGCTGTCGACTATGGCGCGGTCATGGGTGCCATGGGCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGCTCAAGTCCAAGAGCCAGGCAACCACCCAGACCAATCAGA  >  NC_002947/1396962‑1397257
                                                                                                                                                  |                                                                                                                                                     
gCGACAAGGCTGTCGACTATGGCGCGGTCATGGGTGCCATGGGCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTg                                                                                                                                                    <  2:2526049/150‑1 (MQ=255)
 cGACAAGGCTGTCGACTATGGCGCGGTCATGGGTGCCATGGGCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTgg                                                                                                                                                   <  2:2493471/150‑1 (MQ=255)
 cGACAAGGCTGTCGACTATGGCGCGGTCATGGGTGCCATGGGCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTgg                                                                                                                                                   <  2:123432/150‑1 (MQ=255)
 cGACAAGGCTGTCGACTATGGCGCGGTCATGGGTGCCATGGGCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTgg                                                                                                                                                   <  2:3913210/150‑1 (MQ=255)
  gACAAGGCTGTCGACTATGGCGCGGTCATGGGTGCCATGGGCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTggg                                                                                                                                                  >  2:146091/1‑150 (MQ=255)
  gACAAGGCTGTCGACTATGGCGCGGTCATGGGTGCCATGGGCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTggg                                                                                                                                                  <  1:3335909/150‑1 (MQ=255)
  gACAAGGCTGTCGACTATGGCGCGGTCATGGGTGCCATGGGCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTggg                                                                                                                                                  >  2:4473009/1‑150 (MQ=255)
     aaGGCTGTCGACTATGGCGCGGTCATGGGTGCCATGGGCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGTGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGACACTGGGCCa                                                                                                                                               <  2:1962592/150‑1 (MQ=255)
     aaGGCTGTCGACTATGGCGCGGTCATGGGTGCCATGGGCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCa                                                                                                                                               >  1:4744552/1‑150 (MQ=255)
      aGGCTGTCGACTATGGCGCGGTCATGGGTGCCATGGGCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCAt                                                                                                                                              <  2:2350808/150‑1 (MQ=255)
       ggCTGTCGACTATGGCGCGGTCATGGGTGCCATGGGCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGACACTGGGCCaga                                                                                                                                             <  1:4240044/150‑3 (MQ=255)
        gCTGTCGACTATGGCGCGGTCATGGGTGCCATGGGCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGACACTGGGCCagaa                                                                                                                                            <  2:1557627/150‑4 (MQ=255)
        gCTGTCGACTATGGCGCGGTCATGGGTGCCATGGGCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGACACTGGGCCagaa                                                                                                                                            <  2:4950074/150‑4 (MQ=255)
        gCTGTCGACTATGGCGCGGTCATGGGTGCCATGGGCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGACACTGGGCCagaa                                                                                                                                            <  1:3079877/150‑4 (MQ=255)
          tGTCGACTATGGCGCGGTCATGGGTGCCATGGGCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGc                                                                                                                                          >  1:1970777/1‑150 (MQ=255)
          tGTCGACTATGGCGCGGTCATGGGTGCCAGGGGCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGc                                                                                                                                          <  2:305171/150‑1 (MQ=255)
                  aTGGCGCGGTCATGGGTGCCATGGGCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCAtgt                                                                                                                                  <  1:350326/150‑1 (MQ=255)
                  aTGGCGCGGTCATGGGTGCCATGGGCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCAtgt                                                                                                                                  <  1:350292/150‑1 (MQ=255)
                   tGGCGCGGTCATGGGTGCCATGGGCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCAtgtg                                                                                                                                 <  1:5107401/150‑1 (MQ=255)
                    ggCGCGGTCATGGGTGCCATGGGCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGc                                                                                                                                >  1:1080783/1‑150 (MQ=255)
                       gcgGTCATGGGTGCCATGGGCGGGTTGCAGAAGGCCGTTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTgg                                                                                                                             >  1:341794/1‑150 (MQ=255)
                        cgGTCATGGGTGCCATGGGCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGgt                                                                                                                            <  2:4547057/150‑1 (MQ=255)
                          gTCATGGGTGCCATGGGCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGgtgt                                                                                                                          >  1:3773569/1‑150 (MQ=255)
                           tCATGGGTGCCATGGGCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGtt                                                                                                                         >  1:4386981/1‑150 (MQ=255)
                            cATGGGTGCCATGGGCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTc                                                                                                                        <  1:1941742/150‑1 (MQ=255)
                                ggTGCCATGGGCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGccc                                                                                                                    <  2:2416644/150‑1 (MQ=255)
                                    gcATGGGCGGGTTGGAGAAGGCCGGTGTCGGTATCGTTGGCCTGCGTACCGAGGCGCCCTGCTGCAACCGCGAGAGCCTTCCGCCTCGGAACGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTGGCCctgct                                                                                                                <  2:366440/149‑1 (MQ=255)
                                     cATGGGCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCctgctg                                                                                                               >  1:2648143/1‑150 (MQ=255)
                                      aTGGGTGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGcc                                                                                                                                         <  1:1711153/123‑1 (MQ=255)
                                      aTGGGTGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGcc                                                                                                                                         >  2:1711153/1‑123 (MQ=255)
                                      aTGGGCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACCTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGt                                                                                                              >  2:4187663/1‑150 (MQ=255)
                                        gggCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTc                                                                                                            >  2:2234310/1‑150 (MQ=255)
                                        gggCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTc                                                                                                            <  1:2094175/150‑1 (MQ=255)
                                        gggCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTc                                                                                                            >  2:2234449/1‑150 (MQ=255)
                                         ggCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCg                                                                                                           >  2:3051069/1‑150 (MQ=255)
                                           cGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTc                                                                                                         >  1:3485666/1‑150 (MQ=255)
                                            gggTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCa                                                                                                        <  2:579860/150‑1 (MQ=255)
                                            gggTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCa                                                                                                        <  1:984157/150‑1 (MQ=255)
                                            gggTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCa                                                                                                        <  1:3478261/150‑1 (MQ=255)
                                              gTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGt                                                                                                      >  1:4289001/1‑150 (MQ=255)
                                              gTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGt                                                                                                      >  1:4289147/1‑150 (MQ=255)
                                                  cAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTg                                                                                                  <  2:1832191/150‑1 (MQ=255)
                                                   aGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGc                                                                                                 >  1:4123216/1‑150 (MQ=255)
                                                     aaGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCa                                                                                                        >  1:2275511/1‑141 (MQ=255)
                                                         ccGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGAc                                                                                           <  1:4463776/150‑1 (MQ=255)
                                                         ccGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGAc                                                                                           <  1:4463830/150‑1 (MQ=255)
                                                          cGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACg                                                                                          >  2:404790/1‑150 (MQ=255)
                                                               tCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGACCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGa                                                                                     <  2:2102601/150‑1 (MQ=255)
                                                               tCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGACCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGa                                                                                     <  2:2102604/150‑1 (MQ=255)
                                                                 ggTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGc                                                                                   >  2:5023336/1‑150 (MQ=255)
                                                                 ggTAACGGTGGCCTGATTACCGAGGCGCCCTGATGCACCAGCGACAGCCATCCGCCGCGGAAAGCTACTTCTGGCCCAGTCTCAGGGCCAACGCCCAGCAAGAGCTGGTCTACGCCCTCCTGTTCGTCAGGATAGCCATGACGCCGGAGc                                                                                   >  2:5023443/1‑150 (MQ=255)
                                                                         tGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCt                                                                           <  1:2116778/150‑1 (MQ=255)
                                                                            ccTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTcc                                                                        >  1:1233677/1‑150 (MQ=255)
                                                                             cTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCa                                                                       >  1:1147984/1‑150 (MQ=255)
                                                                                 ttACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGcc                                                                   >  2:4357818/1‑150 (MQ=255)
                                                                                   aCCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGa                                                                 <  1:1405835/150‑1 (MQ=255)
                                                                                   aCCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGa                                                                 <  1:4184841/150‑1 (MQ=255)
                                                                                    ccGAGGCGCCCTGATGCACCAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCctg                                                                                                                  <  2:4919313/100‑1 (MQ=255)
                                                                                    ccGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCctg                                                                                                                  >  1:4919313/1‑100 (MQ=255)
                                                                                       aGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGt                                                             >  2:2076043/1‑150 (MQ=255)
                                                                                         gcgcCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTTCTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTTGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTc                                                           >  2:191243/1‑150 (MQ=255)
                                                                                           gcCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAg                                                         >  1:626616/1‑150 (MQ=255)
                                                                                           gcCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAg                                                         >  1:2273072/1‑150 (MQ=255)
                                                                                            cccTGATGCAACAGCGAGAGTCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAgg                                                        <  1:2347590/150‑1 (MQ=255)
                                                                                                    cAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTg                                                >  1:2592302/1‑150 (MQ=255)
                                                                                                     aaCAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGACCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGt                                               >  2:1151204/1‑150 (MQ=255)
                                                                                                        aGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTAcc                                            >  2:3239817/1‑150 (MQ=255)
                                                                                                         gCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCa                                           <  2:1368669/150‑1 (MQ=255)
                                                                                                          cGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAg                                          <  2:116552/150‑1 (MQ=255)
                                                                                                           gagaGCCATCCGCCTCGGAAAGCTACTTCTGGCCCATTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGACCTGCTGTTCGTCAGTTTTGCCAAGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGATCAGGCTACCATGTACCAGc                                         >  2:505367/1‑150 (MQ=255)
                                                                                                           gagaGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGc                                         <  1:1387453/150‑1 (MQ=255)
                                                                                                            agagCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGCt                                        >  2:4499432/1‑150 (MQ=255)
                                                                                                            agagCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGCt                                        >  2:4499576/1‑150 (MQ=255)
                                                                                                                ccATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGCTCAAg                                    >  1:364181/1‑150 (MQ=255)
                                                                                                                ccATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGCTCAAg                                    >  1:3284894/1‑150 (MQ=255)
                                                                                                                 cATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGCTCAAGt                                   >  2:5072003/1‑150 (MQ=255)
                                                                                                                 cATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGCTCAAGt                                   >  1:2526361/1‑150 (MQ=255)
                                                                                                                   tCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGCTCAAGTcc                                 <  1:805658/150‑1 (MQ=255)
                                                                                                                   tCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGCTCAAGTcc                                 <  1:4473009/150‑1 (MQ=255)
                                                                                                                   tCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCAGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGCTCAAGTcc                                 >  1:1208897/1‑150 (MQ=255)
                                                                                                                       ccTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGCTCAAGTCCAaga                             <  1:3880888/150‑1 (MQ=255)
                                                                                                                         tCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGCTCAAGTCCAAGAGc                           <  2:2681822/150‑1 (MQ=255)
                                                                                                                         tCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGCTCAAGTCCAAGAGc                           >  1:2873797/1‑150 (MQ=255)
                                                                                                                            gAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGCTCAAGTCCAAGAGCCAg                        <  1:4770574/150‑1 (MQ=255)
                                                                                                                                 cTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGCTCAAGTCCAAGAGCCAGGCAac                   <  2:1725461/150‑1 (MQ=255)
                                                                                                                                 aggccgatGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGCTCAAGTCCAAGAGCCAGGCAac                   >  2:304050/8‑150 (MQ=255)
                                                                                                                                     cgatGGCCCAGACACTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGCTCAAGTCCAAGAGCCAGGCAaccacc               >  1:951124/4‑150 (MQ=255)
                                                                                                                                      tCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGCTCAAGTCCAAGAGCCAGGCAACCAccc              <  1:3003804/150‑1 (MQ=255)
                                                                                                                                      tCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGCTCAAGTCCAAGAGCCAGGCAACCAccc              <  1:432313/150‑1 (MQ=255)
                                                                                                                                      gatGGCCCAGACACTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGCTCAAGTCCAAGAGCCAGGCAACCAccc              >  1:1172821/3‑150 (MQ=255)
                                                                                                                                        tGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGCTCAAGTCCAAGAGCCAGGCAACCACCCAg            >  1:1437805/1‑150 (MQ=255)
                                                                                                                                         ggCCCAGACACTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGCTCAAGTCCAAGAGCCAGGCAACCACCCAGa           >  1:4190643/1‑150 (MQ=255)
                                                                                                                                          gCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGCTCAAGTCCAAGAGCCAGGCGACCACCCAGAc          >  1:1880371/1‑150 (MQ=255)
                                                                                                                                          gCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGCTCAAGTCCAAGAGCCAGGCAACCACCCAGAc          <  1:2419570/150‑1 (MQ=255)
                                                                                                                                           cccAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGCTCAAGTCCAAGAGCCAGGCAACCACCCAGAcc         >  1:3927613/1‑150 (MQ=255)
                                                                                                                                            ccAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGCTCAAGTCCAAGAGCCAGGCAACCACCCAGACCa        >  1:4085076/1‑150 (MQ=255)
                                                                                                                                              aGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGCTCAAGTCCAAGAGCCAGGCAACCACCCAGACCAAt      <  2:2648143/150‑1 (MQ=255)
                                                                                                                                              aGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGCTCAAGTCCAAGAGCCAGGCAACCACCCAGACCAAt      <  1:4187663/150‑1 (MQ=255)
                                                                                                                                                  tCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGCTCAAGTCCAAGAGCCAGGCAACCACCCAGACCAATCaga  >  2:4449224/1‑150 (MQ=255)
                                                                                                                                                  tCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGCTCAAGTCCAAGAGCCAGGCAACCACCCAGACCAATCaga  >  1:2647998/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCGACAAGGCTGTCGACTATGGCGCGGTCATGGGTGCCATGGGCGGGTTGCAGAAGGCCGGTGTCGGTAACGTTGGCCTGATTACCGAGGCGCCCTGATGCAACAGCGAGAGCCATCCGCCTCGGAAAGCTACTTCTGGCCCAGTGTCTGGGCCATCGGCCTGCATGTGCTGGTGTTCGCCCTGCTGTTCGTCAGTTTTGCCATGACGCCGGAGCTGCCGCCTTCCAAGCCGATCGTTCAGGCTACCCTGTACCAGCTCAAGTCCAAGAGCCAGGCAACCACCCAGACCAATCAGA  >  NC_002947/1396962‑1397257

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: