Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 172,913 G→A 6.2% M484I (ATG→ATA fhuB → iron(3+)‑hydroxamate import ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913172,9130GA 6.2% 2391.2 / 94.8 755M484I (ATG→ATAfhuBiron(3+)‑hydroxamate import ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (363/343);  new base A (26/21);  total (391/364)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.78e-01

Only 100 of 766 total aligned reads displayed.

TGGTGTTGATGCTGTTTCTGGTGCCGGGTAATGCCTTTGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGAATTGTGATGGTGATTT  >  NC_000913/172765‑173061
                                                                                                                                                    |                                                                                                                                                    
          gCTGTTTCTGGTGCCGGGTAATGCCTTTGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCa                                                                                                                                            >  2:2997340/1‑149 (MQ=255)
          gCTGTTTCTGGTGCCGGGTAATGCCTTTGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCa                                                                                                                                            >  2:378400/1‑149 (MQ=255)
          gCTGTTTCTGGTGCCGGGTAATGCCTTTGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCa                                                                                                                                            >  4:547308/1‑149 (MQ=255)
            tGTTTCTGGTGCCGGGTAATGCCTTTGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCAcc                                                                                                                                          >  1:6797377/1‑149 (MQ=255)
              tttCTGGTGCCGGGTAATGCCTTTGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACcgc                                                                                                                                        >  3:1258660/1‑149 (MQ=255)
              tttCTGGTGCCGGGTAATGCCTTTGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACcgc                                                                                                                                        >  2:5501832/1‑149 (MQ=255)
               ttCTGGTGCCGGGTAATGCCTTTGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACcgcg                                                                                                                                       >  3:2853999/1‑149 (MQ=255)
               ttCTGGTGCCGGGTAATGCCTTTGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACcgcg                                                                                                                                       >  2:3677416/1‑149 (MQ=255)
                           gTAATGCCTTTGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCtt                                                                                                                           >  2:3368587/1‑149 (MQ=255)
                            tAATGCCTTTGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCttt                                                                                                                          >  1:848877/1‑149 (MQ=255)
                               tGCCTTTGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTtga                                                                                                                       >  4:3646874/1‑149 (MQ=255)
                               tGCCTTTGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTtga                                                                                                                       <  3:193470/149‑1 (MQ=255)
                               tGCCTTTGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATAGCGTTAAGCACCGCGTTCACCATGCTTTtga                                                                                                                       <  3:5265512/149‑1 (MQ=255)
                                gCCTTTGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGc                                                                                                                                             <  2:1601499/126‑1 (MQ=255)
                                gCCTTTGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTtgat                                                                                                                      >  1:6260637/1‑149 (MQ=255)
                                  cTTTGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTc                                                                                                                                    <  4:1797143/133‑1 (MQ=255)
                                  cTTTGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTtgatga                                                                                                                    <  4:1258660/149‑1 (MQ=255)
                                  cTTTGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGAGTTCACCATGCTTTtgatga                                                                                                                    >  3:5856978/1‑149 (MQ=255)
                                    ttGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGatgatg                                                                                                                  <  2:6337174/149‑1 (MQ=255)
                                    ttGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGatgatg                                                                                                                  <  2:5196608/149‑1 (MQ=255)
                                       gctggctgTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTg                                                                                                               >  4:385054/1‑149 (MQ=255)
                                        ctggctgTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATATTGc                                                                                                              >  4:6163451/1‑149 (MQ=255)
                                         tggctgTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGtt                                                                                                                <  1:1701948/146‑1 (MQ=255)
                                         tggctgTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCa                                                                                                             >  1:971253/1‑149 (MQ=255)
                                         tggctgTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATAGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCa                                                                                                             >  2:5105195/1‑149 (MQ=255)
                                          ggctgTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAg                                                                                                            <  3:2569108/149‑1 (MQ=255)
                                          ggctgTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAg                                                                                                            <  3:4991985/149‑1 (MQ=255)
                                             tgTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCa                                                                                                         >  1:1249192/1‑149 (MQ=255)
                                             tgTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCa                                                                                                         <  1:1543414/149‑1 (MQ=255)
                                              gTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCaa                                                                                                        <  3:5575602/149‑1 (MQ=255)
                                                 aCCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGttga                                                                                                                                               >  3:639740/1‑105 (MQ=255)
                                                  ccTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAgtgg                                                                                                    <  1:351657/149‑1 (MQ=255)
                                                  ccTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAgt                                                                                                      >  3:3756008/1‑147 (MQ=255)
                                                   cTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAgtggt                                                                                                   <  3:2304140/149‑1 (MQ=255)
                                                   cTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAgt                                                                                                      <  1:6612183/146‑1 (MQ=255)
                                                     gcaggcagTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGtt                                                                                                                <  4:3064777/134‑1 (MQ=255)
                                                        ggcagTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGAccc                                                                                              <  4:2149125/149‑1 (MQ=255)
                                                        ggcagTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGAccc                                                                                              >  3:7118850/1‑149 (MQ=255)
                                                         gcagTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCcg                                                                                             >  2:1180307/1‑149 (MQ=255)
                                                             tctcGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGaa                                                                                         <  4:4966249/149‑1 (MQ=255)
                                                             tctcGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGaa                                                                                         <  1:263462/149‑1 (MQ=255)
                                                              ctcGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAAt                                                                                        <  4:2541754/149‑1 (MQ=255)
                                                               tcGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATg                                                                                       <  1:5186942/149‑1 (MQ=255)
                                                               tcGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATg                                                                                       <  2:4717139/149‑1 (MQ=255)
                                                                 ggcgcggcgGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGgc                                                                                     <  2:2234795/149‑1 (MQ=255)
                                                                          gTGACGCTGTTGATCAATATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCt                                                                            <  4:5790990/149‑1 (MQ=255)
                                                                           tGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTg                                                                           <  1:5846074/149‑1 (MQ=255)
                                                                             aCGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGAc                                                                         >  4:897320/1‑149 (MQ=255)
                                                                               gCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCt                                                                       >  1:7144830/1‑149 (MQ=255)
                                                                                 tGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTgg                                                                     >  2:6673091/1‑149 (MQ=255)
                                                                                   ttGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGAt                                                                   <  3:6644317/149‑1 (MQ=255)
                                                                                    tGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGAtt                                                                  <  4:6000476/149‑1 (MQ=255)
                                                                                     gATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCcgc                                                                                            >  4:6683799/1‑122 (MQ=255)
                                                                                      aTCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGc                                                                             >  3:996099/1‑136 (MQ=255)
                                                                                      aTCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTc                                                                <  2:1925348/149‑1 (MQ=255)
                                                                                         aTTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCgg                                                             >  4:4081706/1‑149 (MQ=255)
                                                                                         aTTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCgg                                                             <  3:317864/149‑1 (MQ=255)
                                                                                          ttATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAgtgg                                                                                                    <  3:880204/109‑1 (MQ=255)
                                                                                          ttATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGt                                                            >  2:7237544/1‑149 (MQ=255)
                                                                                          ttATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGt                                                            >  3:350675/1‑149 (MQ=255)
                                                                                            aTGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCTCCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTc                                                          >  2:623000/1‑149 (MQ=255)
                                                                                            aTGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATAGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTc                                                          <  2:4010128/149‑1 (MQ=255)
                                                                                                  gccgccGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTa                                                    <  4:4750572/149‑1 (MQ=255)
                                                                                                      ccggccgCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAAc                                                <  3:4636155/149‑1 (MQ=255)
                                                                                                       cggccgCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAAcg                                               >  1:2959138/1‑149 (MQ=255)
                                                                                                           cgcgGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATAGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGAc                                           >  1:4538723/1‑149 (MQ=255)
                                                                                                              ggTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGa                                        >  1:2388387/1‑149 (MQ=255)
                                                                                                               gTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGc                                                                                                              <  3:1684586/78‑1 (MQ=255)
                                                                                                               gTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGAt                                       >  3:1824488/1‑149 (MQ=255)
                                                                                                                    ttttCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCa                                  >  2:4158140/1‑149 (MQ=255)
                                                                                                                    ttttCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCa                                  >  1:3185183/1‑149 (MQ=255)
                                                                                                                    ttttCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCa                                  >  1:5557228/1‑149 (MQ=255)
                                                                                                                     ttttCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTcc                                                               >  1:3913257/1‑119 (MQ=255)
                                                                                                                       tCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGt                               <  1:4199919/149‑1 (MQ=255)
                                                                                                                         cccACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCt                             >  3:4206272/1‑149 (MQ=255)
                                                                                                                          ccACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGAt                                                                   <  2:4350630/110‑1 (MQ=255)
                                                                                                                          ccACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTg                            >  3:6194285/1‑149 (MQ=255)
                                                                                                                          ccACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTg                            >  4:3303764/1‑149 (MQ=255)
                                                                                                                          ccACACCGTATGTTACTGGCGGGGATAGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGa                                            >  1:1862448/1‑133 (MQ=255)
                                                                                                                          ccAAACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTg                            >  1:1026322/1‑149 (MQ=255)
                                                                                                                           cacaCCGTATGTTACTGGCGGGGATAGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCCGGTCTgg                           >  2:1372633/1‑149 (MQ=255)
                                                                                                                              acCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGgcgc                        >  4:6280717/1‑149 (MQ=255)
                                                                                                                                cGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGCGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCAc                      >  1:3761520/1‑149 (MQ=255)
                                                                                                                                 gTATGTTACTGGCGGGGATAGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCAcc                     >  4:5242395/1‑149 (MQ=255)
                                                                                                                                    tGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGa                  >  1:4911142/1‑149 (MQ=255)
                                                                                                                                     gTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCggga                 >  2:1761231/1‑147 (MQ=255)
                                                                                                                                     gTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGaa                 <  2:2541319/149‑1 (MQ=255)
                                                                                                                                     gTTACTGGCGGGGATAGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCaa                                                                                 <  2:3153334/85‑1 (MQ=255)
                                                                                                                                       tACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGAAtt               <  2:855074/149‑1 (MQ=255)
                                                                                                                                         cTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGAATtgt             <  4:1511225/149‑1 (MQ=255)
                                                                                                                                         cTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGAATtgt             <  3:1491074/149‑1 (MQ=255)
                                                                                                                                         cTGGCGGGGATAGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGAATtgt             <  4:647916/149‑1 (MQ=255)
                                                                                                                                          tGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGAATtgtg            >  2:6000978/1‑149 (MQ=255)
                                                                                                                                           ggCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGAATTGTGa           <  3:6008815/149‑1 (MQ=255)
                                                                                                                                            gCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGAATTGTGAt          <  1:2004735/149‑1 (MQ=255)
                                                                                                                                             cGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGAATTGTGATg         <  2:7045253/149‑1 (MQ=255)
                                                                                                                                             cGGGGATAGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATCTTGCAGGCAAGTGGTGACCCGCGAAGGGCGCAATTGCTGACCTGGATTTGCGGCTCGACCTACAACGCGACCGATGCGCAGGTATGGCGCACCGGAATTgtggtg         <  4:4434152/149‑1 (MQ=255)
                                                                                                                                               gggATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGAATTGTGATGGt       <  4:6825963/149‑1 (MQ=255)
                                                                                                                                                 gATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTc                                                          >  3:6809557/1‑96 (MQ=255)
                                                                                                                                                   tGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGAATTGTGATGGTGAtt   >  2:6808845/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGGTGTTGATGCTGTTTCTGGTGCCGGGTAATGCCTTTGGCTGGCTGTTACCTGCAGGCAGTCTCGGCGCGGCGGTGACGCTGTTGATCATTATGATCGCCGCCGGCCGCGGTGGATTTTCCCCACACCGTATGTTACTGGCGGGGATGGCGTTAAGCACCGCGTTCACCATGCTTTTGATGATGTTGCAGGCAAGTGGTGACCCGCGAATGGCGCAAGTGCTGACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGAATTGTGATGGTGATTT  >  NC_000913/172765‑173061

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