Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 613,255 A→G E173E (GAA→GAG fes → enterobactin/ferrienterobactin esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913613,2550AG100.0% 264.6 / NA 69E173E (GAA→GAGfesenterobactin/ferrienterobactin esterase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (41/28);  total (41/28)

TAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAACGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCA  >  NC_000913/613107‑613401
                                                                                                                                                    |                                                                                                                                                  
tAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCg                                                                                                                                                  >  1:848858/1‑151 (MQ=255)
 aGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCgg                                                                                                                                                 <  1:590195/151‑1 (MQ=255)
 aGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCgg                                                                                                                                                 >  2:66530/1‑151 (MQ=255)
     gATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCTAAAGAAATTATCTGGAAAAGTGAGCGGTTg                                                                                                                                              >  1:521840/1‑150 (MQ=255)
     gATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGa                                                                                                                                             <  1:1172630/151‑1 (MQ=255)
       tCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGaa                                                                                                                                            >  2:674244/1‑150 (MQ=255)
        ccGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGaaaa                                                                                                                                          <  2:600056/151‑1 (MQ=255)
          gCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAg                                                                                                                                        <  2:1017787/151‑1 (MQ=255)
                   aCAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCg                                                                                                                                <  1:24951/150‑1 (MQ=255)
                   aCAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACGCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGc                                                                                                                                 <  1:1230148/149‑1 (MQ=255)
                      aaGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTa                                                                                                                            <  1:574989/151‑1 (MQ=255)
                          tGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATgg                                                                                                                         >  1:1000752/1‑150 (MQ=255)
                                ggcgggcgAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTa                                                                                                                  <  2:895876/151‑1 (MQ=255)
                                  tgggcgAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGGttttttc                                                                                                                 >  1:1077012/2‑148 (MQ=255)
                                  cgggcgAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGGTTTTTacc                                                                                                                >  2:225778/1‑151 (MQ=255)
                                   gggcgAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTacc                                                                                                                >  1:1110785/1‑150 (MQ=255)
                                    ggcgAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTacca                                                                                                               >  2:813049/1‑150 (MQ=255)
                                            cACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCggt                                                                                                       >  1:486251/1‑148 (MQ=255)
                                                 tGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAAc                                                                                                  >  1:1166880/1‑150 (MQ=255)
                                                 tGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTAACACCGGCGATGCaa                                                                                                   >  1:90441/1‑149 (MQ=255)
                                                  gTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAAcag                                                                                                <  1:1031696/151‑1 (MQ=255)
                                                       tGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAAcagca                                                                                              >  2:528283/1‑148 (MQ=255)
                                                        gCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCaga                                                                                            <  2:427478/149‑1 (MQ=255)
                                                        gCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGAAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCaga                                                                                            <  1:907452/149‑1 (MQ=255)
                                                         cACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGGGGCAACACGagaag                                                                                          >  1:595134/1‑150 (MQ=255)
                                                         cACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCagaaga                                                                                         <  2:296528/151‑1 (MQ=255)
                                                              gAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGcc                                                                                    <  2:848896/151‑1 (MQ=255)
                                                                aaTGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGGGATGCAACAACAGAAGAAAGccc                                                                                   >  1:970752/1‑150 (MQ=255)
                                                                   gCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTg                                                                               <  1:393663/151‑1 (MQ=255)
                                                                    ccGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTgg                                                                              <  1:450029/151‑1 (MQ=255)
                                                                      gCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCa                                                                            <  1:814464/151‑1 (MQ=255)
                                                                           cgcCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGGTTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGtttt                                                                       >  2:944414/1‑151 (MQ=255)
                                                                            gcCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGtttt                                                                       >  2:216162/1‑150 (MQ=255)
                                                                                        gggatgggatTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGGTTTTTTCCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCCGGTTTGCTCGATGGCg                                                            >  2:941537/1‑149 (MQ=255)
                                                                                          gatgggatTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCTCCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAAtt                                                        <  2:328698/151‑1 (MQ=255)
                                                                                          gatgggatTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAAtt                                                        <  1:756197/151‑1 (MQ=255)
                                                                                             gggatTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAAtttt                                                      <  2:404240/150‑1 (MQ=255)
                                                                                                     ccGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGGATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAAAGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCaa                                              >  2:87286/1‑150 (MQ=255)
                                                                                                       gCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGGTGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCaaa                                             >  1:983146/1‑149 (MQ=255)
                                                                                                       gCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCCCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCaaa                                             >  2:724679/1‑149 (MQ=255)
                                                                                                        cAAGCGCCAGAAATACCTGCCAAAGAAACTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGGTTTTTACCAACGGGGATGCAACAGCAGAAGAAGGCCCGCGGGCAGTTTTGCTCGGTGGGGGATTTTGGGCGGaaa                                             >  1:240083/1‑148 (MQ=255)
                                                                                                         aaGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGTAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTAt                                         <  1:1089882/151‑1 (MQ=255)
                                                                                                         aaGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGGATTGATTTTTACCAACGGCGGAGCAAAAACAGGAGAACGCCCGCCGGCAGGTTTGGTCGATGGGGAAATTTGGGGGCAAAGTa                                          >  1:701110/1‑150 (MQ=255)
                                                                                                           gcgcCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGCGCAAAGTAAGcc                                      >  1:151737/1‑151 (MQ=255)
                                                                                                               cAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCg                                    >  1:1079347/1‑150 (MQ=255)
                                                                                                                aGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCg                                    >  2:347430/1‑149 (MQ=255)
                                                                                                                  aaaTACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTTGGCGCAAAGTATGCCCGTCt                                 >  2:280691/1‑150 (MQ=255)
                                                                                                                   aaTACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTgg                               <  2:189220/151‑1 (MQ=255)
                                                                                                                    aTACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGGTTTTTAACACCGGCGGTGGCAAAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAAATTTGGGCGGAAAATATGCCCGTCTgg                               >  2:174820/1‑150 (MQ=255)
                                                                                                                     tACCTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGGATGGAATTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCCGGCAGTTTTGCTCGATGGCGGATTTTTGGCGCAAAGTATGCCCGTCTGGcc                             >  2:245163/1‑151 (MQ=255)
                                                                                                                       ccTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCa                            <  2:364438/150‑1 (MQ=255)
                                                                                                                       ccTGCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGGATGCCCGTCTGGcc                             >  1:682481/1‑149 (MQ=255)
                                                                                                                          gCCAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGc                        <  2:413093/151‑1 (MQ=255)
                                                                                                                            cAAAGAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGGATTGATTTTTACCACCGGGGATGCAAAAGAAGAAGAACGCCCGCTGGCGGGTTTGCTCGATTGCGAATTTTGGGGGCAAAGTATGCCCCTCTGGCCCGt                          >  1:653953/1‑147 (MQ=255)
                                                                                                                                gAAATTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCCGCGATGCAACAGCAAAAAAACGCCCCCTGGCAGTTTTTCTCGATGGGGAATTTTTGGCGCAAAGTATGCCCGTCCGGCCAGTGCTGACtt                  >  1:676655/1‑151 (MQ=255)
                                                                                                                                 acaTTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGGATTTTACCACCCGCGATGCAACAGCAGAAGAACGCCCGGTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGAAAGCCCGGCTGGCCAGTGCTGACtt                  >  1:224802/3‑150 (MQ=255)
                                                                                                                                 aaaTTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGTGATGCAACAGCAGAAGAACGCCCGCCGGCAGTTTTGATCCGTCGCGAATTTTAGGCGGAAAATATGCCCGTCCGGCCAGGGCTGACtt                  >  2:445345/1‑150 (MQ=255)
                                                                                                                                  aaTTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTTCCACCGGGGATGCAAAAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACtt                  >  2:896960/1‑149 (MQ=255)
                                                                                                                                  aaTTATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGGATGCCCGTCTGGCCCGTGCTGACt                   >  1:526840/1‑148 (MQ=255)
                                                                                                                                     tATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCt              <  2:1186979/150‑1 (MQ=255)
                                                                                                                                     tATCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTGTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCt              <  2:387830/150‑1 (MQ=255)
                                                                                                                                      aTCTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTg             <  1:225774/150‑1 (MQ=255)
                                                                                                                                        cTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTTCCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGGTTTTCTCGATGGCGAATTTTGGGGGCAAAGTATGCCCGGCTGGCCAGTGCTGACTTCGCTGAcc          >  2:259371/1‑151 (MQ=255)
                                                                                                                                        cTGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTAACCCCGGCGGTGCAACAGCAGAAGAACGCCCCCTGGCAGTTTTTACCGCTGGCGGAATTTGGGCGCCAAGGATGCCCGTCTGGGCCGTTCTGACTTTGCTTAcc          >  1:1095914/1‑151 (MQ=255)
                                                                                                                                         tGGAAAAGTGAGCGGTTGAAAAAGTCACGGCGTGGATGGGGTTTTACCACCGGCGGTGGAACAGAAGAAGAACGCCCGCTGGCAGTTTTTCTCGATGGCGAATTTTGGGCGCCAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGAccc         >  1:508913/1‑151 (MQ=255)
                                                                                                                                            aaaaGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGGAATTTGGGCGCCAAGTATGCCCGTCTGGCCCGTGCTGACTTCCCTGGCCCATc      >  1:778762/1‑151 (MQ=255)
                                                                                                                                            aaaaGTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATc      <  2:1147770/151‑1 (MQ=255)
                                                                                                                                                gTGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTTCCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTTCTCGATGGCGAATTTTTGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGt    >  1:558460/1‑149 (MQ=255)
                                                                                                                                                 tGAGCGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCa  <  2:213100/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAACGGTTGAAAAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCA  >  NC_000913/613107‑613401

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