Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,229,253:1 +T 6.2% coding (1776/2391 nt) PP_4606 ← ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,229,2531.T 6.2% 606.7 / 33.5 144coding (1776/2391 nt)PP_4606ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base . (75/60);  new base T (4/5);  total (79/65)
Fisher's exact test for biased strand distribution p-value = 7.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 157 total aligned reads displayed.

CTGTAGACCGTCTCGCTTGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGG  >  NC_002947/5229104‑5229401
                                                                                                                                                      |                                                                                                                                                    
 tGTAGACCGTCTCGCTTGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCg                                                                                                                                                     >  1:3335711/1‑150 (MQ=255)
  gTAGACCGTCTCGCTTGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCga                                                                                                                                                    <  2:17823/150‑1 (MQ=255)
   tAGACCGTCTCGCTTGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCgag                                                                                                                                                   <  2:2598552/150‑1 (MQ=255)
   tAGACCGTCTCGCTTGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCgag                                                                                                                                                   >  2:1666133/1‑150 (MQ=255)
   tAGACCGTCTCGCTTGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCgag                                                                                                                                                   <  1:5002143/150‑1 (MQ=255)
       ccGTCTCGCTTGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAAGGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCGGTGGCACCTgaggagg                                                                                                                                               <  2:2873946/150‑2 (MQ=255)
          tctcGCTTGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGa                                                                                                                                            <  1:4549454/150‑1 (MQ=255)
          tctcGCTTGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTCGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGa                                                                                                                                            >  2:2490671/1‑150 (MQ=255)
              gCTTGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCAc                                                                                                                                        >  2:4908038/1‑150 (MQ=255)
              gCTTGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCAc                                                                                                                                        >  2:3757182/1‑150 (MQ=255)
                ttGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGa                                                                                                                                      >  1:1215283/1‑150 (MQ=255)
                ttGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGa                                                                                                                                      >  1:1215277/1‑150 (MQ=255)
                  gtgACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGAGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGAcc                                                                                                                                    >  2:3803523/1‑149 (MQ=255)
                        gCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCg                                                                                                                              <  1:3388012/150‑1 (MQ=255)
                        gCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCg                                                                                                                              <  1:3937371/150‑1 (MQ=255)
                          gTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGt                                                                                                                            <  1:3520530/150‑1 (MQ=255)
                                   cGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGgcgc                                                                                                                   <  1:5026246/150‑1 (MQ=255)
                                   cGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGgcgc                                                                                                                   <  1:5027508/150‑1 (MQ=255)
                                    gagaCCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGcgcg                                                                                                                  >  1:411462/1‑150 (MQ=255)
                                       aCCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATg                                                                                                               >  2:2522683/1‑150 (MQ=255)
                                         cAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGt                                                                                                             >  1:1083659/1‑150 (MQ=255)
                                          aGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTc                                                                                                            <  1:3926868/150‑1 (MQ=255)
                                           gTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCg                                                                                                           >  2:2893845/1‑150 (MQ=255)
                                              cATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGt                                                                                                        >  2:4484570/1‑150 (MQ=255)
                                               aTTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGtt                                                                                                       <  1:3089184/150‑1 (MQ=255)
                                               aTTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGtt                                                                                                       >  2:3264637/1‑150 (MQ=255)
                                               aTTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGtt                                                                                                       >  1:1159375/1‑150 (MQ=255)
                                                ttGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTg                                                                                                      >  2:4799847/1‑150 (MQ=255)
                                                  gcgcTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTggg                                                                                                    <  1:3812397/150‑1 (MQ=255)
                                                   cgcTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGc                                                                                                   <  2:2813047/150‑1 (MQ=255)
                                                   cgcTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGc                                                                                                   <  2:838989/150‑1 (MQ=255)
                                                   cgcTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGc                                                                                                   <  2:839136/150‑1 (MQ=255)
                                                   cgcTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGc                                                                                                   >  2:1039868/1‑150 (MQ=255)
                                                    gcTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCt                                                                                                  >  1:4724035/1‑150 (MQ=255)
                                                     cTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCtt                                                                                                 <  1:2111246/150‑1 (MQ=255)
                                                      ttGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTg                                                                                                >  2:3767845/1‑150 (MQ=255)
                                                      ttGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTg                                                                                                >  2:4020531/1‑150 (MQ=255)
                                                        ggTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGaa                                                                                              <  1:4593083/150‑1 (MQ=255)
                                                         gTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAAc                                                                                             <  2:4335005/150‑1 (MQ=255)
                                                            atatCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGAc                                                                                          >  2:4066955/1‑150 (MQ=255)
                                                             tatCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACt                                                                                         >  1:2500647/1‑150 (MQ=255)
                                                               tCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCg                                                                                       >  2:1914663/1‑150 (MQ=255)
                                                               tCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCg                                                                                       >  2:2452860/1‑150 (MQ=255)
                                                               tCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCg                                                                                       >  2:3850283/1‑150 (MQ=255)
                                                               tCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCg                                                                                       >  1:1176750/1‑150 (MQ=255)
                                                                 gAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCt                                                                                     <  2:4686469/150‑1 (MQ=255)
                                                                   aCAATGCCAGGGGGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCACCCACGCCTTGCCCCCCTCCGGTGCCACCGTGGTGTCGAGCACGACTTTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGt                                                                                   <  1:446672/150‑1 (MQ=255)
                                                                        gCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTg                                                                              >  1:4161063/1‑150 (MQ=255)
                                                                            gggTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCAt                                                                          <  1:89530/150‑1 (MQ=255)
                                                                               tGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGa                                                                       >  2:2782249/1‑150 (MQ=255)
                                                                               tGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGa                                                                       >  2:4145920/1‑150 (MQ=255)
                                                                                gCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGaa                                                                      <  2:3525550/150‑1 (MQ=255)
                                                                                  cTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACg                                                                    >  1:3398998/1‑150 (MQ=255)
                                                                                     tgATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACgaga                                                                 <  2:4602068/150‑1 (MQ=255)
                                                                                     tgATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACgaga                                                                 >  1:1486411/1‑150 (MQ=255)
                                                                                      gATGCGGCCTGGCAGGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCCGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGagag                                                                <  1:1458770/150‑1 (MQ=255)
                                                                                       aTGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGc                                                               >  1:4067520/1‑150 (MQ=255)
                                                                                         gCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTg                                                             <  1:1634319/150‑1 (MQ=255)
                                                                                         gCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTg                                                             <  2:1881129/150‑1 (MQ=255)
                                                                                            gCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGgtc                                                          <  1:209146/150‑1 (MQ=255)
                                                                                             ccTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGgtcg                                                         >  1:2309248/1‑150 (MQ=255)
                                                                                               tGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGgtcgtc                                                       >  1:4778485/1‑150 (MQ=255)
                                                                                                ggCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGtcgtcg                                                      >  1:3648198/1‑150 (MQ=255)
                                                                                                 gCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGa                                                     <  2:2217983/150‑1 (MQ=255)
                                                                                                     gTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACttt                                                 >  2:2049410/1‑150 (MQ=255)
                                                                                                     gTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACttt                                                 <  1:2277451/150‑1 (MQ=255)
                                                                                                        gAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTa                                              >  2:5029300/1‑150 (MQ=255)
                                                                                                          gCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTaaa                                            <  2:4826973/150‑1 (MQ=255)
                                                                                                          gCTGGGCGCCCAGCTCCCACGCCTTGCCCTCCTCCGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTaaa                                            <  2:1945984/150‑1 (MQ=255)
                                                                                                           cTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAAc                                           <  2:207108/150‑1 (MQ=255)
                                                                                                                cgcCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATAc                                      >  1:1492838/1‑150 (MQ=255)
                                                                                                                   ccAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCgg                                   >  2:1667344/1‑150 (MQ=255)
                                                                                                                   ccAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCgg                                   >  2:231452/1‑150 (MQ=255)
                                                                                                                   ccAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCgg                                   >  1:1940793/1‑150 (MQ=255)
                                                                                                                      gCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGt                                >  2:3661269/1‑150 (MQ=255)
                                                                                                                      gCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGt                                <  2:369676/150‑1 (MQ=255)
                                                                                                                      gCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGt                                <  1:2723829/150‑1 (MQ=255)
                                                                                                                      gCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGt                                <  2:855303/150‑1 (MQ=255)
                                                                                                                       cTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTg                               >  1:4990279/1‑150 (MQ=255)
                                                                                                                       cTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTg                               >  1:4644759/1‑150 (MQ=255)
                                                                                                                        tCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGa                              <  2:2335014/150‑1 (MQ=255)
                                                                                                                        tCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGa                              <  1:2189109/150‑1 (MQ=255)
                                                                                                                           cACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGt                           >  2:827475/1‑150 (MQ=255)
                                                                                                                           cACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGt                           >  1:103938/1‑150 (MQ=255)
                                                                                                                                 ttGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAg                     <  1:3994212/150‑1 (MQ=255)
                                                                                                                                 ttGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAg                     <  1:1666133/150‑1 (MQ=255)
                                                                                                                                  tcgactcctcGGTGGCACCTGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAgg                    >  2:1986081/5‑150 (MQ=255)
                                                                                                                                    gactcctcGGTGGCACCTGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGcc                  >  2:2359075/3‑150 (MQ=255)
                                                                                                                                    gactcctcGGTGGCACCTGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGcc                  >  1:767955/3‑150 (MQ=255)
                                                                                                                                      ctcctcTGGGGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTg                <  1:3133562/150‑1 (MQ=255)
                                                                                                                                      ctcctcAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTg                >  1:2482539/1‑150 (MQ=255)
                                                                                                                                        tcctcggTGGCACCTGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTgg               >  2:1223920/6‑150 (MQ=255)
                                                                                                                                            aGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACt          >  1:1933030/1‑150 (MQ=255)
                                                                                                                                            aGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACt          >  2:3125364/1‑150 (MQ=255)
                                                                                                                                             ggTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTg         <  1:2923371/150‑1 (MQ=255)
                                                                                                                                              gTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGa        >  2:1643844/1‑150 (MQ=255)
                                                                                                                                                   aCCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCg   <  2:1572748/150‑1 (MQ=255)
                                                                                                                                                   aCCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCg   >  1:4581284/1‑150 (MQ=255)
                                                                                                                                                    ccGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCgg  <  2:1196577/150‑1 (MQ=255)
                                                                                                                                                    ccGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGAGGTCgg  >  1:4940543/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CTGTAGACCGTCTCGCTTGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGG  >  NC_002947/5229104‑5229401

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: