Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 619,419 G→A 9.7% L261F (CTC→TTC)  fepC ← ferrienterobactin ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913619,4190GA 9.7% 267.8 / 13.3 93L261F (CTC→TTC) fepCferrienterobactin ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (37/46);  new base A (3/6);  total (41/52)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.61e-01

ATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAATGGCTGG  >  NC_000913/619271‑619568
                                                                                                                                                    |                                                                                                                                                     
attattGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGc                                                                                                                                                      <  2:1290515/150‑1 (MQ=255)
     tgggATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGAAGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGCTCGGACGAACGGcac                                                                                                                                                 <  1:155709/147‑1 (MQ=255)
      gtgATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGcacc                                                                                                                                                <  2:1514862/150‑1 (MQ=255)
      gtgATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGCTCGGACGAACGGcacc                                                                                                                                                <  2:625130/150‑1 (MQ=255)
      gtgATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTAGCGCTATGAGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTCAAGGTGCATTTCCTCCGATCGGcacc                                                                                                                                                >  1:32095/1‑150 (MQ=255)
      gtgATCCTTTCCGCGTTGATCGGCGGGATGGTGGCGTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGCTCGGACGAACGGcacc                                                                                                                                                <  1:760413/150‑1 (MQ=255)
       tgATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGcacca                                                                                                                                               >  1:505979/1‑150 (MQ=255)
       tgATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGCTCGGACGAACGGcacct                                                                                                                                               <  2:603385/150‑2 (MQ=255)
       tgATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGCTCGGACGAACGGcacct                                                                                                                                               <  2:160980/150‑2 (MQ=255)
          tCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACaa                                                                                                                                            >  1:1202353/1‑150 (MQ=255)
              ttCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGGGCTGTTGCGCTATGCGTTGTCTTCAATAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAgcgg                                                                                                                                        <  2:470565/150‑1 (MQ=255)
               tCCGCGTTGATCGGCGGGATGGTGGCTTGTGGGGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAgcggc                                                                                                                                       <  2:550598/150‑1 (MQ=255)
                ccGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAgcggcg                                                                                                                                      >  2:1296767/1‑150 (MQ=255)
                 cgcgTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGt                                                                                                                                     >  2:1610129/1‑150 (MQ=255)
                   cgTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTc                                                                                                                                   >  2:679504/1‑150 (MQ=255)
                   cgTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTc                                                                                                                                   >  2:841028/1‑150 (MQ=255)
                     gtGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCg                                                                                                                                 <  1:1644545/149‑1 (MQ=255)
                       gATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGc                                                                                                                               >  2:138550/1‑150 (MQ=255)
                        aTCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGACGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCa                                                                                                                              >  1:897058/1‑150 (MQ=255)
                               ggATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTAAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCACCTGGAt                                                                                                                       >  1:1564506/1‑150 (MQ=255)
                                 aTGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCg                                                                                                                     <  1:1573317/150‑1 (MQ=255)
                                   ggtggCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTc                                                                                                                   <  2:1700141/150‑1 (MQ=255)
                                        cTTGTGGTGGCGTGCTGTTGCGCTATGAGATGGCATACAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGa                                                                                                              >  1:1128921/1‑150 (MQ=255)
                                             ggtggCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATg                                                                                                         >  2:1004729/1‑150 (MQ=255)
                                             ggtggCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATg                                                                                                         >  2:1410884/1‑150 (MQ=255)
                                             ggtggCGTGATGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCATACCACTTATTTAAACTATTTGAGGTTGAAGGTGCCGATCGTCCGAGCGGAACCACAAGCGGAGTTCCGTCAACTGGATCGTCAATGATCAtc                                                                                                         >  2:90666/1‑149 (MQ=255)
                                                ggCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAg                                                                                                      >  2:1459273/1‑150 (MQ=255)
                                                   gTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAgcgc                                                                                                   <  2:1620286/150‑1 (MQ=255)
                                                   gTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAgcgc                                                                                                   >  2:121506/1‑150 (MQ=255)
                                                     gCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAg                                                                                                 <  2:669643/150‑1 (MQ=255)
                                                        gTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGAcc                                                                                              >  1:83730/1‑150 (MQ=255)
                                                         ttGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCa                                                                                             <  1:886663/150‑1 (MQ=255)
                                                           gcgcTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGTGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATa                                                                                           <  1:411612/150‑1 (MQ=255)
                                                            cgcTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATaa                                                                                          >  1:943520/1‑150 (MQ=255)
                                                              cTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGATGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATaaaa                                                                                        >  1:1437477/1‑149 (MQ=255)
                                                               tATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATg                                                                                       <  2:622749/150‑1 (MQ=255)
                                                                   cGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCt                                                                                   >  2:633935/1‑150 (MQ=255)
                                                                          aTCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAg                                                                            <  2:51034/150‑1 (MQ=255)
                                                                          aTCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAg                                                                            >  2:1718958/1‑150 (MQ=255)
                                                                            ccAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGtt                                                                          >  1:1339680/1‑150 (MQ=255)
                                                                            ccAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGtt                                                                          <  1:1647710/150‑1 (MQ=255)
                                                                             cAGAAAACAGGATGCCATGATGGCAGCCCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTCCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATATATGCGCTCAATCAGTTc                                                                         >  2:952365/1‑150 (MQ=255)
                                                                             cAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGtta                                                                         >  1:346753/1‑149 (MQ=255)
                                                                             cAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTc                                                                         >  2:378534/1‑150 (MQ=255)
                                                                             cAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTc                                                                         >  2:1190723/1‑150 (MQ=255)
                                                                             cAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTc                                                                         <  2:63640/150‑1 (MQ=255)
                                                                                aaaaCAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTcagc                                                                      <  1:955712/150‑1 (MQ=255)
                                                                                 aaaCAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGCGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTcagca                                                                     >  1:1200945/1‑150 (MQ=255)
                                                                                      ggATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGAc                                                                <  1:919326/150‑1 (MQ=255)
                                                                                      ggATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGAc                                                                <  1:1662237/150‑1 (MQ=255)
                                                                                      ggATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGAc                                                                >  1:365333/1‑150 (MQ=255)
                                                                                         tGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAAt                                                             >  1:139203/1‑150 (MQ=255)
                                                                                         tGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAACGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAAt                                                             <  2:104639/150‑1 (MQ=255)
                                                                                           ccATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGGGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAAtct                                                           <  1:656376/150‑1 (MQ=255)
                                                                                            cATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAAtctc                                                          >  2:1359728/1‑150 (MQ=255)
                                                                                             atgatgGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTcc                                                         <  2:505979/150‑1 (MQ=255)
                                                                                              tgatgGCAGCCCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGCCCATAAATGCGGTAAATCAGTTCAGCAGTGACAATATCCt                                                        >  2:1239886/1‑150 (MQ=255)
                                                                                                   gCAGACCACTTAGGTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGAGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGgc                                                   <  1:290565/150‑1 (MQ=255)
                                                                                                    cAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGgcg                                                  <  2:1078803/150‑1 (MQ=255)
                                                                                                      taCCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCt                                                <  2:20062/149‑1 (MQ=255)
                                                                                                          aCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCt                                            <  1:1515746/150‑1 (MQ=255)
                                                                                                           cTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTg                                           <  2:401481/150‑1 (MQ=255)
                                                                                                            ttAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCTCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGa                                          <  2:1128921/150‑1 (MQ=255)
                                                                                                            ttAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGa                                          <  1:1273593/150‑1 (MQ=255)
                                                                                                              aGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGc                                        <  2:1202353/150‑1 (MQ=255)
                                                                                                              aGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGc                                        <  1:382985/150‑1 (MQ=255)
                                                                                                                 ttAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGcaac                                     >  1:467382/1‑150 (MQ=255)
                                                                                                                  tAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGcaaca                                    <  1:121506/150‑1 (MQ=255)
                                                                                                                   aaCTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTc                                                                         <  2:3602/112‑1 (MQ=255)
                                                                                                                   aaCTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTc                                                                         >  1:3602/1‑112 (MQ=255)
                                                                                                                   aaCTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGcaacaa                                   >  2:428581/1‑150 (MQ=255)
                                                                                                                   aaCTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGcaacaa                                   <  1:633935/150‑1 (MQ=255)
                                                                                                                        tttGCTGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCAGGCAACTGGATCGTCAATGTTCATGAAGCGCAGACCATAAATGCGCTCAATCAGTTGTGCAGTGACAATCTCCTTCTGCGCTCCCTGAGCAACAATTttc                              <  2:583934/150‑1 (MQ=255)
                                                                                                                        tttGCGGTTGAAGGTGCCGTTCGTCCAAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTttc                              <  2:146353/150‑1 (MQ=255)
                                                                                                                        tttGCGGTTGAAGGGGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTGCGGCGCTCCCTGAGCAACAATTttc                              <  1:1011287/150‑1 (MQ=255)
                                                                                                                         ttGCGGTTGACGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGGTCCGGCAACTGGATCGTCAATGATCATTCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTttcc                             <  1:290957/150‑1 (MQ=255)
                                                                                                                         ttGCGGTTGAAGGTGCCGTTCGTCCGCGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATAATCATGCAGCGCAGACCATAAATGCGCTCAATCAGATCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTttcc                             <  1:138550/150‑1 (MQ=255)
                                                                                                                         ttGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTttcc                             <  1:1086013/150‑1 (MQ=255)
                                                                                                                         ttGCGGGTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTttcc                             <  1:1146876/150‑1 (MQ=255)
                                                                                                                           gCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTttccct                           <  1:181174/150‑1 (MQ=255)
                                                                                                                             cgcttgtggTGCCGCTCGGACGAACGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTttcccttc                         >  1:1322723/8‑150 (MQ=255)
                                                                                                                              gTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTttcccttcc                        >  1:1717651/1‑150 (MQ=255)
                                                                                                                              gTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTttcccttcc                        >  1:1717649/1‑150 (MQ=255)
                                                                                                                              gTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTttcccttcc                        >  2:712033/1‑150 (MQ=255)
                                                                                                                                 gAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCACTCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCgc                     >  1:944392/1‑150 (MQ=255)
                                                                                                                                   tggTGCCGCTCGGACGAACGGCACCACAAGCGGAGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATTAATGCGCTCAATCAGTTCAGCAGAGCCAATCTCCTTCGGCGCACCCTGAGCAACACTTTTCCCTTCCCgcaa                   >  1:351161/2‑150 (MQ=255)
                                                                                                                                     gTGCCGCTCGGACGAACGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGCCCATAAATGCGCGCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCgcaatg                 >  1:1184971/1‑150 (MQ=255)
                                                                                                                                      tGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCgcaatgc                <  2:1339680/150‑1 (MQ=255)
                                                                                                                                      tGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCgcaatgc                <  2:706059/150‑1 (MQ=255)
                                                                                                                                             cGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCaaa         <  2:83730/150‑1 (MQ=255)
                                                                                                                                              gTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAAt        <  1:1459273/150‑1 (MQ=255)
                                                                                                                                               tCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAAtg       <  1:378534/150‑1 (MQ=255)
                                                                                                                                                  gAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCCTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAAtggct    <  1:413582/150‑1 (MQ=255)
                                                                                                                                                    gCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAAtggctgg  <  1:1323734/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAATGGCTGG  >  NC_000913/619271‑619568

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