Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,138,936 A→T 6.9% T187S (ACA→TCA)  flgK → flagellar hook‑filament junction protein 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,138,9360AT 6.9% 461.7 / 18.1 131T187S (ACA→TCA) flgKflagellar hook‑filament junction protein 1
Reads supporting (aligned to +/- strand):  ref base A (60/61);  new base T (4/5);  total (65/66)
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.45e-01

Only 100 of 133 total aligned reads displayed.

GAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACT  >  NC_000913/1138789‑1139081
                                                                                                                                                   |                                                                                                                                                 
gAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGa                                                                                                                                                   >  4:851283/1‑148 (MQ=255)
        ggTGAATCAGTTTAAAACCACCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTg                                                                                                                                           <  3:851283/148‑1 (MQ=255)
        ggTGAATCAGTTTAAAACCACCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCACGCCTGTCAGGCgcg                                                                                                                                           <  3:1129059/148‑3 (MQ=255)
                 gTTTAAAACCACCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAggg                                                                                                                                  >  1:788821/1‑148 (MQ=255)
                                tCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGt                                                                                                                                            >  1:1011972/1‑123 (MQ=255)
                                tCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGt                                                                                                                                            <  2:1011972/123‑1 (MQ=255)
                                tCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTg                                                                                                                                           <  1:793162/124‑1 (MQ=255)
                                tCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTg                                                                                                                                           >  2:793162/1‑124 (MQ=255)
                                    tatCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTgg                                                                                                                                          >  1:592606/1‑121 (MQ=255)
                                    tatCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTgg                                                                                                                                          <  2:592606/121‑1 (MQ=255)
                                       cTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAgg                                                                                                                                               >  4:1077502/1‑113 (MQ=255)
                                       cTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTgg                                                                                                            >  1:63936/1‑148 (MQ=255)
                                       cTGCGCGACCAGGACAAACAGGTCAATATAGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAgg                                                                                                                                               <  3:1077502/113‑1 (MQ=255)
                                                  ggACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATctg                                                                                                                >  2:904558/1‑133 (MQ=255)
                                                  ggACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATctg                                                                                                                <  1:904558/133‑1 (MQ=255)
                                                    aCAAACAGGTCACTATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAgg                                                                                                                                   <  3:942580/112‑1 (MQ=255)
                                                    aCAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAgg                                                                                                                                   >  4:942580/1‑112 (MQ=255)
                                                    aCAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCa                                                                                               >  4:648104/1‑148 (MQ=255)
                                                         cAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATctgc                                                                                                               >  2:163851/1‑127 (MQ=255)
                                                         cAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATctgc                                                                                                               <  1:163851/127‑1 (MQ=255)
                                                         cAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTgg                                                                                          <  2:788821/148‑1 (MQ=255)
                                                             tCAATATCGCGATAGGTGGCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAg                                                                                      >  2:312416/1‑148 (MQ=255)
                                                             tCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAAcgcg                                                                                                   >  2:145793/1‑135 (MQ=255)
                                                             tCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAAcgcg                                                                                                   <  1:145793/135‑1 (MQ=255)
                                                             tCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAg                                                                                      >  2:1072263/1‑148 (MQ=255)
                                                              cAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATc                                                                                                                  >  4:315271/1‑119 (MQ=255)
                                                              cAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATc                                                                                                                  <  3:315271/119‑1 (MQ=255)
                                                                  atCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAAt                                                                                 >  2:952206/1‑148 (MQ=255)
                                                                        aTAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAAcc                                                                           <  1:312416/148‑1 (MQ=255)
                                                                        aTAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAAcc                                                                           <  1:1072263/148‑1 (MQ=255)
                                                                          aGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAg                                                                         <  3:648104/148‑1 (MQ=255)
                                                                           ggTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGa                                                                        <  2:63936/148‑1 (MQ=255)
                                                                             tGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTg                                                                                                                                           >  3:240313/1‑79 (MQ=255)
                                                                             tGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTg                                                                                                                                           <  4:240313/79‑1 (MQ=255)
                                                                              gCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATctgct                                                                                                              >  2:563637/1‑107 (MQ=255)
                                                                              gCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATctgct                                                                                                              <  1:563637/107‑1 (MQ=255)
                                                                                cAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTa                                                                                                                        <  1:575906/95‑1 (MQ=255)
                                                                                cAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTa                                                                                                                        >  2:575906/1‑95 (MQ=255)
                                                                                  gCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAgg                                                                                                                                   >  2:339789/1‑82 (MQ=255)
                                                                                  gCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAgg                                                                                                                                   <  1:339789/82‑1 (MQ=255)
                                                                                  gCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCa                                                                                                                            >  3:1058499/1‑89 (MQ=255)
                                                                                     ttGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGAttgtt                                                                   >  4:134906/1‑143 (MQ=255)
                                                                                         tcagatcaACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTACCAATct                                                                                                                 <  2:548863/93‑1 (MQ=255)
                                                                                         tcagatcaACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATct                                                                                                                 >  1:548863/1‑93 (MQ=255)
                                                                                          cagatcaACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGa                                                                        >  3:326098/1‑133 (MQ=255)
                                                                                            gatcaACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAAt                                                                                                                   >  3:143282/1‑88 (MQ=255)
                                                                                            gatcaACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAAt                                                                                                                   <  4:143282/88‑1 (MQ=255)
                                                                                             atcaACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAg                                                                                                                                    >  3:22111/1‑70 (MQ=255)
                                                                                             atcaACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGa                                                                                   >  4:831544/1‑119 (MQ=255)
                                                                                             atcaACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGa                                                           <  3:522093/143‑1 (MQ=255)
                                                                                             atcaACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGa                                                           >  4:522093/1‑143 (MQ=255)
                                                                                              tcaACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAgg                                                                                                                                   <  2:978728/70‑1 (MQ=255)
                                                                                              tcaACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAgg                                                                                                                                   >  1:978728/1‑70 (MQ=255)
                                                                                                 acaacTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCg                                                                                    <  3:576277/114‑1 (MQ=255)
                                                                                                 acaacTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAAATGGTGAGCg                                                                                    >  4:576277/1‑114 (MQ=255)
                                                                                                  caacTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGtt                                                 <  1:348048/148‑1 (MQ=255)
                                                                                                  caacTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGtt                                                 >  2:348048/1‑148 (MQ=255)
                                                                                                   aacTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGAc                                                                                                                                                  >  1:1075358/1‑50 (MQ=255)
                                                                                                   aacTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGAc                                                                                                                                                  <  2:1075358/50‑1 (MQ=255)
                                                                                                     cTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGAttgt                                                                    <  2:643604/126‑1 (MQ=255)
                                                                                                     cTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGAttgt                                                                    >  1:643604/1‑126 (MQ=255)
                                                                                                          cTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTaa                                                                              <  1:1005953/111‑1 (MQ=255)
                                                                                                          cTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGAttg                                                                     >  4:590935/1‑120 (MQ=255)
                                                                                                            aaaCAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATggcg                                       >  2:892801/1‑148 (MQ=255)
                                                                                                              aCAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGa                                                           <  1:1035935/126‑1 (MQ=255)
                                                                                                              aCAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGa                                                           >  2:1035935/1‑126 (MQ=255)
                                                                                                               cAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAg                                                     <  1:1082156/131‑1 (MQ=255)
                                                                                                               cAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAg                                                     >  2:1082156/1‑131 (MQ=255)
                                                                                                                 aaTTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGAt                                                                                                 >  4:310621/1‑85 (MQ=255)
                                                                                                                 aaTTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACt                                  >  3:899127/1‑148 (MQ=255)
                                                                                                                   ttGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTta                                >  4:341819/1‑148 (MQ=255)
                                                                                                                   ttGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTAAGCGTTCAGGATGGCAGCACTta                                >  2:102894/1‑148 (MQ=255)
                                                                                                                    tGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTtat                               <  1:952206/148‑1 (MQ=255)
                                                                                                                         gCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAAcgcg                                                                                                   >  2:249659/1‑75 (MQ=255)
                                                                                                                         gCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAAcgcg                                                                                                   <  1:249659/75‑1 (MQ=255)
                                                                                                                             gAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACg                      <  3:839869/148‑1 (MQ=255)
                                                                                                                               aCGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGAttgtt                                                                   <  2:877662/101‑1 (MQ=255)
                                                                                                                               aCGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGAttgtt                                                                   >  1:877662/1‑101 (MQ=255)
                                                                                                                               aCGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAg                                                     >  2:266889/1‑115 (MQ=255)
                                                                                                                               aCGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAg                                                     <  1:266889/115‑1 (MQ=255)
                                                                                                                                  aTCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATctgct                                                                                                              <  3:943460/55‑1 (MQ=255)
                                                                                                                                  aTCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGaa                                                                                  >  4:48502/1‑83 (MQ=255)
                                                                                                                                  aTCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGaa                                                                                  <  3:48502/83‑1 (MQ=255)
                                                                                                                                    cAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTg                                                                                           >  1:799010/1‑72 (MQ=255)
                                                                                                                                    cAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTg                                                                                           <  2:799010/72‑1 (MQ=255)
                                                                                                                                    cAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGa                                                           >  1:40266/1‑104 (MQ=255)
                                                                                                                                    cAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGa                                                           <  2:40266/104‑1 (MQ=255)
                                                                                                                                     aaaTCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCaa              <  4:1021098/148‑1 (MQ=255)
                                                                                                                                        tctcGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATgg           <  1:892801/148‑1 (MQ=255)
                                                                                                                                           cgcgCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTa        <  1:102894/148‑1 (MQ=255)
                                                                                                                                           cgcgCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTa        <  3:341819/148‑1 (MQ=255)
                                                                                                                                           cacgcCTGTCAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGAttgt                                                                    <  4:362960/86‑1 (MQ=255)
                                                                                                                                           cacgcCTGTCAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGAt                                                                       <  2:478858/83‑1 (MQ=255)
                                                                                                                                           cacgcCTGTCAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGAt                                                                       >  1:478858/3‑85 (MQ=255)
                                                                                                                                           cacgcCTGTCAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGAt                                                                       <  1:98078/83‑1 (MQ=255)
                                                                                                                                           cacgcCTGTCAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGAt                                                                       >  2:98078/3‑85 (MQ=255)
                                                                                                                                            gcgcCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTAc       <  4:899127/148‑1 (MQ=255)
                                                                                                                                               ccTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGa                                                                                                  <  4:741016/54‑1 (MQ=255)
                                                                                                                                               ccTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGa                                                                                                  >  3:741016/1‑54 (MQ=255)
                                                                                                                                               ccTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATg                   >  3:637043/1‑133 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACT  >  NC_000913/1138789‑1139081

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