Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 619,445 A→G V252A (GTT→GCT)  fepC ← ferrienterobactin ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913619,4450AG99.0% 385.2 / ‑3.9 102V252A (GTT→GCT) fepCferrienterobactin ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (42/59);  total (42/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 102 total aligned reads displayed.

CGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAATGGCTGGCGTAACGACAGGCCTGATTAA  >  NC_000913/619297‑619589
                                                                                                                                                    |                                                                                                                                                
cggcggGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCt                                                                                                                                                >  1:794142/1‑151 (MQ=255)
  gcggGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTgg                                                                                                                                              <  2:1235461/151‑1 (MQ=255)
   cggGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTgg                                                                                                                                              >  1:1228524/1‑150 (MQ=255)
   cggGATGGTGGCTTGTGGGGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTgg                                                                                                                                              >  2:620568/1‑150 (MQ=255)
     ggATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATc                                                                                                                                           >  2:860352/1‑151 (MQ=255)
     ggATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGCTCCGGCAGCTGGATc                                                                                                                                           <  2:226549/151‑1 (MQ=255)
         ggtggCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGTTCGTCa                                                                                                                                       <  1:385819/151‑1 (MQ=255)
                tgtgGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGAt                                                                                                                                  <  2:745650/149‑1 (MQ=255)
                  tggtggCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATg                                                                                                                              >  2:1209598/1‑151 (MQ=255)
                   ggtggCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATAATg                                                                                                                              <  2:658455/150‑1 (MQ=255)
                       gCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAgcg                                                                                                                         >  1:657351/1‑151 (MQ=255)
                        cGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAgcgc                                                                                                                        <  2:892250/151‑1 (MQ=255)
                         gTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCa                                                                                                                       <  2:1323718/151‑1 (MQ=255)
                           gCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAg                                                                                                                      <  1:1393582/150‑1 (MQ=255)
                            cTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGAc                                                                                                                    >  2:241584/1‑151 (MQ=255)
                              gTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGAcc                                                                                                                   >  2:1234416/1‑150 (MQ=255)
                              gTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGAc                                                                                                                    <  2:1587070/149‑1 (MQ=255)
                                 gcgcTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATaa                                                                                                               <  2:414753/151‑1 (MQ=255)
                                 gcgcTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATa                                                                                                                <  1:826251/150‑1 (MQ=255)
                                  cgcTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATaa                                                                                                               >  2:1524488/1‑150 (MQ=255)
                                      aTGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATgcg                                                                                                          <  2:1014023/151‑1 (MQ=255)
                                      aTGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATg                                                                                                            <  2:1259814/149‑1 (MQ=255)
                                        gCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCt                                                                                                        <  2:538494/151‑1 (MQ=255)
                                          gATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCa                                                                                                      <  1:151743/151‑1 (MQ=255)
                                          gATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCa                                                                                                      <  1:1502145/151‑1 (MQ=255)
                                           aTGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCaa                                                                                                     <  1:247567/151‑1 (MQ=255)
                                             ggCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAAt                                                                                                    <  2:666616/150‑1 (MQ=255)
                                              gCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCa                                                                                                  <  1:117334/151‑1 (MQ=255)
                                               cATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATc                                                                                                   >  1:378731/1‑149 (MQ=255)
                                                aTCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGt                                                                                                <  1:883261/151‑1 (MQ=255)
                                                   cAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTc                                                                                              >  2:1538547/1‑150 (MQ=255)
                                                     gAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTcagc                                                                                           <  1:833893/151‑1 (MQ=255)
                                                      aaaaCAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTcagca                                                                                          <  2:1248148/151‑1 (MQ=255)
                                                      aaaaCAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTcagca                                                                                          <  1:1626654/151‑1 (MQ=255)
                                                      aaaaCAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTcagc                                                                                           >  1:462068/1‑150 (MQ=255)
                                                      aaaaCAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTcagc                                                                                           >  2:507585/1‑150 (MQ=255)
                                                       aaaCAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTcagcag                                                                                         <  1:253995/151‑1 (MQ=255)
                                                       aaaCAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTcagc                                                                                           >  2:355200/1‑149 (MQ=255)
                                                        aaCAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGt                                                                                        <  1:388249/151‑1 (MQ=255)
                                                         aCAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTg                                                                                       >  2:681887/1‑151 (MQ=255)
                                                            ggATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACa                                                                                    <  1:706332/151‑1 (MQ=255)
                                                            ggATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACa                                                                                    <  1:534634/151‑1 (MQ=255)
                                                             gATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACaa                                                                                   <  1:1582955/151‑1 (MQ=255)
                                                             gATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACaa                                                                                   <  1:600600/151‑1 (MQ=255)
                                                              aTGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAAt                                                                                  <  1:860284/151‑1 (MQ=255)
                                                                gCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAAtct                                                                                <  2:758243/151‑1 (MQ=255)
                                                                   atgatgGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCt                                                                             >  1:1069612/1‑151 (MQ=255)
                                                                   atgatgGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCt                                                                             <  1:1012636/151‑1 (MQ=255)
                                                                    tgatgGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCtt                                                                            >  2:1592434/1‑151 (MQ=255)
                                                                    tgatgGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCt                                                                             <  1:967684/150‑1 (MQ=255)
                                                                        ggCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCgg                                                                         <  1:685060/150‑1 (MQ=255)
                                                                        ggCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCgg                                                                         <  2:149411/150‑1 (MQ=255)
                                                                              ccACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTccc                                                                  >  1:1541655/1‑151 (MQ=255)
                                                                               cACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCt                                                                 <  2:1431599/151‑1 (MQ=255)
                                                                                 cTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTg                                                                >  2:707326/1‑150 (MQ=255)
                                                                                   tAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGa                                                               <  2:731817/149‑1 (MQ=255)
                                                                                   tAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGa                                                               <  2:342264/149‑1 (MQ=255)
                                                                                   tAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAg                                                              <  1:606532/150‑1 (MQ=255)
                                                                                    aGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGca                                                            <  1:546426/151‑1 (MQ=255)
                                                                                    aGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGca                                                            <  2:896351/151‑1 (MQ=255)
                                                                                      tttAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGcaac                                                          >  1:847399/1‑151 (MQ=255)
                                                                                         aaCTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAAt                                                       >  1:1097119/1‑151 (MQ=255)
                                                                                          aCTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAAtt                                                      >  1:390483/1‑151 (MQ=255)
                                                                                          aCTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAAtg                                                      <  1:241562/151‑2 (MQ=255)
                                                                                             aTTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAAtttt                                                    <  2:1185196/150‑1 (MQ=255)
                                                                                             aTTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTttc                                                   >  2:1449022/1‑151 (MQ=255)
                                                                                             aTTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTttc                                                   >  1:260363/1‑151 (MQ=255)
                                                                                             aTTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTttc                                                   >  2:211346/1‑151 (MQ=255)
                                                                                               ttGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTttccc                                                 >  2:35457/1‑151 (MQ=255)
                                                                                                   ggTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTttcccttcc                                             >  2:485032/1‑151 (MQ=255)
                                                                                                   ggTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTttcccttcc                                             <  1:181676/151‑1 (MQ=255)
                                                                                                     ttGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCg                                           >  1:1651015/1‑151 (MQ=255)
                                                                                                       gAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCgca                                         <  1:399043/151‑1 (MQ=255)
                                                                                                       gAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCgc                                          <  2:657410/150‑1 (MQ=255)
                                                                                                          tgTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCgcaatg                                      >  2:1044028/2‑151 (MQ=255)
                                                                                                            tGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCgcaatgc                                     >  1:1169640/1‑150 (MQ=255)
                                                                                                               cGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATc                                 >  2:1561889/1‑151 (MQ=255)
                                                                                                                gTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCa                                <  1:589367/151‑1 (MQ=255)
                                                                                                                 ttCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCa                                >  2:180112/1‑150 (MQ=255)
                                                                                                                 ttCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCa                                >  1:373542/1‑150 (MQ=255)
                                                                                                                   cGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCaaa                              <  1:426666/150‑1 (MQ=255)
                                                                                                                    gTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAAtg                            >  2:1500666/1‑151 (MQ=255)
                                                                                                                    gTCCGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAAtg                            >  2:359431/1‑151 (MQ=255)
                                                                                                                      ccGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAAtggc                          >  1:143880/1‑151 (MQ=255)
                                                                                                                       cGAGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAAtggct                         <  1:699661/151‑1 (MQ=255)
                                                                                                                         aGCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAAtggctgg                       <  2:687000/151‑1 (MQ=255)
                                                                                                                          gCGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAAtggctggc                      <  1:1506494/151‑1 (MQ=255)
                                                                                                                           cGGCACCACAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAATGGCTGGCg                     >  1:1340395/1‑151 (MQ=255)
                                                                                                                              caccacAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAATGGCTGGCGTaa                  <  2:1317171/151‑1 (MQ=255)
                                                                                                                              caccacAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAATGGCTGGCGTa                   >  2:1499504/1‑150 (MQ=255)
                                                                                                                              caccacAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAATGGCTGGCGTa                   >  1:1702156/1‑150 (MQ=255)
                                                                                                                               accacAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAATGGCTGGCGTAAc                 <  2:177548/151‑1 (MQ=255)
                                                                                                                                ccacAAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAATGGCTGGCGTAACg                >  2:759974/1‑151 (MQ=255)
                                                                                                                                 cacaAGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAATGGCTGGCGTAACGa               <  1:553458/151‑1 (MQ=255)
                                                                                                                                    aaGCGGCGTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAATGGCTGGCGTAACGACAg            <  1:1209505/151‑1 (MQ=255)
                                                                                                                                       cggcgTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAATGGCTGGCGTAACGACAGGcc         <  2:1541751/151‑1 (MQ=255)
                                                                                                                                           gTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAATGGCTGGCGTAACGACAGGCCTGAt     >  1:758723/1‑151 (MQ=255)
                                                                                                                                           gTTCCGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAATGGCTGGCGTAACGACAGGCCTGAt     <  1:1327197/151‑1 (MQ=255)
                                                                                                                                              ccGGCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAATGGCTGGCGTAACGACAGGCCTGATTaa  <  2:1499898/151‑1 (MQ=255)
                                                                                                                                                ggCAGCTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAATGGCTGGCGTAACGACAGGCCTGATTa   >  2:908286/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAATGGCTGGCGTAACGACAGGCCTGATTAA  >  NC_000913/619297‑619589

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