Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 622,455 G→C 6.2% 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.2% 466.2 / 16.9 130W52C (TGG→TGCentSenterobactin exporter, iron‑regulated
Reads supporting (aligned to +/- strand):  ref base G (58/64);  new base C (4/4);  total (62/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.93e-01

Only 100 of 136 total aligned reads displayed.

AACAATCCTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACTGTGCCTTA  >  NC_000913/622307‑622603
                                                                                                                                                     |                                                                                                                                                    
aaCAATCCTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTgg                                                                                                                                                      <  4:1136715/149‑1 (MQ=255)
    aTCCTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAgg                                                                                                                                                  >  4:1026849/1‑149 (MQ=255)
    aTCCTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAgg                                                                                                                                                  >  2:1366262/1‑149 (MQ=255)
       cTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTgg                                                                                                                                               <  1:106459/149‑1 (MQ=255)
       cTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTgg                                                                                                                                               <  3:1026849/149‑1 (MQ=255)
       cTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTgg                                                                                                                                               <  4:52004/149‑1 (MQ=255)
             tctTAAACTCCACCTGTTGAAAACGCACCCGGCGTTTCGGGCCGTATTCCTCGCTCGGTTCATTTCCATTGTGGCTCCGGGGTTTGCCGGCGTCGCGGGGGCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCttt                                                                                                                                         <  2:904761/148‑1 (MQ=255)
               ttAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTgg                                                                                                                                               >  4:406800/1‑141 (MQ=255)
               ttAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTgg                                                                                                                                               <  3:406800/141‑1 (MQ=255)
               ttAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCg                                                                                                                                       >  2:932662/1‑149 (MQ=255)
               ttAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCg                                                                                                                                       <  4:1437919/149‑1 (MQ=255)
                 aaCCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGt                                                                                                                                     <  3:220403/149‑1 (MQ=255)
                    cTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGAc                                                                                                                                  >  1:366796/1‑149 (MQ=255)
                        gCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTg                                                                                                                              <  1:1366262/149‑1 (MQ=255)
                          cTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGAc                                                                                                                            >  3:427099/1‑149 (MQ=255)
                                 aaaCGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCggcggc                                                                                                                     >  3:1065554/1‑149 (MQ=255)
                                  aaCGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGgcggcg                                                                                                                    <  1:932662/149‑1 (MQ=255)
                                   aCGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGgcgc                                                                                                                   >  2:364440/1‑149 (MQ=255)
                                     gCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGCCAGGTggaat                                                                                                                                            >  4:464869/1‑119 (MQ=255)
                                          cGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGa                                                                                                                 <  4:708923/144‑1 (MQ=255)
                                          cGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGa                                                                                                                 >  3:708923/1‑144 (MQ=255)
                                           ggCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAtgtttg                                                                                                           <  4:1102503/149‑1 (MQ=255)
                                           ggCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAtgtttg                                                                                                           <  2:510752/149‑1 (MQ=255)
                                             cGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAtgtttgtt                                                                                                         <  4:1067388/149‑1 (MQ=255)
                                                    cgcAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGa                                                                                                  <  3:190372/149‑1 (MQ=255)
                                                     gcAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCTGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGAt                                                                                                 >  1:948335/1‑149 (MQ=255)
                                                         tATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTttgttg                                                                                                        <  3:1182744/138‑1 (MQ=255)
                                                              ctcgctcgTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTcggcgg                                                                                        >  3:48249/1‑149 (MQ=255)
                                                              ctcgctcgTTTCATCTCAATTGTGTCTCTGAGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAtgtttg                                                                                                           >  1:1372515/1‑130 (MQ=255)
                                                              ctcgctcgTTTCATCTCAATTGTGTCTCTGAGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAtgtttg                                                                                                           <  2:1372515/130‑1 (MQ=255)
                                                                   tcgTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGc                                                                                   <  4:427099/149‑1 (MQ=255)
                                                                      tttCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTgg                                                                                >  2:809413/1‑149 (MQ=255)
                                                                          aTCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGa                                                                            >  1:1230380/1‑149 (MQ=255)
                                                                           tctcAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGCCAGGTggaa                                                                                                                                             <  4:743800/83‑3 (MQ=255)
                                                                           tctcAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGCCAGGTggaa                                                                                                                                             >  3:743800/1‑81 (MQ=255)
                                                                            ctcAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATc                                                                          >  4:909204/1‑149 (MQ=255)
                                                                               aaTTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCt                                                                       <  1:364440/149‑1 (MQ=255)
                                                                               aaTTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCt                                                                       <  4:681470/149‑1 (MQ=255)
                                                                                 ttGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTAt                                                                     <  2:366796/149‑1 (MQ=255)
                                                                                   gtgtCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAtgtttgtt                                                                                                         >  2:1359526/1‑111 (MQ=255)
                                                                                   gtgtCTCTGGGATTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGcc                                                                                                     >  1:1085966/1‑115 (MQ=255)
                                                                                        tctGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCg                                                                         >  4:20928/1‑138 (MQ=255)
                                                                                        tctGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCg                                                                         <  3:20928/138‑1 (MQ=255)
                                                                                            gggTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTg                                                                                                                                    >  2:232653/1‑76 (MQ=255)
                                                                                              gTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCaaaaaa                                                         <  3:909204/149‑1 (MQ=255)
                                                                                              gTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCaaaaaa                                                         <  3:68511/149‑1 (MQ=255)
                                                                                               tttGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAg                                                        >  4:1035802/1‑149 (MQ=255)
                                                                                                  gCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAtgtttgt                                                                                                          <  1:1154005/96‑1 (MQ=255)
                                                                                                  gCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAtgtttgt                                                                                                          >  2:1154005/1‑96 (MQ=255)
                                                                                                  gCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAtgtttgt                                                                                                          <  3:147245/96‑1 (MQ=255)
                                                                                                   cTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGAt                                                    >  3:457312/1‑149 (MQ=255)
                                                                                                    tCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTgg                                                                                                                                                      <  2:1254349/50‑1 (MQ=255)
                                                                                                    tCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAtgtttgtt                                                                                                         <  3:497367/95‑1 (MQ=255)
                                                                                                    tCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTAt                                                                     <  2:328361/131‑1 (MQ=255)
                                                                                                    tCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTAt                                                                     >  1:328361/1‑131 (MQ=255)
                                                                                                    tCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGAtt                                                   <  2:948335/149‑1 (MQ=255)
                                                                                                    tCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGAtt                                                   >  3:830572/1‑149 (MQ=255)
                                                                                                             cggtgccggtgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGAt                                                                           <  3:958196/116‑1 (MQ=255)
                                                                                                             cggtgccggtgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGAt                                                                           >  4:958196/1‑116 (MQ=255)
                                                                                                              ggtgccggtgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGAt                                                    >  1:915742/1‑138 (MQ=255)
                                                                                                                     gtgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCaaaaaa                                                         >  3:610844/1‑126 (MQ=255)
                                                                                                                     gtgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCaaaaaa                                                         <  4:610844/126‑1 (MQ=255)
                                                                                                                     gtgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGgcgcgc                                       >  4:292614/1‑144 (MQ=255)
                                                                                                                     gtgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCAc                                  <  2:880766/149‑1 (MQ=255)
                                                                                                                        cAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAtgtttgt                                                                                                          <  1:1292069/74‑1 (MQ=255)
                                                                                                                        cAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAtgtttgt                                                                                                          >  2:1292069/1‑74 (MQ=255)
                                                                                                                           aTCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTgg                            <  4:48249/149‑1 (MQ=255)
                                                                                                                               aGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGgt                                                                                       <  2:247070/86‑1 (MQ=255)
                                                                                                                               aGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCaaaaaa                                                         <  4:799609/116‑1 (MQ=255)
                                                                                                                               aGATGATGACACATTCCACCTGGCAGGAGGGGCTTTCGGTGACGTTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGgt                                                                                       >  1:247070/1‑86 (MQ=255)
                                                                                                                                gatgatGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGc                                               <  4:66397/125‑1 (MQ=255)
                                                                                                                                gatgatGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGc                                               >  3:66397/1‑125 (MQ=255)
                                                                                                                                 atgatgACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTgg                      <  1:1282494/149‑1 (MQ=255)
                                                                                                                                   gatgaCACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCt                    >  4:13407/1‑149 (MQ=255)
                                                                                                                                     tgaCACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGa                                                     <  2:233786/114‑1 (MQ=255)
                                                                                                                                     tgaCACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGa                                                     >  1:233786/1‑114 (MQ=255)
                                                                                                                                     tgaCACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTc                  >  4:655268/1‑149 (MQ=255)
                                                                                                                                       acacATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCAt                >  3:164430/1‑149 (MQ=255)
                                                                                                                                       acacATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGACGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCAt                >  1:1010091/1‑149 (MQ=255)
                                                                                                                                         acaTTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTg              >  1:500222/1‑149 (MQ=255)
                                                                                                                                           gccccACCTGCCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGc                           <  4:1355026/131‑1 (MQ=255)
                                                                                                                                           gccccACCTGCCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGc                           >  3:1355026/4‑134 (MQ=255)
                                                                                                                                            ttCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGAc           >  2:1345241/1‑149 (MQ=255)
                                                                                                                                             tCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCAtt                        >  1:1149396/1‑135 (MQ=255)
                                                                                                                                             tCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCAtt                        >  1:1346689/1‑135 (MQ=255)
                                                                                                                                             tCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCAtt                        <  2:1149396/135‑1 (MQ=255)
                                                                                                                                             tCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACt          >  3:972602/1‑149 (MQ=255)
                                                                                                                                             tCCACCTGGCAGATTGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACt          >  3:106583/1‑149 (MQ=255)
                                                                                                                                               cACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCaaaaa                                                          <  4:1168145/99‑1 (MQ=255)
                                                                                                                                               cACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCaaaaa                                                          >  3:1168145/1‑99 (MQ=255)
                                                                                                                                                aCCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGa            >  4:1087584/1‑144 (MQ=255)
                                                                                                                                                aCCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGa            <  3:1087584/144‑1 (MQ=255)
                                                                                                                                                aCCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACtgtg       <  1:1381491/149‑1 (MQ=255)
                                                                                                                                                aCCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTGTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGa            <  1:748043/144‑1 (MQ=255)
                                                                                                                                                aCCTGGCAGGTGGGGCTTTCGGTAACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGa            >  2:748043/1‑144 (MQ=255)
                                                                                                                                                  cTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACTGTGcc     <  4:164430/149‑1 (MQ=255)
                                                                                                                                                   tGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGa                                                     >  4:172104/1‑100 (MQ=255)
                                                                                                                                                    ggCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGa                                                     >  4:669587/1‑99 (MQ=255)
                                                                                                                                                    ggCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGa                                                     <  3:669587/99‑1 (MQ=255)
                                                                                                                                                     gCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACTGTGCCTTa  >  1:613435/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 < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: