Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,512,350:1 +G 5.1% coding (65/957 nt) fecD ← ferric citrate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,512,3501.G 5.1% 606.5 / 26.6 136coding (65/957 nt)fecDferric citrate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base . (67/62);  new base G (4/3);  total (71/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

Only 100 of 145 total aligned reads displayed.

GCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGAT  >  NC_000913/4512078‑4512604
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                              
gcGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATgg                                                                                                                                                                                                                                                              <  2:415291/275‑1 (MQ=255)
                    aGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGcc                                                                                                                                                                                                                                          >  1:120641/1‑275 (MQ=255)
                                                 agagGCCAGGCTGGCGGCATGGTTAACTCCGGGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAAc                                                                                                                                                                                                               <  2:427874/273‑1 (MQ=255)
                                                     gCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACCTAATAATGCCCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCCCGCCAAGGCACGGGGATCCCCCCCAAAGGGCGTGATAACAGCCCAACACCCGCCAGGGCAAGGGTGACGAAAAACACCAgag                                                                                                                                                                                                          >  4:390879/1‑272 (MQ=255)
                                                                               gagaAGATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACaa                                                                                                                                                                                <  2:182482/274‑1 (MQ=255)
                                                                                agaATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACCGTGCCAGCCGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATCCTCGTGTCCGGCCTGCCAGTCGGTCACCAGCGCGCGCCACGGCACGGGGAACACTCCCAACTGGGGTGGTACCAGCGCACAGCCCCCCCGCTCCAGGGTGATCAAAATACCCAGGCCAATTTTCCTCCTCtcctcccaccacac                                                                                                                                                                              >  4:19753/1‑267 (MQ=255)
                                                                                agaATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTACACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCCCAGCCCGCCAGGGCAAGGGTGCTGAAAATACCCAGCTCAATTTTCATCCTCGCCTCCTCCCAAGc                                                                                                                                                                              >  2:94206/1‑275 (MQ=255)
                                                                                                                      ggggaacAGCACGCCCGGCAGGGGGGGGGCTACACCGCAAAAACGGGCGAGCAGCAAGCGGGGCAGTCGGTACTCANTCAATACATAATAAGGCGCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGAACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTg                                                                                                                                         <  3:120641/268‑1 (MQ=255)
                                                                                                                                   cccgccACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGAAg                                                                                                                            >  2:69492/1‑274 (MQ=255)
                                                                                                                                      gccACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATa                                                                                                                                                                                                                                                       <  6:201774/148‑1 (MQ=255)
                                                                                                                                                  gCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGcac                                                                                                                                                                                                                                              <  5:852045/145‑1 (MQ=255)
                                                                                                                                                     gCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCg                                                                                                                                                                                                                                        >  8:26994/1‑148 (MQ=255)
                                                                                                                                                      cACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCACATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGAAAAAGCAAGGGCAGCCAATCAGCGCCAGCACCGCCACCGCCGGCAGAACCCAGGGGAAGGCCAgac                                                                                                        >  4:114593/1‑273 (MQ=255)
                                                                                                                                                      cACCGAAAAACAGTGCCGCGCGAAAGGGTGGCAGTTGGTACCACCGCATTGCAAAAGAATGCCCGGTGCCGGGCTGCCAGTGGGACCGCAGCGGGCGCCACGGGACGGGGATCACTCCCATATGGGGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAgcg                                                                                                        <  4:261558/275‑1 (MQ=255)
                                                                                                                                                            cAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAAc                                                                                                                                                                                                                                                     >  7:1121795/1‑128 (MQ=255)
                                                                                                                                                                   tGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCaa                                                                                                                                                                                                                               >  7:316823/1‑143 (MQ=255)
                                                                                                                                                                   tGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGCTCTCCGGGGACGGCCAGCGCGCGTGCGAGTa                                                                                            >  2:186941/1‑274 (MQ=255)
                                                                                                                                                                   tGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTAc                                                                                           >  2:453435/1‑275 (MQ=255)
                                                                                                                                                                       agcagcAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGtgatga                                                                                                                                                                                                                      >  7:748011/1‑148 (MQ=255)
                                                                                                                                                                        gcagcaAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGa                                                                                                                                                                                                                                                             >  5:143761/1‑108 (MQ=255)
                                                                                                                                                                        gcagcaAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGa                                                                                                                                                                                                                                                             <  6:143761/108‑1 (MQ=255)
                                                                                                                                                                            caAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAg                                                                                                                                                                                                                              <  7:664013/135‑1 (MQ=255)
                                                                                                                                                                             aaGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCTCCAGGGCAAGGGTGATGAAAATAACCAGNGCAATTTTCCTCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCNNATc                                                                                                                                                      >  4:543026/1‑206 (MQ=255)
                                                                                                                                                                               gcgcGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAAcc                                                                                                                                                                                                              >  6:507752/1‑148 (MQ=255)
                                                                                                                                                                                cgcgGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCATATgggagt                                                                                                                                                                                                                                                          >  5:66756/1‑98 (MQ=255)
                                                                                                                                                                                     cAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTGGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCa                                                                                                                                                                                                        >  6:963128/1‑148 (MQ=255)
                                                                                                                                                                                        tCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAAtt                                                                                                                                                                                                     <  7:26994/148‑1 (MQ=255)
                                                                                                                                                                                           gTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTg                                                                     <  2:483950/273‑1 (MQ=255)
                                                                                                                                                                                            tACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCaa                                                                                                                                                                                                                               <  7:783691/118‑1 (MQ=255)
                                                                                                                                                                                            tACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAAc                                                                                                                                                                                                               >  6:345446/1‑134 (MQ=255)
                                                                                                                                                                                            tACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAAc                                                                                                                                                                                                               <  5:345446/134‑1 (MQ=255)
                                                                                                                                                                                            tACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCa                                                                                                                                                                                                 <  8:748011/148‑1 (MQ=255)
                                                                                                                                                                                                 cATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCCCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCCTCAGCGTGGCCCCcag                                                             >  2:473064/1‑275 (MQ=255)
                                                                                                                                                                                                  aTCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCaa                                                                                                                                                                                                       >  6:1054701/1‑136 (MQ=255)
                                                                                                                                                                                                  aTCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCaa                                                                                                                                                                                                       <  5:1054701/136‑1 (MQ=255)
                                                                                                                                                                                                  aTCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAAtttt                                                                                                                                                                                                   >  5:1072030/1‑140 (MQ=255)
                                                                                                                                                                                                   tCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGc                                                                                                                                                                                          >  6:249509/1‑148 (MQ=255)
                                                                                                                                                                                                    cAATACATAATAATGCTGGTGTACGGCCTGCCAGGCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCcagca                                                           <  2:401848/274‑1 (MQ=255)
                                                                                                                                                                                                    cAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGcc                                                                                                                                                                                         >  8:424575/1‑148 (MQ=255)
                                                                                                                                                                                                      aTACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGcctc                                                                                                                                                                                       >  6:752692/1‑148 (MQ=255)
                                                                                                                                                                                                          ataataATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCt                                                                                                                                    <  2:259395/195‑1 (MQ=255)
                                                                                                                                                                                                          ataataATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTg                                                                                                                                   <  2:518175/196‑1 (MQ=255)
                                                                                                                                                                                                           taataaTGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTcac                                                                                                                                                                                  >  8:207133/1‑148 (MQ=255)
                                                                                                                                                                                                           taataaTGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAg                                                                                                                                                                     >  4:20042/1‑161 (MQ=255)
                                                                                                                                                                                                           taataaTGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCATATggga                                                                                                                                                                                                                                                            >  8:329654/1‑71 (MQ=255)
                                                                                                                                                                                                           taataaTGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCATATggga                                                                                                                                                                                                                                                            <  7:329654/73‑3 (MQ=255)
                                                                                                                                                                                                            aataatGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCga                                                                                                                                                                                                                                       <  5:292043/94‑2 (MQ=255)
                                                                                                                                                                                                             ctaatGCTCTTGTCCGGCCTGCCACACGGTCAGCAGCGCGAGCTACGGCACGGGGATCACTCCCATATGGAGTGCTAACAGCGCAAAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCAc                                                 <  2:293618/274‑1 (MQ=255)
                                                                                                                                                                                                               aatGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAg                                                                                                                                                                 <  2:20042/161‑1 (MQ=255)
                                                                                                                                                                                                                 tGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAg                                                                                                                                                                                                                              <  7:695801/98‑1 (MQ=255)
                                                                                                                                                                                                                 tGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAg                                                                                                                                                                                                                              >  8:695801/1‑98 (MQ=255)
                                                                                                                                                                                                                 tGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCCGCACTGCGCCTGCGGGCCGATCTCCGGGGAAGGCCAGCGCCCGTGCGCGTACATCTGCCCGCAGCCTCAGCGTCGCCCCCCGCCGCATGCTCAc                                                 >  1:147510/1‑271 (MQ=255)
                                                                                                                                                                                                                 tGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGGGTGCGAGTACACCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCCCGGGc                                             >  2:458486/1‑275 (MQ=255)
                                                                                                                                                                                                                 tGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCCCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTCCGAGTACAACTGCCAGCAGCATCAGGGTGGCCCCCAGCAGACTGCCCACTGGc                                             >  4:501464/1‑275 (MQ=255)
                                                                                                                                                                                                                   cTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGa                                                                                                                                                                          >  5:44665/1‑148 (MQ=255)
                                                                                                                                                                                                                   cTCGTGTCCGGCCTGCCAGTCGGTCAGAAGCGCGCGCCACGGCACGGGGATAACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACACAGCAAGGGCTGCCAAGCAGCGCCAGCCCTGCGCCTGCGGGAAGATCTACGGGGAAGGCCAGCGCGCGTGCGCGTACATCTGACAGCAGCATCAGCGGGGCACCCACCCGCCTGCAAACTggga                                            >  1:132023/1‑272 (MQ=255)
                                                                                                                                                                                                                       tgtCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAg                                                                                                                                                                                                                              >  5:460229/1‑92 (MQ=255)
                                                                                                                                                                                                                          ccGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCa                                                                                                                                                                   <  5:963128/148‑1 (MQ=255)
                                                                                                                                                                                                                          ccGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCa                                                                                                                                                                   >  7:995934/1‑148 (MQ=255)
                                                                                                                                                                                                                            ggCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCa                                                                                                                                                                                                 <  6:944605/116‑1 (MQ=255)
                                                                                                                                                                                                                            ggCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATcc                                                                                                                                                                                              >  8:872629/1‑119 (MQ=255)
                                                                                                                                                                                                                            ggCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATcc                                                                                                                                                                                              <  7:872629/119‑1 (MQ=255)
                                                                                                                                                                                                                             gCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGcc                                                                                                                                                                             >  8:428756/1‑135 (MQ=255)
                                                                                                                                                                                                                             gCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGcc                                                                                                                                                                             <  7:428756/135‑1 (MQ=255)
                                                                                                                                                                                                                             gCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCa                                                                                                                                                                   >  7:568895/1‑145 (MQ=255)
                                                                                                                                                                                                                             gCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGCACTCCGGGGAAGGCCAGCGCGCGTTCGAGTAAATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTgcc                                  >  2:289370/1‑273 (MQ=255)
                                                                                                                                                                                                                             gCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACACCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTTCTCACAAGCCAGACAAAGCa                                                                                                                                                                   <  8:568895/145‑1 (MQ=255)
                                                                                                                                                                                                                                    aGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCaa                                                                                                                                                                                                       >  8:942256/1‑102 (MQ=255)
                                                                                                                                                                                                                                    aGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCaa                                                                                                                                                                                                       <  7:942256/102‑1 (MQ=255)
                                                                                                                                                                                                                                    aGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTg                                                                                                                                                            <  6:144443/145‑1 (MQ=255)
                                                                                                                                                                                                                                         gTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAg                                                                                                                                                    <  6:44665/148‑1 (MQ=255)
                                                                                                                                                                                                                                           cagcagCGCGCGCCAGGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGtga                                                                                                                                                                                                                         <  6:737816/77‑1 (MQ=255)
                                                                                                                                                                                                                                           cagcagCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGtga                                                                                                                                                                                                                         >  5:737816/1‑77 (MQ=255)
                                                                                                                                                                                                                                                  gcgcgcCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAg            <  2:543027/275‑1 (MQ=255)
                                                                                                                                                                                                                                                   cgcgcCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCaa                                                                                                                                                                                                                               <  5:358719/63‑1 (MQ=255)
                                                                                                                                                                                                                                                       ccACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACaaa                                                                                                                                                                      >  7:203706/1‑116 (MQ=255)
                                                                                                                                                                                                                                                       ccACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACaaa                                                                                                                                                                      <  8:203706/116‑1 (MQ=255)
                                                                                                                                                                                                                                                         agggcagggGGAGCACCCCCATAGTGGGGGATAAAAGCCAACAGCCCGCAAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCGGGCCTCCTCACAAGCCAGACAAGGCAAGGGCGGCCAATCAGCGCCAGCACGGCGCCGGCGGGCAGATGTCCGGGGAAGGCAAGAGCGCGGGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGcc      <  4:134634/272‑1 (MQ=255)
                                                                                                                                                                                                                                                         acggcacggGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTg                                                                                                                                         >  7:258515/1‑143 (MQ=255)
                                                                                                                                                                                                                                                         acggcacggGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTg                                                                                                                                         <  8:258515/143‑1 (MQ=255)
                                                                                                                                                                                                                                                         acggcacggGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTg                                                                                                                                         >  5:981440/1‑143 (MQ=255)
                                                                                                                                                                                                                                                         acggcacggGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTg                                                                                                                                         <  6:981440/143‑1 (MQ=255)
                                                                                                                                                                                                                                                               cggGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCaa                                                                                                                                                        <  8:771584/122‑1 (MQ=255)
                                                                                                                                                                                                                                                                ggggATCACTCACATATGGAGTGATAACAGCGCACAGCACGCCAGGACAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAg                                                                                                                                                                 <  5:190992/112‑1 (MQ=255)
                                                                                                                                                                                                                                                                 gggATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTgcg                                                                                                                                       >  6:32407/1‑137 (MQ=255)
                                                                                                                                                                                                                                                                 gggATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTgcg                                                                                                                                       <  5:32407/137‑1 (MQ=255)
                                                                                                                                                                                                                                                                 gggATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAg                                                                                                                            <  7:246404/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                   taTCACTCCATATGGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTg                                                                                                                                                            >  8:1020284/2‑114 (MQ=255)
                                                                                                                                                                                                                                                                   gATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGcc                                                                                                                                                                             >  7:739016/1‑97 (MQ=255)
                                                                                                                                                                                                                                                                   gATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGcc                                                                                                                                                                             <  6:1011259/97‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aTCACTCCATATGGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTAATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAg                                                                                                                                                                 >  8:705612/1‑108 (MQ=255)
                                                                                                                                                                                                                                                                     tCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAAt                                                                                                                                                                                                      >  6:417676/1‑70 (MQ=255)
                                                                                                                                                                                                                                                                      cACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCa                                                                                                                                                     >  6:984340/1‑118 (MQ=255)
                                                                                                                                                                                                                                                                      cACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGAtctc                                                                                                                       >  8:546685/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                        cTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCa                                                                                                                                                                   >  5:829570/1‑102 (MQ=255)
                                                                                                                                                                                                                                                                        cTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCa                                                                                                                                                                   <  6:829570/102‑1 (MQ=255)
                                                                                                                                                                                                                                                                        actccATATGGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGcc                                                                                                                                                                             >  5:185386/4‑93 (MQ=255)
                                                                                                                                                                                                                                                                        actccATATGGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGcc                                                                                                                                                                             <  6:185386/90‑1 (MQ=255)
                                                                                                                                                                                                                                                                          cccATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAggg                                                                                                                                                                                                                            >  6:509241/1‑43 (MQ=255)
                                                                                                                                                                                                                                                                          cccATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAggg                                                                                                                                                                                                                            <  5:509241/43‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                              
GCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGAT  >  NC_000913/4512078‑4512604

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: