Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 622,455 G→C 6.1% W52C (TGG→TGC entS → enterobactin exporter, iron‑regulated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913622,4550GC 6.1% 347.3 / 11.9 98W52C (TGG→TGCentSenterobactin exporter, iron‑regulated
Reads supporting (aligned to +/- strand):  ref base G (45/47);  new base C (3/3);  total (48/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

ATCCTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACTGTGCCTT  >  NC_000913/622311‑622602
                                                                                                                                                |                                                                                                                                                   
aTCCTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCACGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAg                                                                                                                                                  >  4:478781/1‑148 (MQ=255)
   cTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGCTGATGACACATTCCACCTGGCAGGTg                                                                                                                                               <  2:972689/148‑1 (MQ=255)
    tGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTgg                                                                                                                                              <  1:758936/148‑1 (MQ=255)
     ggCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTggg                                                                                                                                             >  3:134205/1‑148 (MQ=255)
      gctgctTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTgggg                                                                                                                                            <  4:134205/148‑1 (MQ=255)
                cTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGa                                                                                                                                  >  3:791752/1‑148 (MQ=255)
                 tCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGAc                                                                                                                                 >  3:802452/1‑148 (MQ=255)
                 tCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGCCAGGTgga                                                                                                                                             <  4:999576/136‑2 (MQ=255)
                 tCAGCCTGTTCAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGCCAGGTgga                                                                                                                                             >  3:999576/1‑135 (MQ=255)
                  cAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAg                                                                                                                                                  <  2:1068337/130‑1 (MQ=255)
                  cAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAg                                                                                                                                                  >  1:1068337/1‑130 (MQ=255)
                      cTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGa                                                                                                                            >  4:697663/1‑148 (MQ=255)
                             aaaCGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACcggcgg                                                                                                                     >  2:106682/1‑148 (MQ=255)
                                    cccGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAtg                                                                                                              >  1:907837/1‑148 (MQ=255)
                                    cccGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAtg                                                                                                              >  1:814027/1‑148 (MQ=255)
                                      cGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAtgtt                                                                                                            <  2:878616/148‑1 (MQ=255)
                                         cGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAtgtttgt                                                                                                         <  2:814027/148‑1 (MQ=255)
                                         cGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAtgtttgt                                                                                                         <  2:907837/148‑1 (MQ=255)
                                         cGTTTCACGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAtgtttgt                                                                                                         <  3:478781/148‑1 (MQ=255)
                                            ttCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGCCAGGTgga                                                                                                                                             <  2:49318/109‑2 (MQ=255)
                                            ttCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGCCAGGTgga                                                                                                                                             >  1:49318/1‑108 (MQ=255)
                                                 gcAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGa                                                                                                 <  4:802452/148‑1 (MQ=255)
                                                 gcAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGa                                                                                                 <  4:791752/148‑1 (MQ=255)
                                                      aTTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAgg                                                                                                                                                 >  2:846894/1‑95 (MQ=255)
                                                      aTTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAgg                                                                                                                                                 <  1:846894/95‑1 (MQ=255)
                                                        tCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTg                                                                                                                                               >  1:512950/1‑95 (MQ=255)
                                                        tCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTg                                                                                                                                               <  2:512950/95‑1 (MQ=255)
                                                          ctcgctcgTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGc                                                                                                     <  1:98300/135‑1 (MQ=255)
                                                          ctcgctcgTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGc                                                                                                     >  2:98300/1‑135 (MQ=255)
                                                            cgctcgTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGCCAGGTgga                                                                                                                                             >  1:1100013/1‑92 (MQ=255)
                                                            cgctcgTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGCCAGGTgga                                                                                                                                             <  2:1100013/93‑2 (MQ=255)
                                                                   ttCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCa                                                                                                                                                   >  1:651618/1‑80 (MQ=255)
                                                                   ttCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCa                                                                                                                                                   <  2:651618/80‑1 (MQ=255)
                                                                       tctcAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTg                                                                                                                                   <  2:1008595/92‑1 (MQ=255)
                                                                       tctcAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTg                                                                                                                                   >  1:1008595/1‑92 (MQ=255)
                                                                        ctcAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGACGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGAt                                                                          >  1:727114/1‑148 (MQ=255)
                                                                         tcAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATc                                                                         >  3:932420/1‑148 (MQ=255)
                                                                          cAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTcggc                                                                                         >  2:858406/1‑131 (MQ=255)
                                                                          cAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTcggc                                                                                         <  1:858406/131‑1 (MQ=255)
                                                                           aaTTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGc                                                                       <  1:106682/148‑1 (MQ=255)
                                                                              tgtgTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGc                                                                                  <  4:19231/134‑1 (MQ=255)
                                                                              tgtgTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGc                                                                                  >  3:19231/1‑134 (MQ=255)
                                                                                   ctctGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAgcg                                                               <  3:697663/148‑1 (MQ=255)
                                                                                     ctGGGTTTGCTCTGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCa                                                             >  2:709028/1‑148 (MQ=255)
                                                                                     ctGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTa                                                                     <  2:78288/140‑1 (MQ=255)
                                                                                     ctGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTa                                                                     >  1:78288/1‑140 (MQ=255)
                                                                                       gggTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAt                                                                                                               <  1:188569/96‑1 (MQ=255)
                                                                                       gggTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAt                                                                                                               >  2:188569/1‑96 (MQ=255)
                                                                                       gggTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTttgttg                                                                                                       >  2:407320/1‑104 (MQ=255)
                                                                                       gggTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTttgttg                                                                                                       <  1:407320/104‑1 (MQ=255)
                                                                                        ggTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTa                                                                     >  3:438971/1‑137 (MQ=255)
                                                                                        ggTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTa                                                                     <  4:438971/137‑1 (MQ=255)
                                                                                        ggTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCaaaa                                                          >  3:731867/1‑148 (MQ=255)
                                                                                         gTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCttt                                                                                                                                        >  3:1141747/1‑69 (MQ=255)
                                                                                         gTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCttt                                                                                                                                        <  4:1141747/69‑1 (MQ=255)
                                                                                         gTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCaaaaa                                                         <  2:727114/148‑1 (MQ=255)
                                                                                          tttGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTc                                                                                                                                       >  3:452012/1‑69 (MQ=255)
                                                                                          tttGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTc                                                                                                                                       <  4:452012/69‑1 (MQ=255)
                                                                                           ttGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCa                                                                                                                                                   <  2:991144/56‑1 (MQ=255)
                                                                                           ttGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCa                                                                                                                                                   >  1:991144/1‑56 (MQ=255)
                                                                                              ctctGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGa                                                    <  1:709028/148‑1 (MQ=255)
                                                                                               tCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTg                                                                                                                                               >  4:329240/1‑56 (MQ=255)
                                                                                               tCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTg                                                                                                                                               <  3:329240/56‑1 (MQ=255)
                                                                                                cGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCg                                                                                                                                      >  1:941232/1‑64 (MQ=255)
                                                                                                cGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCg                                                                                                                                      <  2:941232/64‑1 (MQ=255)
                                                                                                     tCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCt                                             <  4:731867/148‑1 (MQ=255)
                                                                                                        cggtgccggtgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGcgcg                                                                                                                 <  2:639162/77‑1 (MQ=255)
                                                                                                        cggtgccggtgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGcgcg                                                                                                                 >  1:639162/1‑77 (MQ=255)
                                                                                                        cggtgccggtgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATg                                                                   >  3:671002/1‑123 (MQ=255)
                                                                                                        cggtgccggtgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATg                                                                   <  4:671002/123‑1 (MQ=255)
                                                                                                         ggtgccggtgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTc                                                                                                                                       >  4:913974/1‑54 (MQ=255)
                                                                                                         ggtgccggtgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTc                                                                                                                                       <  3:913974/54‑1 (MQ=255)
                                                                                                             ccggtgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGcgcgcg                                     <  4:932420/148‑1 (MQ=255)
                                                                                                                 tgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTg                                                                                <  1:24374/101‑1 (MQ=255)
                                                                                                                 tgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTg                                                                                >  2:24374/1‑101 (MQ=255)
                                                                                                                      aTCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCtgtg                            >  3:282989/1‑148 (MQ=255)
                                                                                                                       tCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGAttt                                                 >  3:52944/1‑126 (MQ=255)
                                                                                                                       tCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGAttt                                                 <  4:52944/126‑1 (MQ=255)
                                                                                                                            atgatgACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCt                                                                      >  3:116036/1‑100 (MQ=255)
                                                                                                                            atgatgACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCt                                                                      <  4:116036/100‑1 (MQ=255)
                                                                                                                            atgatgACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCAt                        <  1:1122758/146‑1 (MQ=255)
                                                                                                                            atgatgACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCAt                        >  2:1122758/1‑146 (MQ=255)
                                                                                                                               atgaCACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGCGCTgg                                                                               <  2:624757/88‑1 (MQ=255)
                                                                                                                               atgaCACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGCGCTgg                                                                               >  1:624757/1‑88 (MQ=255)
                                                                                                                                   cacaTTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAtgtttgt                                                                                                         <  4:102014/58‑1 (MQ=255)
                                                                                                                                   cacaTTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAtgtttgt                                                                                                         >  3:102014/1‑58 (MQ=255)
                                                                                                                                   cacaTTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCAt               >  2:427122/1‑148 (MQ=255)
                                                                                                                                     caTTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAtgt                                                                                                             >  3:917990/1‑52 (MQ=255)
                                                                                                                                     caTTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAtgt                                                                                                             <  4:917990/52‑1 (MQ=255)
                                                                                                                                       ttCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGa           <  1:427122/148‑1 (MQ=255)
                                                                                                                                        tCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTg                                                     <  1:688410/105‑1 (MQ=255)
                                                                                                                                        tCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTg                                                     >  2:688410/1‑105 (MQ=255)
                                                                                                                                            ccTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATTGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACtgtg      >  1:937703/1‑148 (MQ=255)
                                                                                                                                            ccTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATTGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACtgtg      >  1:19700/1‑148 (MQ=255)
                                                                                                                                            ccTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAgcg                                                               <  4:349814/91‑1 (MQ=255)
                                                                                                                                            ccTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAgcg                                                               >  3:349814/1‑91 (MQ=255)
                                                                                                                                             cTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATTGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACTGTGc     <  2:937703/148‑1 (MQ=255)
                                                                                                                                             cTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATTGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACTGTGc     <  2:19700/148‑1 (MQ=255)
                                                                                                                                                gCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACTGTGCCtt  >  3:417409/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
ATCCTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACTGTGCCTT  >  NC_000913/622311‑622602

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