Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 622,455 G→C 6.8% 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.8% 359.2 / 14.1 103W52C (TGG→TGCentSenterobactin exporter, iron‑regulated
Reads supporting (aligned to +/- strand):  ref base G (50/46);  new base C (3/4);  total (53/50)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.68e-01

Only 100 of 107 total aligned reads displayed.

AACAATCCTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACTGTGCCTTA  >  NC_000913/622307‑622603
                                                                                                                                                    |                                                                                                                                                    
aaCAATCCTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTgg                                                                                                                                                      >  4:714041/1‑149 (MQ=255)
 aCAATCCTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTgg                                                                                                                                                      >  2:589188/1‑148 (MQ=255)
 aCAATCCTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGc                                                                                                                                                     >  4:547467/1‑149 (MQ=255)
   aaTCCTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAg                                                                                                                                                   <  4:287872/149‑1 (MQ=255)
      ccTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGt                                                                                                                                                 >  2:101950/1‑148 (MQ=255)
       cTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTg                                                                                                                                                <  2:795989/148‑1 (MQ=255)
          gctgctTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGCCAGGTggaa                                                                                                                                             <  2:206758/148‑3 (MQ=255)
             gctTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCtt                                                                                                                                          >  1:125835/1‑148 (MQ=255)
              ctTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCttt                                                                                                                                         <  1:338520/148‑1 (MQ=255)
               ttAACCTCAGCGTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCTGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCg                                                                                                                                       >  3:794272/1‑149 (MQ=255)
               ttAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTc                                                                                                                                        <  2:360049/148‑1 (MQ=255)
               ttAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCCTTCg                                                                                                                                       <  3:414922/149‑1 (MQ=255)
                tAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCg                                                                                                                                       <  2:125835/148‑1 (MQ=255)
                   ccTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGa                                                                                                                                   >  3:682894/1‑149 (MQ=255)
                     tCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGAc                                                                                                                                  >  2:423610/1‑148 (MQ=255)
                      cAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACg                                                                                                                                 >  1:442354/1‑148 (MQ=255)
                        gCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCAGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTg                                                                                                                              <  4:191036/149‑1 (MQ=255)
                             ttGGAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACcgg                                                                                                                         >  4:393752/1‑149 (MQ=255)
                                   aCGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGgcgc                                                                                                                   >  3:580233/1‑149 (MQ=255)
                                     gCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGcgcg                                                                                                                  <  1:671456/148‑1 (MQ=255)
                                     gCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGa                                                                                                                 <  3:547467/149‑1 (MQ=255)
                                      cACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGa                                                                                                                 >  2:387752/1‑148 (MQ=255)
                                      cACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGa                                                                                                                 <  2:430753/148‑1 (MQ=255)
                                            gCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGCtgtttg                                                                                                           <  2:499411/148‑1 (MQ=255)
                                             cGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAtgtttgt                                                                                                          <  1:675940/148‑1 (MQ=255)
                                             cGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAtgtttgt                                                                                                          <  2:162822/148‑1 (MQ=255)
                                               tttCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGCCAGGTggaa                                                                                                                                             <  1:455154/111‑3 (MQ=255)
                                               ttgcgcgCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGCCAGGTggaa                                                                                                                                             >  2:455154/4‑109 (MQ=255)
                                                   gcgcAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTg                                                                                                   >  4:267769/1‑149 (MQ=255)
                                                       aGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATg                                                                                                >  1:437939/1‑148 (MQ=255)
                                                              ctcgctcgTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTcggcgg                                                                                        >  3:83679/1‑149 (MQ=255)
                                                               tcgctcgTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGCCAGGTggaa                                                                                                                                             <  3:247037/95‑3 (MQ=255)
                                                                            ctcAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGAt                                                                           >  1:888472/1‑148 (MQ=255)
                                                                             tcAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGa                                                                                                                                   <  3:683276/91‑1 (MQ=255)
                                                                             tcAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGa                                                                                                                                   >  4:683276/1‑91 (MQ=255)
                                                                             tcAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCg                                                                         >  3:840648/1‑149 (MQ=255)
                                                                             tcAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCg                                                                         >  3:30437/1‑149 (MQ=255)
                                                                              cAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTg                                                                                                                                    >  2:265289/1‑89 (MQ=255)
                                                                              cAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTg                                                                                                                                    <  1:265289/89‑1 (MQ=255)
                                                                              cAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCg                                                                         >  1:247175/1‑148 (MQ=255)
                                                                              cAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCg                                                                         >  1:749713/1‑148 (MQ=255)
                                                                                 ttGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTa                                                                      <  2:749713/148‑1 (MQ=255)
                                                                                 ttGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTAt                                                                     >  4:764871/1‑149 (MQ=255)
                                                                                   gtgtCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGa                                                                   <  4:840648/149‑1 (MQ=255)
                                                                                   gtgtCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGa                                                                   <  4:207682/149‑1 (MQ=255)
                                                                                   gtgtCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGa                                                                   <  3:139281/149‑1 (MQ=255)
                                                                                     gtCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAgc                                                                 <  3:418088/149‑1 (MQ=255)
                                                                                     gtCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAg                                                                  <  1:167563/148‑1 (MQ=255)
                                                                                        tctGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCa                                                              >  4:517657/1‑149 (MQ=255)
                                                                                           gggTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCaaa                                                            >  1:191066/1‑148 (MQ=255)
                                                                                             gTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCaaaaa                                                          <  1:101950/148‑1 (MQ=255)
                                                                                             gTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCaaaaa                                                          <  1:240638/148‑1 (MQ=255)
                                                                                             gTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCaaaaa                                                          <  2:91276/148‑1 (MQ=255)
                                                                                                 gCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGa                                                     <  3:159359/149‑1 (MQ=255)
                                                                                                    cGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGAtt                                                   >  1:422897/1‑148 (MQ=255)
                                                                                                         tCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTg                                             <  4:83679/149‑1 (MQ=255)
                                                                                                           gcgGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGgc                                           <  4:146972/149‑1 (MQ=255)
                                                                                                            cggtgccggtgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGgc                                           >  1:158693/1‑148 (MQ=255)
                                                                                                             ggtgccggtgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGgcg                                          <  2:422897/148‑1 (MQ=255)
                                                                                                             ggtgccggtgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGgcg                                          <  2:185512/148‑1 (MQ=255)
                                                                                                             ggtgccggtgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGgcg                                          <  1:589188/148‑1 (MQ=255)
                                                                                                                gccggtgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGcgcgcg                                      >  3:743210/1‑149 (MQ=255)
                                                                                                                  cggtgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGAtttt                                                 >  2:825558/1‑136 (MQ=255)
                                                                                                                  cggtgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGAtttt                                                 <  1:825558/136‑1 (MQ=255)
                                                                                                                   ggtgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGgcg                                          >  2:316544/1‑142 (MQ=255)
                                                                                                                   ggtgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGgcg                                          <  1:316544/142‑1 (MQ=255)
                                                                                                                    gtgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCAc                                  <  3:279612/149‑1 (MQ=255)
                                                                                                                      gcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCt                                <  3:393752/149‑1 (MQ=255)
                                                                                                                          aTCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCaaa                                                            >  4:290339/1‑117 (MQ=255)
                                                                                                                          aTCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCaaa                                                            <  3:290339/117‑1 (MQ=255)
                                                                                                                          aTCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTgg                            >  4:227630/1‑149 (MQ=255)
                                                                                                                           tCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTgg                            >  2:830303/1‑148 (MQ=255)
                                                                                                                            ccAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGc                           >  2:223255/1‑148 (MQ=255)
                                                                                                                              aGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCAtt                        >  3:482254/1‑149 (MQ=255)
                                                                                                                               gatgatGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTg                       >  4:232959/1‑149 (MQ=255)
                                                                                                                                atgatgACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTg                       <  1:223255/148‑1 (MQ=255)
                                                                                                                                   atgaCACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCtt                   <  3:770021/149‑1 (MQ=255)
                                                                                                                                    tgaCACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCtt                   >  2:791521/1‑148 (MQ=255)
                                                                                                                                    tgaCACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTc                  >  3:645564/1‑149 (MQ=255)
                                                                                                                                     gaCACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTc                  >  2:466691/1‑148 (MQ=255)
                                                                                                                                     gaCACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTc                  >  2:202747/1‑148 (MQ=255)
                                                                                                                                     gaCACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTc                  >  1:167922/1‑148 (MQ=255)
                                                                                                                                     gaCACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCa                 >  3:323972/1‑149 (MQ=255)
                                                                                                                                       cacaTTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCAtt               >  4:774827/1‑149 (MQ=255)
                                                                                                                                        acaTTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCAtt               <  1:791521/148‑1 (MQ=255)
                                                                                                                                        acaTTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTg              >  3:354892/1‑149 (MQ=255)
                                                                                                                                         caTTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTg              >  2:631680/1‑148 (MQ=255)
                                                                                                                                         caTTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTg              >  1:152414/1‑148 (MQ=255)
                                                                                                                                          aTTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTgg             <  2:191066/148‑1 (MQ=255)
                                                                                                                                            cccACCTGCCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGAtttt                                                 >  3:280080/2‑110 (MQ=255)
                                                                                                                                            cccACCTGCCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGAtttt                                                 <  4:280080/109‑1 (MQ=255)
                                                                                                                                              cACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACtg         <  2:158693/148‑1 (MQ=255)
                                                                                                                                               aCCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACtgtg       >  3:490967/1‑149 (MQ=255)
                                                                                                                                                ccTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACtgta       <  1:202747/148‑2 (MQ=255)
                                                                                                                                                 cTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACTGTGcc     >  3:673559/1‑149 (MQ=255)
                                                                                                                                                 cTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACTGTGc      >  1:618164/1‑148 (MQ=255)
                                                                                                                                                  tGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGc                                               <  1:152384/106‑1 (MQ=255)
                                                                                                                                                  tGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACTGTGcc     <  2:442354/148‑1 (MQ=255)
                                                                                                                                                   ggCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACTGTGCCtt   >  3:527073/1‑149 (MQ=255)
                                                                                                                                                    gCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACTGTGCCTTa  >  4:57879/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AACAATCCTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACTGTGCCTTA  >  NC_000913/622307‑622603

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: