Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,767,863 C→T 8.4% intergenic (+32/+37) fpvA → / ← PP_4218 TonB‑dependent outer membrane ferripyoverdine receptor FpvA/lipase/esterase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,767,8630CT 8.4% 263.6 / 12.1 95intergenic (+32/+37)fpvA/PP_4218TonB‑dependent outer membrane ferripyoverdine receptor FpvA/lipase/esterase family protein
Reads supporting (aligned to +/- strand):  ref base C (45/42);  new base T (5/3);  total (50/45)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

CAACGTATTCGACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGC  >  NC_002947/4767717‑4768011
                                                                                                                                                  |                                                                                                                                                     
caacGTATTCGACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAgg                                                                                                                                                    >  1:1465617/1‑150 (MQ=255)
       ttCGACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGtt                                                                                                                                             >  2:2262196/1‑150 (MQ=255)
       ttCGACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGtt                                                                                                                                             >  2:438224/1‑150 (MQ=255)
          gACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGttttt                                                                                                                                          >  1:413394/1‑150 (MQ=255)
          gACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGttttt                                                                                                                                          >  2:2663/1‑150 (MQ=255)
           aCAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCTGGGTGCGGtttttt                                                                                                                                         <  1:1828911/150‑1 (MQ=255)
           aCAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCTGGGTGCGGtt                                                                                                                                             <  2:499178/146‑1 (MQ=255)
           aCAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCTGGGTGCGGtt                                                                                                                                             >  1:499178/1‑146 (MQ=255)
           aCAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGtttttt                                                                                                                                         >  1:80922/1‑150 (MQ=255)
           aCAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGtttttt                                                                                                                                         <  2:1307400/150‑1 (MQ=255)
            cAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTc                                                                                                                                        >  2:295690/1‑150 (MQ=255)
             aagaagTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCa                                                                                                                                       >  2:1913701/1‑150 (MQ=255)
             aagaagTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCa                                                                                                                                       <  2:1236389/150‑1 (MQ=255)
             aagaagTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTGTTTTCa                                                                                                                                       >  1:447477/1‑150 (MQ=255)
              agaagTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCAt                                                                                                                                      <  2:1744182/150‑1 (MQ=255)
                aagTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTAAGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATgg                                                                                                                                    <  1:870651/150‑1 (MQ=255)
                   tactacCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATg                                                                                                                                 >  1:614575/1‑150 (MQ=255)
                        ccAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCa                                                                                                                                       <  2:213731/139‑1 (MQ=255)
                        ccAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCa                                                                                                                                       >  1:213731/1‑139 (MQ=255)
                              cATTGCCAGCAACACCTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCt                                                                                                                      <  1:2088683/150‑1 (MQ=255)
                              cATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCt                                                                                                                              <  1:1431831/142‑1 (MQ=255)
                              cATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCt                                                                                                                              >  2:1431831/1‑142 (MQ=255)
                                ttGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTcc                                                                                                                    <  1:2218616/150‑1 (MQ=255)
                                   ccAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCTCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCtt                                                                                                                 >  1:1098860/1‑150 (MQ=255)
                                    cAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTa                                                                                                                >  2:2396726/1‑150 (MQ=255)
                                    cAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTa                                                                                                                >  1:1458766/1‑150 (MQ=255)
                                      gCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACa                                                                                                              <  2:614575/150‑1 (MQ=255)
                                           aacTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtgc                                                                                                         <  2:11037/150‑1 (MQ=255)
                                              tCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtgctgc                                                                                                      <  2:756091/150‑1 (MQ=255)
                                               ccTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtgctgct                                                                                                     <  1:1633690/150‑1 (MQ=255)
                                                cTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtgctgctg                                                                                                    <  1:438224/150‑1 (MQ=255)
                                                     gCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGGTCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTc                                                                                                                     <  2:2197357/128‑1 (MQ=255)
                                                     gCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGGTCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTc                                                                                                                     >  1:2197357/1‑128 (MQ=255)
                                                         gTACGACGTGTATGGTGACCCATGCGACGTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTTCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTc                                                                                          >  1:328909/1‑150 (MQ=255)
                                                         gTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCTCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTc                                                                                          <  2:1098860/150‑1 (MQ=255)
                                                           acgatgTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGAGTGTTTGTTTTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGa                                                                                        <  2:328909/150‑1 (MQ=255)
                                                                  gtATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGaaa                                                                                                                          >  2:1577552/1‑110 (MQ=255)
                                                                  gtATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGaaa                                                                                                                          <  1:1577552/110‑1 (MQ=255)
                                                                   tATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCt                                                                                                                      <  2:966315/113‑1 (MQ=255)
                                                                   tATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCt                                                                                                                      >  1:966315/1‑113 (MQ=255)
                                                                   tATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATc                                                                                >  2:2182907/1‑150 (MQ=255)
                                                                       gTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTccc                                                                                                                   <  1:1214369/112‑1 (MQ=255)
                                                                       gTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTccc                                                                                                                   >  2:1214369/1‑112 (MQ=255)
                                                                       gTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtgctgct                                                                                                     <  1:1274500/126‑1 (MQ=255)
                                                                       gTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtgctgct                                                                                                     >  2:1274500/1‑126 (MQ=255)
                                                                        tGACCCACTCAACTGCAGCATCAGCGCAAAGTACAGCTTCTGATCGCCTGACGTTGAACGCACAAAACCGCACCCAGGTGCGGTTTTTTCATGGATTCCTGGAATGCTTCCTTGCAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAg                                                                           <  2:413394/150‑1 (MQ=255)
                                                                         gACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAgg                                                                          <  2:1997850/150‑1 (MQ=255)
                                                                          aCCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGaaa                                                                                                                          >  2:1971867/1‑102 (MQ=255)
                                                                          aCCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGaaa                                                                                                                          <  1:1971867/102‑1 (MQ=255)
                                                                          aCCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTcc                                                                                                                    <  1:1313433/108‑1 (MQ=255)
                                                                          aCCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTcc                                                                                                                    >  2:1313433/1‑108 (MQ=255)
                                                                           cccACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACACAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATTCTGCTGCAGTGCCTCGACCACCATCGCCAGGcc                                                                        >  1:275634/1‑150 (MQ=255)
                                                                            ccACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCt                                                                       >  1:1509819/1‑150 (MQ=255)
                                                                             cACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCtt                                                                      >  2:1491504/1‑150 (MQ=255)
                                                                              aCGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCtcc                                                                     >  2:1979607/1‑148 (MQ=255)
                                                                                  aaCTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCg                                                                 <  1:2182907/150‑1 (MQ=255)
                                                                                  aaCTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCg                                                                 <  2:1465617/150‑1 (MQ=255)
                                                                                  aaCTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCg                                                                 <  2:169459/150‑1 (MQ=255)
                                                                                  aaCTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGATTCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCg                                                                 >  1:1109839/1‑150 (MQ=255)
                                                                                    cTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGc                                                               <  2:1621800/150‑1 (MQ=255)
                                                                                     ttCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtg                                                                                                          >  2:1950805/1‑107 (MQ=255)
                                                                                     ttCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGCTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtg                                                                                                          <  1:1950805/107‑1 (MQ=255)
                                                                                      tCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAg                                                                                                                                                     <  2:443929/63‑1 (MQ=255)
                                                                                      tCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAg                                                                                                                                                     >  1:443929/1‑63 (MQ=255)
                                                                                        aCCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCa                                                                                                                                       <  2:237035/75‑1 (MQ=255)
                                                                                        aCCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCa                                                                                                                                       >  1:237035/1‑75 (MQ=255)
                                                                                        aCCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGcctcc                                                           <  2:1458766/150‑1 (MQ=255)
                                                                                         ccATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGcctcct                                                          <  1:1626134/150‑1 (MQ=255)
                                                                                              aCCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCgatg                                                     >  1:2380865/1‑150 (MQ=255)
                                                                                                 gCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATc                                                  <  2:2380865/150‑1 (MQ=255)
                                                                                                          aGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATc                                                  <  1:1040134/141‑1 (MQ=255)
                                                                                                          aGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATc                                                  >  2:1040134/1‑141 (MQ=255)
                                                                                                          aGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATa                                         >  2:1410032/1‑150 (MQ=255)
                                                                                                           gCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATAc                                        <  1:2663/150‑1 (MQ=255)
                                                                                                            cTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACg                                       >  2:96082/1‑150 (MQ=255)
                                                                                                                     gCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTc                              <  1:2262196/150‑1 (MQ=255)
                                                                                                                        tGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCgat                                                      <  2:2398235/123‑1 (MQ=255)
                                                                                                                        tGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCgat                                                      >  1:2398235/1‑123 (MQ=255)
                                                                                                                        tGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATg                           >  1:2171058/1‑150 (MQ=255)
                                                                                                                         gACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACCACCTTTCCACTCTGGCTGCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATgg                          <  2:695046/150‑1 (MQ=255)
                                                                                                                                  cGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTccc                 >  1:1182861/1‑150 (MQ=255)
                                                                                                                                     aaaaaaCCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATgg                          >  2:2141988/1‑138 (MQ=255)
                                                                                                                                     aaaaaaCCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATgg                          <  1:2141988/138‑1 (MQ=255)
                                                                                                                                      aaaaaCCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGGCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGa             >  1:451886/1‑150 (MQ=255)
                                                                                                                                       aaaaCCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATgg                          <  2:2359111/136‑1 (MQ=255)
                                                                                                                                       aaaaCCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATgg                          >  1:2359111/1‑136 (MQ=255)
                                                                                                                                       aaaaCCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCg                >  1:459710/1‑146 (MQ=255)
                                                                                                                                       aaaaCCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCg                <  2:459710/146‑1 (MQ=255)
                                                                                                                                        aaaCCGCACCTGGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATc           >  2:753084/1‑150 (MQ=255)
                                                                                                                                        aaaCCGCACCTGGGTGCGGTTTTTTCANGGATGCCTGAAANGCTCCCNTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATc           >  1:1907170/1‑150 (MQ=255)
                                                                                                                                         aaCCGCACCTGGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCa          >  2:320350/1‑150 (MQ=255)
                                                                                                                                         aaCCGCACCTGGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCa          >  2:2196264/1‑150 (MQ=255)
                                                                                                                                         aaCCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATcc          >  1:2369906/1‑149 (MQ=255)
                                                                                                                                           ccGCACCTGGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCAcc        <  1:753084/150‑1 (MQ=255)
                                                                                                                                           ccGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCAcc        >  2:1995914/1‑150 (MQ=255)
                                                                                                                                                 ccAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTccc                 <  2:2249893/135‑1 (MQ=255)
                                                                                                                                                 ccAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTccc                 >  1:2249893/1‑135 (MQ=255)
                                                                                                                                                 ccAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGc  >  2:2384057/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CAACGTATTCGACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGC  >  NC_002947/4767717‑4768011

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