Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,559,989 C→A 6.6% intergenic (+14/‑50) fepA → / → PP_2243 ferric enterobactin transport system outer membrane subunit/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,559,9890CA 6.6% 216.9 / 4.3 69intergenic (+14/‑50)fepA/PP_2243ferric enterobactin transport system outer membrane subunit/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (38/26);  new base A (2/3);  total (40/29)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

CGTGCGGGTAGGCGTGAACAACGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATG  >  NC_002947/2559846‑2560135
                                                                                                                                                |                                                                                                                                                  
cgTGCGGGTAGGCGTGAACAACGTCTTCGACAAGCAGATCCTGCGCGGCGGCACCGCCAGCAGCTCGGGGGCAAAGACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCCCTGAATATCAACTTCTCATCCCAAGGGCCGCCCTGCg                                                                                                                                               >  3:303342/1‑149 (MQ=255)
  tGCGGGTAGGCGTGAACAACGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGc                                                                                                                                             <  8:225124/149‑1 (MQ=255)
       gTAGGCGTGAACAACGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACACCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCAGGGCGGcccttg                                                                                                                                        <  8:320041/149‑4 (MQ=255)
        tAGGCGTGAACAACGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCAGGGCGGcccttgg                                                                                                                                       <  5:164125/149‑5 (MQ=255)
           gCGTGAACAACGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATcgcg                                                                                                                                    <  7:56702/149‑1 (MQ=255)
             gTGAACAACGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGt                                                                                                                                  >  6:119152/1‑149 (MQ=255)
                aacaacGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTc                                                                                                                                >  4:227753/1‑148 (MQ=255)
                   aacGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATc                                                                                                                            >  3:384808/1‑149 (MQ=255)
                         ttCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTcg                                                                                                                      >  1:311554/1‑149 (MQ=255)
                            gACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTcgcg                                                                                                                    >  1:388174/1‑148 (MQ=255)
                             aCAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGaa                                                                                                                  >  5:318155/1‑149 (MQ=255)
                                 gCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGc                                                                                                              >  5:369160/1‑149 (MQ=255)
                                  cAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGcc                                                                                                             <  2:3588/149‑1 (MQ=255)
                                     aTCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCcgac                                                                                                          >  8:41030/1‑149 (MQ=255)
                                     aTCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCcgac                                                                                                          >  3:84960/1‑149 (MQ=255)
                                      tCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCcgacg                                                                                                         <  6:318155/149‑1 (MQ=255)
                                           cgcgGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACAtt                                                                                                    <  4:384808/149‑1 (MQ=255)
                                               gcggcAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTg                                                                                                >  6:298958/1‑149 (MQ=255)
                                                cggcAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGc                                                                                               <  5:16674/149‑1 (MQ=255)
                                                  gcAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAg                                                                                             >  5:203890/1‑149 (MQ=255)
                                                     aCGCCAGCAGCTCGGGGGCAAACACGTATAAGCCAGCCGGGCAGGGCGGTGTTTGATGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTCCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGc                                                                                            >  1:41584/1‑148 (MQ=255)
                                                       gCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGCGGACATGCTTGCAGCCaaa                                                                                        >  6:122898/1‑149 (MQ=255)
                                                        ccAGCAGCTCGGGGGCAAGCACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAAc                                                                                       >  5:270033/1‑149 (MQ=255)
                                                        ccAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAAc                                                                                       <  3:227753/149‑1 (MQ=255)
                                                        ccAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAAc                                                                                       <  5:119152/149‑1 (MQ=255)
                                                        ccAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTACAGCCAAAc                                                                                       <  7:275437/149‑1 (MQ=255)
                                                          agcagcTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGt                                                                                     <  1:179339/149‑1 (MQ=255)
                                                                  gggggCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCGGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCa                                                                               <  2:41584/147‑1 (MQ=255)
                                                                    gggCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGccc                                                                                                                                                   <  3:358275/77‑1 (MQ=255)
                                                                    gggCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGccc                                                                                                                                                   >  4:358275/1‑77 (MQ=255)
                                                                    gggCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGt                                                                           <  4:275962/149‑1 (MQ=255)
                                                                    cggCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCCACTTCAGATCCCAAGGGCCGCCCTGCGGCCCAAACGCGGGACCATCGGGACGCGAATGGCCGACGACATTCTTGCAGGCAAACGTCTAGCAGCgg                                                                           >  8:150406/2‑148 (MQ=255)
                                                                      gCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCcgacga                                                                                                        <  3:67023/118‑1 (MQ=255)
                                                                      gCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCcgacga                                                                                                        >  4:67023/1‑118 (MQ=255)
                                                                      gCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGGGc                                                                         >  3:339428/1‑149 (MQ=255)
                                                                          acacGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCCGCCAAACGTCTCGCAGCGTGCTaaa                                                                     >  7:184532/1‑149 (MQ=255)
                                                                          acacGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGcc                                                                                           <  6:69795/127‑1 (MQ=255)
                                                                          acacGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGcc                                                                                           >  5:69795/1‑127 (MQ=255)
                                                                          acacGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTaaa                                                                     >  1:181944/1‑149 (MQ=255)
                                                                              gTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCAt                                                                 >  4:377638/1‑149 (MQ=255)
                                                                              gTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCAt                                                                 >  2:282055/1‑149 (MQ=255)
                                                                              gTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCAGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCAt                                                                 >  3:353498/1‑149 (MQ=255)
                                                                               tataACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATg                                                                <  2:388174/149‑1 (MQ=255)
                                                                               tataACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGAAGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATg                                                                >  3:220695/1‑149 (MQ=255)
                                                                                ataACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATg                                                                >  1:15375/1‑148 (MQ=255)
                                                                                       aGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCAccc                                                           >  4:330826/1‑147 (MQ=255)
                                                                                          cGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGcc                                                      >  7:85366/1‑149 (MQ=255)
                                                                                                  gCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCga                                              <  2:181944/149‑1 (MQ=255)
                                                                                                        tttGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATg                                        >  8:30601/1‑149 (MQ=255)
                                                                                                            ctgcGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTg                                    <  3:377638/149‑1 (MQ=255)
                                                                                                            ctgcGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTg                                    <  2:311554/149‑1 (MQ=255)
                                                                                                              gcgcTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGt                                  >  6:140425/1‑149 (MQ=255)
                                                                                                                  tGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGaaa                              >  4:261619/1‑149 (MQ=255)
                                                                                                                  tGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGa                                >  5:14137/1‑147 (MQ=255)
                                                                                                                    aaTATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGaaa                              <  2:15375/147‑1 (MQ=255)
                                                                                                                        tCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTg                                    <  3:150217/137‑1 (MQ=255)
                                                                                                                        tCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTg                                    >  4:150217/1‑137 (MQ=255)
                                                                                                                        tCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATc                        >  7:3979/1‑149 (MQ=255)
                                                                                                                        tCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATc                        <  7:30601/149‑1 (MQ=255)
                                                                                                                             ttCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCt                   >  4:281812/1‑149 (MQ=255)
                                                                                                                             ttCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCt                   >  4:276479/1‑149 (MQ=255)
                                                                                                                               cTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCt                   <  8:85366/147‑1 (MQ=255)
                                                                                                                               cTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCt                   <  6:203890/147‑1 (MQ=255)
                                                                                                                                  aTCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCCGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGAccc              >  1:16011/1‑149 (MQ=255)
                                                                                                                                    cccAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCt             >  5:18166/1‑148 (MQ=255)
                                                                                                                                     gattgggCCGCAGGGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGAGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGa           >  2:329235/5‑149 (MQ=255)
                                                                                                                                     ccAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGa           <  3:330826/149‑1 (MQ=255)
                                                                                                                                      attgggCCGCAGGGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTg                                                 <  2:21158/107‑1 (MQ=255)
                                                                                                                                      attgggCCGCAGGGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTg                                                 >  1:21158/4‑110 (MQ=255)
                                                                                                                                        aGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCt        <  6:14137/149‑1 (MQ=255)
                                                                                                                                        aGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCt        <  5:298958/149‑1 (MQ=255)
                                                                                                                                              gccCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATg  >  3:362859/1‑149 (MQ=255)
                                                                                                                                              gccCTGCGGCCCAATCGCGGTGACATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATg  <  6:270033/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CGTGCGGGTAGGCGTGAACAACGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATG  >  NC_002947/2559846‑2560135

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