Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,018,986 A→T 8.5% intergenic (+49/+27) matP → / ← ompA Ter macrodomain organizer matS‑binding protein/outer membrane protein A (3a;II*;G;d)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,018,9860AT 8.5% 465.4 / 25.3 139intergenic (+49/+27)matP/ompATer macrodomain organizer matS‑binding protein/outer membrane protein A (3a;II*;G;d)
Reads supporting (aligned to +/- strand):  ref base A (70/57);  new base T (7/5);  total (77/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

Only 100 of 146 total aligned reads displayed.

GATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCA  >  NC_000913/1018838‑1019133
                                                                                                                                                     |                                                                                                                                                   
gATTGTTGAGCTGATTGAAGAAGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCgcagca                                                                                                                                                     >  2:85696/1‑149 (MQ=255)
   tGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCg                                                                                                                                                  <  2:961473/149‑1 (MQ=255)
        aGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCTGCTGCGGGGtt                                                                                                                                             <  1:326125/149‑1 (MQ=255)
        aGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCTGCTGCGGGGtt                                                                                                                                             <  3:991069/149‑1 (MQ=255)
         gCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGA‑AAAAAAAACCCCGCAGCAGCGGGGttt                                                                                                                                            <  3:154931/149‑1 (MQ=255)
         gCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGCTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGttt                                                                                                                                            >  4:640447/1‑149 (MQ=255)
         gCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGAAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGttt                                                                                                                                            >  3:1034710/1‑149 (MQ=255)
          cTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCTGCTGCGGGGtttt                                                                                                                                           <  4:587587/149‑1 (MQ=255)
           tGATTGAAGAGGCGGAAAACAGAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTGAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGttttt                                                                                                                                          <  3:231993/149‑1 (MQ=255)
               tGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTAc                                                                                                                                      <  1:197294/149‑1 (MQ=255)
                 aaGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCa                                                                                                                                    >  4:949566/1‑149 (MQ=255)
                 aaGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCa                                                                                                                                    >  2:716959/1‑149 (MQ=255)
                      gCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGCCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACga                                                                                                                               <  3:624549/149‑1 (MQ=255)
                         gAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGaa                                                                                                                            <  1:779368/149‑1 (MQ=255)
                         gAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGaa                                                                                                                            >  4:944669/1‑149 (MQ=255)
                           aaaCAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACt                                                                                                                          <  3:277053/149‑1 (MQ=255)
                            aaCAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCTGCTGCGGGGtt                                                                                                                                             >  3:897105/1‑129 (MQ=255)
                            aaCAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCTGCTGCGGGGtt                                                                                                                                             <  4:897105/129‑1 (MQ=255)
                              cAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTaa                                                                                                                       <  1:181951/149‑1 (MQ=255)
                                aaGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGc                                                                                                                     <  1:178407/149‑1 (MQ=255)
                                 agagAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTc                                                                                                                                         >  1:426658/1‑128 (MQ=255)
                                     aaaTACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCg                                                                                                                <  3:52275/149‑1 (MQ=255)
                                         aCGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGCTTTTCTACCAGACGAGCACTTAAGCCTGCGGCTg                                                                                                            >  3:416894/1‑149 (MQ=255)
                                          cgcgAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGa                                                                                                           <  1:716959/149‑1 (MQ=255)
                                          cgcgAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGa                                                                                                           <  1:151445/149‑1 (MQ=255)
                                           gcgAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAg                                                                                                          >  3:1076852/1‑149 (MQ=255)
                                              aaTAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTa                                                                                                       <  2:1097814/149‑1 (MQ=255)
                                              aaTAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTa                                                                                                       >  1:229153/1‑149 (MQ=255)
                                               aTAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTAc                                                                                                      <  3:334256/149‑1 (MQ=255)
                                                 aaaaTGTCTTCTATGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAATTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACaa                                                                                                    <  4:927643/149‑1 (MQ=255)
                                                  aaaTGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACaa                                                                                                    <  2:152848/149‑1 (MQ=255)
                                                   aaTGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACg                                                                                                  <  3:223326/149‑1 (MQ=255)
                                                     tGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTc                                                                                                <  4:1076852/149‑1 (MQ=255)
                                                        cttctCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGGATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCttt                                                                                             <  1:522619/149‑1 (MQ=255)
                                                        cttctCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCttt                                                                                             <  1:554464/149‑1 (MQ=255)
                                                         ttctCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTg                                                                                            >  4:876876/1‑149 (MQ=255)
                                                         ttctCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTg                                                                                            >  3:278147/1‑149 (MQ=255)
                                                         ttctCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTg                                                                                            >  4:169162/1‑149 (MQ=255)
                                                            tctGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCgg                                                                                                               <  2:745143/127‑1 (MQ=255)
                                                            tctGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCgg                                                                                                               >  1:745143/1‑127 (MQ=255)
                                                            tctGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATa                                                                                         >  1:108139/1‑149 (MQ=255)
                                                             ctGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATAc                                                                                        <  3:576964/149‑1 (MQ=255)
                                                               gAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGttttt                                                                                                                                          >  1:350800/1‑97 (MQ=255)
                                                               gAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGttttt                                                                                                                                          <  2:350800/97‑1 (MQ=255)
                                                                   cAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGCCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTaa                                                                                  <  1:85696/149‑1 (MQ=255)
                                                                     ggATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACt                                                                                >  4:208350/1‑149 (MQ=255)
                                                                         cTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCt                                                                                                                   >  3:646522/1‑110 (MQ=255)
                                                                          tGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGAt                                                                           <  2:386292/149‑1 (MQ=255)
                                                                               gCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTa                                                                      >  3:1041580/1‑149 (MQ=255)
                                                                                              gAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGGTACCTTTAACTTCGATCTCTACGCGACGATCCGGAg                                                       >  3:352891/1‑149 (MQ=255)
                                                                                              gAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAg                                                       >  3:584656/1‑149 (MQ=255)
                                                                                              gAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGCTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAg                                                       >  4:809153/1‑149 (MQ=255)
                                                                                               aaTAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTAc                                                                     >  4:240498/1‑134 (MQ=255)
                                                                                               aaTAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTAc                                                                     <  3:240498/134‑1 (MQ=255)
                                                                                                aTAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGcc                                                     >  3:886488/1‑149 (MQ=255)
                                                                                                aTAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGcc                                                     >  2:276533/1‑149 (MQ=255)
                                                                                                aTAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGcc                                                     >  1:893073/1‑149 (MQ=255)
                                                                                                 tAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCa                                                    <  4:685677/149‑1 (MQ=255)
                                                                                                    cTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGc                                                 >  4:1037793/1‑149 (MQ=255)
                                                                                                     tGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCa                                                <  3:547399/149‑1 (MQ=255)
                                                                                                       aCGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGt                                              >  4:270352/1‑149 (MQ=255)
                                                                                                         gAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCg                                            >  3:416162/1‑149 (MQ=255)
                                                                                                          aaaGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGa                                           <  1:780287/149‑1 (MQ=255)
                                                                                                           aaGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGAt                                          <  3:934749/149‑1 (MQ=255)
                                                                                                             gTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCa                                        >  4:460236/1‑149 (MQ=255)
                                                                                                                  ttCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTc                                             <  4:132212/139‑1 (MQ=255)
                                                                                                                  ttCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTc                                             >  3:132212/1‑139 (MQ=255)
                                                                                                                  ttCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTgca                                   >  2:772167/1‑149 (MQ=255)
                                                                                                                  ttCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTgca                                   >  2:223759/1‑149 (MQ=255)
                                                                                                                       tttACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGTTTTTTCTACGAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTCTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACg                              >  2:197455/1‑149 (MQ=255)
                                                                                                                       tttACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACg                              >  4:601318/1‑149 (MQ=255)
                                                                                                                        ttACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGc                             <  4:827494/149‑1 (MQ=255)
                                                                                                                        ttACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGCTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTAGATACCTTTAACTTCGATCTCTACGCGACGATCCGGATCCAGGCAGTCGATCAGTGCAGCACGc                             >  1:122119/1‑149 (MQ=255)
                                                                                                                            tAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATcc                                                           >  3:940156/1‑115 (MQ=255)
                                                                                                                            tAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGtt                         >  2:1029421/1‑149 (MQ=255)
                                                                                                                            tAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGtt                         >  4:455017/1‑149 (MQ=255)
                                                                                                                             aaaGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGttt                        >  2:898954/1‑149 (MQ=255)
                                                                                                                                 gC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGt                                      >  3:1012371/1‑131 (MQ=255)
                                                                                                                                 gC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGt                                      <  4:1012371/131‑1 (MQ=255)
                                                                                                                                  c‑AAAAAAAACCCTTCAGCAGCGGGGTTTTTCAACCAGACGAGAACTTAAGCCTGCTGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCACCACGCTGTTTCACGt                   <  2:851471/149‑1 (MQ=255)
                                                                                                                                    aaaaaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGtt                  >  1:729700/1‑149 (MQ=255)
                                                                                                                                     aaaaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTg                 >  2:767108/1‑149 (MQ=255)
                                                                                                                                      gaaaaaCCCCGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGt                >  2:815104/2‑149 (MQ=255)
                                                                                                                                      gaaaaaCCCCGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCAGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTg                                     >  3:207671/2‑128 (MQ=255)
                                                                                                                                      aaaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGt                <  3:455017/149‑1 (MQ=255)
                                                                                                                                      aaaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGt                >  2:591676/1‑149 (MQ=255)
                                                                                                                                      aaaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGt                <  1:591676/149‑1 (MQ=255)
                                                                                                                                         aaaCCCCGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTcac             >  4:914182/1‑149 (MQ=255)
                                                                                                                                         aaaCCCCGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTcac             >  4:36026/1‑149 (MQ=255)
                                                                                                                                          aaCCCCGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTcaca            >  2:909129/1‑149 (MQ=255)
                                                                                                                                          aaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGTGGCTGAGATACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGCCGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTcaca            >  4:699488/1‑149 (MQ=255)
                                                                                                                                           aCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAg           <  3:949566/149‑1 (MQ=255)
                                                                                                                                             cccGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGgt         >  1:666469/1‑149 (MQ=255)
                                                                                                                                             cccGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGgt         >  2:411133/1‑149 (MQ=255)
                                                                                                                                              ccGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGgtg        >  1:852438/1‑149 (MQ=255)
                                                                                                                                                gcagcaGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGtt      <  1:411133/149‑1 (MQ=255)
                                                                                                                                                gcagcaGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGtt      <  1:898954/149‑1 (MQ=255)
                                                                                                                                                  agcagcGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGc    <  1:551852/149‑1 (MQ=255)
                                                                                                                                                  agcagcGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGc    <  2:47846/149‑1 (MQ=255)
                                                                                                                                                    cagcGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCa  <  1:619827/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCA  >  NC_000913/1018838‑1019133

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