Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,546,758 A→G T555A (ACG→GCG)  fimD → type I fimbriae usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,546,7580AG98.5% 249.2 / ‑4.3 68T555A (ACG→GCG) fimDtype I fimbriae usher protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (32/35);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.34e-01

GTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCA  >  NC_000913/4546626‑4546890
                                                                                                                                    |                                                                                                                                    
gTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTaa                                                                                                                               >  1:87025/1‑140 (MQ=255)
     gCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGTAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCg                                                                                                                          >  2:1040676/1‑140 (MQ=255)
     gCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTaa                                                                                                                               <  1:514666/135‑1 (MQ=255)
     gCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTaa                                                                                                                               >  2:514666/1‑135 (MQ=255)
     gCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCg                                                                                                                          >  2:711263/1‑140 (MQ=255)
     gCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCg                                                                                                                          >  2:581329/1‑140 (MQ=255)
     gCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCg                                                                                                                          <  1:561429/140‑1 (MQ=255)
          aaTTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAg                                                                                                                     >  1:180514/1‑140 (MQ=255)
          aaTTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAg                                                                                                                     >  1:554653/1‑140 (MQ=255)
          aaTTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAg                                                                                                                     <  1:277029/140‑1 (MQ=255)
           aTTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGc                                                                                                                    <  1:866475/140‑1 (MQ=255)
               aCCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAAtt                                                                                                                <  2:852140/140‑1 (MQ=255)
               aCCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAAtt                                                                                                                <  2:396235/140‑1 (MQ=255)
                ccGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTc                                                                                                               <  2:180514/140‑1 (MQ=255)
                    cTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAgg                                                                                                           <  2:234163/140‑1 (MQ=255)
                    cTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAgg                                                                                                           <  1:633358/140‑1 (MQ=255)
                     tATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGc                                                                                                          <  2:250518/140‑1 (MQ=255)
                      aTTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCt                                                                                                         <  1:931434/140‑1 (MQ=255)
                      aTTACAACCTCGCTTATAACAAACGCGGGAAACTACAACTTACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTCTTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCt                                                                                                         <  2:937127/140‑1 (MQ=255)
                        tACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTgg                                                                                                       <  1:581329/140‑1 (MQ=255)
                             ccTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTaa                                                                                                  >  2:1021342/1‑140 (MQ=255)
                               tCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCt                                                                                                         <  2:258329/131‑1 (MQ=255)
                               tCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCt                                                                                                         >  1:258329/1‑131 (MQ=255)
                               tCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAAt                                                                                                >  2:49502/1‑140 (MQ=255)
                               tCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAAt                                                                                                >  2:1022462/1‑140 (MQ=255)
                                            aCGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAg                                                                                   <  1:32628/140‑1 (MQ=255)
                                               cgGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGata                                                                                >  1:166137/1‑140 (MQ=255)
                                                 ggAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATc                                                                              >  2:283369/1‑140 (MQ=255)
                                                 ggAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATc                                                                              <  2:166137/140‑1 (MQ=255)
                                                        aCAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGa                                                                       <  1:1021342/140‑1 (MQ=255)
                                                        aCAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGa                                                                       <  1:49502/140‑1 (MQ=255)
                                                            cTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCt                                                                   >  2:1000907/1‑140 (MQ=255)
                                                                      cTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGc                                                         >  2:594583/1‑140 (MQ=255)
                                                                         aGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCt                                                                   >  1:757752/1‑127 (MQ=255)
                                                                         aGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCt                                                                   <  2:757752/127‑1 (MQ=255)
                                                                           cAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCt                                                                   >  2:1043765/1‑125 (MQ=255)
                                                                           cAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCt                                                                   <  1:1043765/125‑1 (MQ=255)
                                                                           cAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGAc                                                    >  2:387516/1‑140 (MQ=255)
                                                                              cTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGaa                                                 >  2:837619/1‑140 (MQ=255)
                                                                              cTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGaa                                                 >  1:407386/1‑140 (MQ=255)
                                                                              cTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGaa                                                 <  2:106765/140‑1 (MQ=255)
                                                                              cTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGaa                                                 >  2:159593/1‑140 (MQ=255)
                                                                                cGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTCGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGaaaa                                               >  1:681042/1‑140 (MQ=255)
                                                                                     gcACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGcc                                          >  2:426792/1‑140 (MQ=255)
                                                                                         aTCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGc                                      <  1:673829/140‑1 (MQ=255)
                                                                                         aTCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGc                                      >  2:54165/1‑140 (MQ=255)
                                                                                         aTCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGc                                      >  2:183470/1‑140 (MQ=255)
                                                                                            aaCACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCaaa                                   <  2:549612/140‑1 (MQ=255)
                                                                                            aaCACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCaaa                                   <  1:223461/140‑1 (MQ=255)
                                                                                                 tGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGa                              >  1:632245/1‑140 (MQ=255)
                                                                                                 tGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCCGACGAAAAACGCCTGGCAAAAAGGa                              <  1:490286/140‑1 (MQ=255)
                                                                                                  gTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGAc                             <  1:837619/140‑1 (MQ=255)
                                                                                                  gTATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGAc                             <  1:426792/140‑1 (MQ=255)
                                                                                                   tATTTGAGTGTTACCCATCACACTTATTGGGGTGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACg                            >  2:1004768/1‑140 (MQ=255)
                                                                                                   tATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACg                            >  2:277950/1‑140 (MQ=255)
                                                                                                   tATTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACg                            >  2:481125/1‑140 (MQ=255)
                                                                                                    aTTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACgg                           <  1:283369/140‑1 (MQ=255)
                                                                                                    aTTTGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACgg                           <  1:277950/140‑1 (MQ=255)
                                                                                                      ttGAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGa                         <  1:594583/140‑1 (MQ=255)
                                                                                                        gAGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATc                       <  1:387516/140‑1 (MQ=255)
                                                                                                         aGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCa                      <  1:323307/140‑1 (MQ=255)
                                                                                                         aGTGGTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCa                      <  1:711263/140‑1 (MQ=255)
                                                                                                             gTAGCCATCAAACTTATTGGGGAGCGAGTAGTGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATg                  >  2:799660/1‑140 (MQ=255)
                                                                                                             gTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATg                  >  1:201889/1‑140 (MQ=255)
                                                                                                             gTAGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATg                  >  1:614698/1‑140 (MQ=255)
                                                                                                               aGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGAc                             <  2:418862/127‑1 (MQ=255)
                                                                                                               aGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGAc                             >  1:418862/1‑127 (MQ=255)
                                                                                                               aGCCATCAAACTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGtt                <  2:201889/140‑1 (MQ=255)
                                                                                                                        aCTTATTGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTaa       <  1:1040676/140‑1 (MQ=255)
                                                                                                                             ttGGGGAGCGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGAGCAGATGTTAGCGCTTAACGTCa  >  1:750327/1‑140 (MQ=255)
                                                                                                                                    |                                                                                                                                    
GTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCA  >  NC_000913/4546626‑4546890

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