Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,742,462 T→C 11.1% intergenic (+49/+22) hutF → / ← PP_5037 formiminoglutamate deiminase/lipocalin family lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,742,4620TC11.1% 245.6 / 16.9 90intergenic (+49/+22)hutF/PP_5037formiminoglutamate deiminase/lipocalin family lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (28/52);  new base C (7/3);  total (35/55)
Fisher's exact test for biased strand distribution p-value = 4.28e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

GATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTGCGCGACA  >  NC_002947/5742314‑5742611
                                                                                                                                                     |                                                                                                                                                     
gATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGtt                                                                                                                                                      <  2:612533/150‑1 (MQ=255)
 aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTg                                                                                                                                                     <  1:78969/150‑1 (MQ=255)
 aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTg                                                                                                                                                     <  2:718463/150‑1 (MQ=255)
 aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTg                                                                                                                                                     <  2:718424/150‑1 (MQ=255)
  tGGTGAAGGGGCACTGGGTGGTGCGCCAGGGGCGGCATGCCCAGGACGATGACCGCGGGCGGGCGTGTGCCACGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgc                                                                                                                                                    <  2:1154015/150‑1 (MQ=255)
    gTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaa                                                                                                                                                  <  2:2152282/150‑1 (MQ=255)
    gTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaa                                                                                                                                                  <  2:4626328/150‑1 (MQ=255)
    gTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaa                                                                                                                                                  <  1:855407/150‑1 (MQ=255)
     tGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAGGGCTgcaac                                                                                                                                                 <  1:4919332/150‑1 (MQ=38)
         tGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaacgcaa                                                                                                                                             <  1:1543578/150‑1 (MQ=255)
          gggCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaacgcaac                                                                                                                                            <  2:4542162/150‑1 (MQ=255)
            gCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTg                                                                                                                                                     <  2:2742283/139‑1 (MQ=255)
            gCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTg                                                                                                                                                     >  1:2742283/1‑139 (MQ=255)
            gCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGGTGCAGCCCCTAGCGGGCTCGCCCGCGCAGGGTTGCAACCCCAccc                                                                                                                                          >  2:3604584/1‑150 (MQ=255)
                  ggtggtGCGCCCGGGGCGGCAGGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCGGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCaaaaa                                                                                                                                    <  2:5407448/150‑1 (MQ=255)
                   gtggtgCGCCTGGGGCGGCTTGCCCAGGAAGATGAACGGGGGCGGGCTTTGGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCaaaaaa                                                                                                                                   <  2:3136458/150‑1 (MQ=255)
                       ggcgcggggggCGGCCTGCCCTGTAATTTTACGGGGGGGGGGGGTTTGCATCGGTGTGGCGCCTTTTTTGGGGCTGTGGGTCGGGGGGTTGGGTGGCCGCCCTGCGGGTGCACGCCCGCGATGGGTTCAACCGCAACCCCAAAAAACTCa                                                                                                                               <  2:1985021/143‑1 (MQ=255)
                       cgcgcCTGGGGGGGCTTGCCCTGTAAGATTTAAGGGGGGGGGGGTTTGCCATGGTGTTGCGCCTGCTTTTGGGCTGAGCGTACGGGGGTTGCGTTGCAGCCCTTCGCGTGAAAGCCCGCGATGGGTTGCAACGCAACCCCAAAAAACTCa                                                                                                                               <  1:4990194/149‑1 (MQ=255)
                        gcgcCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGGGGTGGAGCCCGTAGCCGGCCAGCCCCCCCACGGGTTGAAAGCAAACCCaaaaaaataaa                                                                                                                              >  1:4536115/1‑145 (MQ=255)
                            cAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAGGGCTg                                                                                                                                                     <  1:4677511/123‑1 (MQ=37)
                            cAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAGGGCTg                                                                                                                                                     >  2:4677511/1‑123 (MQ=37)
                              ggggCGTCAGGCCCAGTATGATGAACGGGGGCGGGGTTTGGCCATGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCAtt                                                                                                                        <  1:5444137/150‑1 (MQ=255)
                                         cccAGTAAGATGACCGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGTGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCTAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTg                                                                                                             <  2:817268/150‑1 (MQ=255)
                                         cccAGGAAGATGAAAGGGGGCGGGCTTTTGCCATGGTGTTGCGCCTGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTg                                                                                                             <  2:817193/150‑1 (MQ=255)
                                             ggAAGAGGACCGGGGGCGGGCGTTGGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCACGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGa                                                                                                         <  2:3725883/150‑1 (MQ=255)
                                                aGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCg                                                                                                      <  2:5321170/150‑1 (MQ=255)
                                                           ggcgggggTTTGCCATGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACGGGGGTTGCGTTGCAGCCCTTGGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCa                                                                                           <  2:800458/150‑1 (MQ=255)
                                                                 cgTTTGCCAAGGTGTTGCGCCTGCTTTTGGGCTGAGCGTACGGGGGTTGCGTTGCCGCCCTTCGCGTGCAAGCCCGCTATGGGTTGCACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCaa                                                                                     <  2:5176645/150‑1 (MQ=255)
                                                                        cacggTGTTGCGCCTGCTTTGGGGCTGAGCGTACCGGGGTGGCGTGGCAGCCCTTGGCGGGCACGCCCGCGAAGGGTTGCAACGCAACCCCAAAACACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTg                                                                              <  2:3427309/147‑1 (MQ=255)
                                                                        cacggTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCGTTCGCGTGCAAGCCCGCTATGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTg                                                                              <  2:1220632/147‑1 (MQ=255)
                                                                          tgtggtTGCGCCCGGTTTGGGGCTGCGCTTACGGGGGTTGGGTTGAAGCCGTTGGGGTGCTAGCCGGGGATGGGTTCCACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGgt                                                                            <  1:304407/146‑1 (MQ=255)
                                                                          tggTGTTGCGCATGTTTTTGGGCTGAGGGTCCGGGGTTTGGGTTGCAGCCCTTCGCGTGCATGCCGGCTATGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGgt                                                                            <  1:1721169/149‑1 (MQ=255)
                                                                          aGGTGTTGGGCCAGCTTTTGGGGTGAGCGTAGCGGGGTTGCGTTGCAGCCCTTCGCGTTCAAGCCCGCGAAGGTTTGCAACGCAACCCCAAAAAACTCAGTGCTTTTTTACTATTTGATGATCGGACGCCCGCCATATCAAGCGGCTGgt                                                                            <  2:4939340/150‑1 (MQ=255)
                                                                          aGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTGGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGgt                                                                            <  2:4939247/150‑1 (MQ=255)
                                                                                   gcCTGCTTTGGGGCTGAGCGTACCGGGGTTGCGTTGCCGCCCTGCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAccc                                                                   <  2:2498380/150‑1 (MQ=255)
                                                                                   gcCCGCTTTGGGGCTGCGCGTCCCGGGGTGGCGTTGCAGCCCTTCGCGTGCACGCCCGCGACGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAccc                                                                   <  2:2318943/150‑1 (MQ=255)
                                                                                     ctgCTTTTGGGCTGAGCGTACGGGGGTTGCGTTGCAGCCCTTCGCTTGCAAGCCCGCGATGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCct                                                                 <  1:5394125/148‑1 (MQ=255)
                                                                                       cattttGCGCTGAGCGTACCGGGGTTCCGTTCCACCCAATCCCGTGAAAGCCCGCGAAGGCTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctgc                                                               <  2:2563550/148‑1 (MQ=255)
                                                                                                   aGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGAAAGCCCGCGACGGGTTGCACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTg                                                   <  1:4451088/150‑1 (MQ=255)
                                                                                                    gCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCACGCCCGCGACGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTgg                                                  <  2:3303853/150‑1 (MQ=255)
                                                                                                         ccGGGGTTGCGTTGCCGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagc                                             <  2:2619313/150‑1 (MQ=255)
                                                                                                          cGGGGTTGGGTTGCAGCCCTTCGCGTGCACGCCCGCTATGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagca                                            <  1:3113118/150‑1 (MQ=255)
                                                                                                          cGGGGTTGGGGTGCAGCCCTTCGCGTGCAACCCCGCGATGGTTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagca                                            <  1:3112805/150‑1 (MQ=255)
                                                                                                                 gcgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtctc                                     <  2:1278326/150‑1 (MQ=255)
                                                                                                                 gcggtgcAGCCCTTCGCGTGAAAGCCCGCGAAGGGTAGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtctc                                     <  2:4257399/150‑1 (MQ=255)
                                                                                                                  cgttgcAGCCCTTGGCGTGCACGCCCGCGATGGTTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCa                                    <  1:842083/150‑1 (MQ=255)
                                                                                                                   gttgcAGCCCTTCGTGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCAc                                   >  1:4383293/1‑150 (MQ=255)
                                                                                                                   gttgcAGCCCTGGGGGTGCAAGCCGGCGATGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCAc                                   <  1:1660698/150‑1 (MQ=255)
                                                                                                                   gttgcAACCCTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAACAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtctctc                                   >  2:3685615/1‑148 (MQ=37)
                                                                                                                      gcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCt                                >  2:2835951/1‑150 (MQ=255)
                                                                                                                       cAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTc                                                                         >  1:675260/1‑108 (MQ=255)
                                                                                                                       cAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTc                                                                         <  2:675260/108‑1 (MQ=255)
                                                                                                                        aacccTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGcc                                                           <  1:1374005/119‑1 (MQ=38)
                                                                                                                        aacccTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGcc                                                           >  2:1374005/3‑121 (MQ=38)
                                                                                                                        aacccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctgc                                                               <  2:148558/115‑1 (MQ=37)
                                                                                                                        aacccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctgc                                                               >  1:148558/3‑117 (MQ=37)
                                                                                                                        aacccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc                              >  2:1088457/3‑150 (MQ=38)
                                                                                                                        aacccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc                              >  1:4849355/3‑150 (MQ=38)
                                                                                                                        aacccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc                              >  2:3176921/3‑150 (MQ=38)
                                                                                                                        aGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc                              >  1:395983/1‑150 (MQ=255)
                                                                                                                         acccTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGgt                                                                            >  2:1392007/2‑103 (MQ=37)
                                                                                                                         acccTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGgt                                                                            <  1:1392007/102‑1 (MQ=37)
                                                                                                                         acccTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTccc                             >  1:4571210/2‑150 (MQ=255)
                                                                                                                          cccTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCa                            >  1:1731629/1‑150 (MQ=255)
                                                                                                                            cTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAATCTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCTTACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg                          >  2:3071646/1‑150 (MQ=39)
                                                                                                                             ttCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc                         >  2:1520641/1‑150 (MQ=255)
                                                                                                                             ttCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc                         >  1:2477800/1‑150 (MQ=255)
                                                                                                                              tCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGACCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGcc                        >  1:2398883/1‑150 (MQ=39)
                                                                                                                                  gTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCa                    >  1:2587536/1‑150 (MQ=255)
                                                                                                                                  gTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGACATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCa                    >  2:1721965/1‑150 (MQ=255)
                                                                                                                                  gTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATGGTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCa                    >  2:3300097/1‑150 (MQ=255)
                                                                                                                                   tGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCt                                                    >  2:3671337/1‑117 (MQ=255)
                                                                                                                                   tGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCt                                                    <  1:3671337/117‑1 (MQ=255)
                                                                                                                                   tGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAg                   >  1:3699730/1‑150 (MQ=255)
                                                                                                                                     cAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAgtg                 >  1:1707209/1‑150 (MQ=255)
                                                                                                                                     cAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAgtg                 >  1:1707291/1‑150 (MQ=255)
                                                                                                                                      aaGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAgtgt                >  2:2019809/1‑150 (MQ=255)
                                                                                                                                      aaGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAgtgt                >  2:1637593/1‑150 (MQ=255)
                                                                                                                                      aaGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAgtgt                >  2:2815807/1‑150 (MQ=255)
                                                                                                                                       aGCCCGCTATGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTc               <  1:4819840/150‑1 (MQ=255)
                                                                                                                                        gCCCGCGATGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCg              <  1:3604584/150‑1 (MQ=255)
                                                                                                                                        gCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCACGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCg              <  2:4505010/150‑1 (MQ=255)
                                                                                                                                         gcacgcgAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCg              <  2:2477800/147‑1 (MQ=255)
                                                                                                                                           ccGCGACGGGTTGCACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgt            <  2:887699/150‑1 (MQ=255)
                                                                                                                                           ccGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGa                                                 <  1:3999941/113‑1 (MQ=255)
                                                                                                                                           ccGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGa                                                 >  2:3999941/1‑113 (MQ=255)
                                                                                                                                           ccGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgt            >  1:885457/1‑150 (MQ=255)
                                                                                                                                           ccGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgt            <  2:887480/150‑1 (MQ=255)
                                                                                                                                             gcgAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgtgt          <  2:1106024/150‑1 (MQ=255)
                                                                                                                                             gcgAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgtgt          >  2:4472419/1‑150 (MQ=255)
                                                                                                                                                 ggggTTGCGTTGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTgcgc      >  1:2248073/2‑150 (MQ=255)
                                                                                                                                                     ttGCACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTGCGCGACa  <  2:2981732/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTGCGCGACA  >  NC_002947/5742314‑5742611

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