Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 623,054 A→T 11.4% N252I (AAC→ATC)  entS → enterobactin exporter, iron‑regulated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913623,0540AT11.4% 218.2 / 18.5 70N252I (AAC→ATC) entSenterobactin exporter, iron‑regulated
Reads supporting (aligned to +/- strand):  ref base A (33/29);  new base T (4/4);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.01e-01

CAGCCGCGTGAGCATCCGTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCTCTCCACGCTGGGATCGTTCCTC  >  NC_000913/622909‑623199
                                                                                                                                                 |                                                                                                                                                 
cAGCCGCGTGAGCATCCGTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACt                                                                                                                                                 >  1:687968/1‑148 (MQ=255)
   ccGCGTGAGCATCCGTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACATCTgcc                                                                                                                                              >  3:689717/1‑146 (MQ=255)
   ccGCGTGAGCATCCGTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCATGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACATCTgcc                                                                                                                                              >  4:77637/1‑146 (MQ=255)
          aGCATCCGTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACATCTGCCAGTTGTc                                                                                                                                       <  4:689717/148‑1 (MQ=255)
          aGCATCCGTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCATGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACCTCTGCCAGTTGTc                                                                                                                                       <  3:77637/148‑1 (MQ=255)
                  ttGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACaa                                                                                                                                                   >  1:11308/1‑128 (MQ=255)
                  ttGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACaa                                                                                                                                                   <  2:11308/128‑1 (MQ=255)
                   tGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGc                                                                                                                                              >  2:660681/1‑132 (MQ=255)
                   tGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGc                                                                                                                                              <  3:633690/132‑1 (MQ=255)
                   tGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGc                                                                                                                                              >  4:633690/1‑132 (MQ=255)
                   tGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGc                                                                                                                                              <  1:660681/132‑1 (MQ=255)
                      aaTCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGAt                                                                                                                           >  2:435827/1‑148 (MQ=255)
                        tCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGt                                                                                                                                        >  2:613755/1‑133 (MQ=255)
                        tCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGt                                                                                                                                        <  1:613755/133‑1 (MQ=255)
                              cTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGc                                                                                                                                              >  1:1099326/1‑121 (MQ=255)
                              cTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGc                                                                                                                                              <  2:1099326/121‑1 (MQ=255)
                                  cAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTa                                                                                                               <  2:364882/148‑1 (MQ=255)
                                        ttCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTtct                                                                                                                  <  2:1009868/139‑1 (MQ=255)
                                        ttCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTtct                                                                                                                  >  3:677595/1‑139 (MQ=255)
                                        ttCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTtct                                                                                                                  <  4:677595/139‑1 (MQ=255)
                                        ttCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTtct                                                                                                                  >  1:1009868/1‑139 (MQ=255)
                                             tttCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGt                                                                                                                       >  1:567896/1‑129 (MQ=255)
                                             tttCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGt                                                                                                                       <  2:567896/129‑1 (MQ=255)
                                                 tGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCt                                                                                                <  1:435827/148‑1 (MQ=255)
                                                               cTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTg                                                                                  >  4:818752/1‑148 (MQ=255)
                                                                       cGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGtttt                                                                                                                    <  4:158668/106‑1 (MQ=255)
                                                                       cGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGtttt                                                                                                                    >  3:158668/1‑106 (MQ=255)
                                                                         ggATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAAc                                                                        >  1:806673/1‑148 (MQ=255)
                                                                            ttGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACgcggcg                                                                                                        >  4:173755/1‑113 (MQ=255)
                                                                            ttGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACgcggcg                                                                                                        <  3:173755/113‑1 (MQ=255)
                                                                            ttGCGCTGCTGGGTGGTTTATGGACGATGGCGAGCGCGGTGCGGGTACTGAATCCGGCGCTGGCTGACAAGTGGCAGATGTCAGCGGCACAGATTGGTTTGCTCTACGAGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAGCCAg                                                                     >  4:710207/1‑148 (MQ=255)
                                                                                 ctgctgGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTgg                                                                                                                        <  4:756125/92‑1 (MQ=255)
                                                                                 cagctgGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTgg                                                                                                                        >  3:756125/3‑92 (MQ=255)
                                                                                      gggTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTg                                                           <  3:818752/148‑1 (MQ=255)
                                                                                        gtggtTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGc                                                         >  3:75748/1‑148 (MQ=255)
                                                                                            tttATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAAcc                                                                       >  3:593104/1‑130 (MQ=255)
                                                                                            tttATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAAcc                                                                       <  4:593104/130‑1 (MQ=255)
                                                                                             ttATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACgcgg                                                                                                          <  4:49442/94‑1 (MQ=255)
                                                                                             ttATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACgcgg                                                                                                          >  3:49442/1‑94 (MQ=255)
                                                                                               aTTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGt                                                  >  1:282640/1‑148 (MQ=255)
                                                                                                 tGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGa                                                                <  3:477353/132‑1 (MQ=255)
                                                                                                 tGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGa                                                                >  4:477353/1‑132 (MQ=255)
                                                                                                  gACGATGGGGAGCGCGGTGCGGGTACTGTATCAGGCGCTGACTGACAACTGGCAGATGTCATCGGCACAGATTGGTTTTCTCTACGCGGCAATCCCGCTCGGCGCGGCTATTCGTGCGTTAACCAGCGGGAAGCTGGCACATAGTgcg                                               <  3:710207/148‑1 (MQ=255)
                                                                                                    cGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACATCTGCCAGTTGTCAg                                                                                                                                     >  3:265555/1‑60 (MQ=37)
                                                                                                    cGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACATCTGCCAGTTGTCAg                                                                                                                                     <  4:265555/60‑1 (MQ=37)
                                                                                                    cGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGttt                                                                                                                     >  3:566250/1‑76 (MQ=255)
                                                                                                    cGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGttt                                                                                                                     <  4:566250/76‑1 (MQ=255)
                                                                                                        ggCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCa                                         <  2:806673/148‑1 (MQ=255)
                                                                                                         gCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGTCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAg                                        <  4:75748/148‑1 (MQ=255)
                                                                                                            aGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGcgcg                                                                                         >  4:757423/1‑96 (MQ=255)
                                                                                                            aGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGcgcg                                                                                         <  3:757423/96‑1 (MQ=255)
                                                                                                                   tGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGgcg                                                                                           <  1:543398/87‑1 (MQ=255)
                                                                                                                   tGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGgcg                                                                                           >  2:543398/1‑87 (MQ=255)
                                                                                                                         tACTGTATCCGGCGCTGACTGACAACTGGCAGATGTCAGCGGCACCGATTGGt                                                                                                                       >  2:954856/1‑53 (MQ=255)
                                                                                                                         tACTGTATCCGGCGCTGACTGACAACTGGCAGATGTCAGCGGCACAGATTGGt                                                                                                                       <  1:954856/53‑1 (MQ=255)
                                                                                                                             gTATCTGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCTCTcc                    >  3:723653/1‑148 (MQ=255)
                                                                                                                               aTCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACAt                                                     >  1:659456/1‑113 (MQ=255)
                                                                                                                               aTCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACAt                                                     >  4:178007/1‑113 (MQ=255)
                                                                                                                               aTCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACAt                                                     <  2:659456/113‑1 (MQ=255)
                                                                                                                               aTCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACAt                                                     <  3:178007/113‑1 (MQ=255)
                                                                                                                               aTCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGt                                                  >  4:1074491/1‑116 (MQ=255)
                                                                                                                               aTCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGt                                                  <  3:1074491/116‑1 (MQ=255)
                                                                                                                                  cGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCTCTCCACGCt               <  2:282640/148‑1 (MQ=255)
                                                                                                                                       ctggctgACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCTCTCCACGCTGGGAt          >  3:967237/1‑148 (MQ=255)
                                                                                                                                       ctggctgACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCTCTCCACGCTGGGAt          >  3:1010670/1‑148 (MQ=255)
                                                                                                                                            tgACATCTGCCAGTTGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCTCTCCACGCTGGGATCGTTc     >  1:44354/1‑148 (MQ=255)
                                                                                                                                              aCATCTGCCAGTTGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCTCTCCACGCTGGGATCGTTCCt   <  2:44354/148‑1 (MQ=255)
                                                                                                                                              aCAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGaa                                                               <  4:620141/88‑1 (MQ=255)
                                                                                                                                              aCAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGaa                                                               >  3:620141/1‑88 (MQ=255)
                                                                                                                                               cAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTgg                                                                                 >  4:529987/1‑69 (MQ=255)
                                                                                                                                               cAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTgg                                                                                 <  3:529987/69‑1 (MQ=255)
                                                                                                                                               cAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCTCTCCACGCTGGGATCGTTCCTc  >  1:138261/1‑148 (MQ=255)
                                                                                                                                                 aCTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGt                                     <  4:367153/111‑1 (MQ=255)
                                                                                                                                                 aCTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGt                                     >  3:367153/1‑111 (MQ=255)
                                                                                                                                                 aCTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCTCTCCACGCTGGGATCGtt      <  1:480622/142‑1 (MQ=255)
                                                                                                                                                 aCTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCTCTCCACGCTGGGATCGtt      >  2:480622/1‑142 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CAGCCGCGTGAGCATCCGTTGAAATCATTACTGGCAGGATTTCGTTTTCTGCTTGCCAGCCCGCTGGTGGGCGGGATTGCGCTGCTGGGTGGTTTATTGACGATGGCGAGCGCGGTGCGGGTACTGTATCCGGCGCTGGCTGACAACTGGCAGATGTCAGCGGCACAGATTGGTTTTCTCTACGCGGCGATCCCGCTCGGCGCGGCTATTGGTGCGTTAACCAGCGGGAAGCTGGCACATAGTGCGCGACCAGGGTTATTGATGCTGCTCTCCACGCTGGGATCGTTCCTC  >  NC_000913/622909‑623199

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: