Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,080,009 G→A 6.7% intergenic (+40/+290) yiiF → / ← fdhE putative thymol sensitivity protein, CopG family putative transcriptional regulator/formate dehydrogenase formation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,080,0090GA 6.7% 392.5 / 12.4 119intergenic (+40/+290)yiiF/fdhEputative thymol sensitivity protein, CopG family putative transcriptional regulator/formate dehydrogenase formation protein
Reads supporting (aligned to +/- strand):  ref base G (53/58);  new base A (4/4);  total (57/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.99e-01

Only 100 of 122 total aligned reads displayed.

TACCTGATAAATCAATCGCAAACAGCAAGAACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTAC  >  NC_000913/4079862‑4080125
                                                                                                                                                   |                                                                                                                    
tACCTGATAAATCAATCGCAAACATCAAGAACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTgcgc                                                                                                                     >  3:1327136/1‑149 (MQ=255)
   cTGATAAATCAATCGCAAACAGCAAGAACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCgtt                                                                                                                  <  4:197374/149‑1 (MQ=255)
     gATAAATCAATCGCAAACAGCAAGAACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGcgcg                                                                                                                    <  3:1370211/145‑1 (MQ=255)
     gATAAATCAATCGCAAACAGCAAGAACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCgttgt                                                                                                                >  3:1414738/1‑149 (MQ=255)
          atcaatcGCAAACAGCAAGAACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATg                                                                                                           >  1:1169162/1‑149 (MQ=255)
           tcaatcGCAAACAGCAAGAACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGc                                                                                                          >  2:489931/1‑149 (MQ=255)
               tcGCAAACAGCAAGAACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGa                                                                                                      >  1:107056/1‑149 (MQ=255)
                   aaaCAGCAAGAACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGAtaag                                                                                                  <  2:107056/149‑4 (MQ=255)
                     aCAGCAAGAACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATgcggc                                                                                                <  4:1414738/149‑1 (MQ=255)
                       aGCAAGAACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGt                                                                                              <  2:1169162/149‑1 (MQ=255)
                          aaGAACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATAAGGCGttca                                                                                           >  2:761148/1‑146 (MQ=255)
                          aaGAACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAAATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATAAGGCGttca                                                                                           >  4:853521/1‑146 (MQ=255)
                           aGAACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGAt                                                                                                     <  2:564139/138‑1 (MQ=255)
                           aGAACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGAt                                                                                                     >  1:564139/1‑138 (MQ=255)
                           aGAACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAAc                                                                                          >  2:1274214/1‑149 (MQ=255)
                             aacactggtTCTTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCATTGTTCATGCCGGATGCGGCGTGAACGc                                                                                        <  3:95624/142‑1 (MQ=255)
                             aaCCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGc                                                                                        <  1:1274214/149‑1 (MQ=255)
                              aCCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                                       >  4:806243/1‑149 (MQ=255)
                                cAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGc                                                                                        >  2:870503/1‑146 (MQ=255)
                                cAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGc                                                                                        <  1:870503/146‑1 (MQ=255)
                                 aGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGt                                                                                              <  1:1373886/139‑1 (MQ=255)
                                  gtgtTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATgcggc                                                                                                <  2:182394/136‑1 (MQ=255)
                                  gtgtTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATgcg                                                                                                  >  1:893157/1‑134 (MQ=255)
                                    gtTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATAAGGCGTTCACGCCGCATcc                                                                                 <  3:853521/149‑1 (MQ=38)
                                       ccTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                              <  2:478099/149‑1 (MQ=255)
                                       ccTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                              <  3:1361138/149‑1 (MQ=255)
                                        cTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                             >  3:665886/1‑149 (MQ=255)
                                         tGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATAAGGCGTTCACGccgcatc                                                                                  >  2:1464103/1‑138 (MQ=38)
                                              tCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaa                                                                       <  4:1327136/149‑1 (MQ=255)
                                              tCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaa                                                                       <  3:806243/149‑1 (MQ=255)
                                                 ggCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTc                                                                    >  1:335664/1‑149 (MQ=255)
                                                  gCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCgtt                                                                                                                  >  1:757670/1‑102 (MQ=255)
                                                  gCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCgtt                                                                                                                  <  2:757670/102‑1 (MQ=255)
                                                   cAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaa                                                                       <  1:600273/144‑1 (MQ=255)
                                                    aaGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGta                                                                 <  3:298606/149‑1 (MQ=255)
                                                         tgtgAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGta                                                                 <  2:123062/144‑1 (MQ=255)
                                                         tgtgAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGta                                                                 >  1:123062/1‑144 (MQ=255)
                                                         tgtgAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAAt                                                            <  4:665886/149‑1 (MQ=255)
                                                           tgAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                              <  4:1246103/129‑1 (MQ=255)
                                                           tgAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                              >  3:1246103/1‑129 (MQ=255)
                                                             aGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCaa                                                        >  3:170627/1‑149 (MQ=255)
                                                              ggAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                <  4:582056/124‑1 (MQ=255)
                                                              ggAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                >  3:582056/1‑124 (MQ=255)
                                                              ggAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAAt                                                       <  2:468553/149‑1 (MQ=255)
                                                                 aGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCgttgt                                                                                                                >  1:1016109/1‑89 (MQ=255)
                                                                 aGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCgttgt                                                                                                                <  2:1016109/89‑1 (MQ=255)
                                                                    tGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTg                                                 <  4:35462/149‑1 (MQ=255)
                                                                      gtgtCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCa                                               <  3:1120815/149‑1 (MQ=255)
                                                                        gtCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGt                                             >  4:249028/1‑149 (MQ=255)
                                                                         tCGAATATCAGCGTGAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAAt                                                       >  3:106721/1‑138 (MQ=255)
                                                                         tCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAAt                                                       <  4:106721/138‑1 (MQ=255)
                                                                         tCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCa                                               >  1:157309/1‑146 (MQ=255)
                                                                         tCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCa                                               <  2:157309/146‑1 (MQ=255)
                                                                            aaTATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTAt                                         >  1:240257/1‑149 (MQ=255)
                                                                            aaTATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGCGCACgttgtca                                                                                                              >  3:1424750/1‑78 (MQ=255)
                                                                            aaTATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGCGCACgttgtca                                                                                                              <  4:1424750/80‑3 (MQ=255)
                                                                               atCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATaaa                                                             <  4:157519/126‑1 (MQ=255)
                                                                                tCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGt                                     >  4:1348125/1‑149 (MQ=255)
                                                                                  aGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAg                                   <  2:335664/149‑1 (MQ=255)
                                                                                  aGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCACGCCGCATCCGGCATg                                                                                             <  1:225236/91‑1 (MQ=37)
                                                                                  aGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCACGCCGCATCCGGCATg                                                                                             >  2:225236/1‑91 (MQ=37)
                                                                                         gCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTg                                       >  1:1064761/1‑138 (MQ=255)
                                                                                         gCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTg                                       <  2:1064761/138‑1 (MQ=255)
                                                                                          cTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGt                                                                     <  3:1017192/107‑1 (MQ=255)
                                                                                          cTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGt                                                                     >  4:1017192/1‑107 (MQ=255)
                                                                                              gcgAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTa                                                    >  3:872428/1‑120 (MQ=255)
                                                                                              gcgAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATGGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTa                                                    <  4:872428/120‑1 (MQ=255)
                                                                                                  ggAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAg                                   >  4:707302/1‑133 (MQ=255)
                                                                                                  ggAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAg                                   <  3:707302/133‑1 (MQ=255)
                                                                                                    aaTGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAagcgcagc                   <  4:1051861/147‑1 (MQ=255)
                                                                                                      tGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCa                                               <  3:1464722/117‑1 (MQ=255)
                                                                                                      tGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCa                                               >  4:1464722/1‑117 (MQ=255)
                                                                                                         tAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCaa                                                        <  3:224409/105‑1 (MQ=255)
                                                                                                         tAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCaa                                                        >  4:224409/1‑105 (MQ=255)
                                                                                                         tAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGCGCACgttgtca                                                                                                              >  3:159964/1‑49 (MQ=37)
                                                                                                         tAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGCGCACgttgtca                                                                                                              <  4:159964/51‑3 (MQ=37)
                                                                                                              ggAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATaa                          >  1:65807/1‑130 (MQ=255)
                                                                                                              ggAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATaa                          <  2:65807/130‑1 (MQ=255)
                                                                                                                aaaaaaTGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGcc                                                                                                         <  2:987987/49‑1 (MQ=255)
                                                                                                                  aaaaTGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTa   >  1:721210/1‑149 (MQ=255)
                                                                                                                     atgatgTCAATCAGATTGCTGACAACGTGCGCgtt                                                                                                                  >  2:1311430/1‑35 (MQ=255)
                                                                                                                     atgatgTCAATCAGATTGCTGACAACGTGCGCgtt                                                                                                                  <  1:1311430/35‑1 (MQ=255)
                                                                                                                          gTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGc                                 <  2:1365692/111‑1 (MQ=255)
                                                                                                                          gTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGc                                 >  1:1365692/1‑111 (MQ=255)
                                                                                                                             aaTCATATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAag                         >  2:844965/1‑116 (MQ=255)
                                                                                                                             aaTCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAag                         <  1:844965/116‑1 (MQ=255)
                                                                                                                             aaTCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTAc  <  2:1266397/139‑1 (MQ=255)
                                                                                                                               tCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTg                                       >  4:1137963/1‑100 (MQ=255)
                                                                                                                               tCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTg                                       <  3:1137963/100‑1 (MQ=255)
                                                                                                                                  gATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                               >  1:402476/1‑57 (MQ=33)
                                                                                                                                    ttGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGt                                     >  3:406017/1‑97 (MQ=255)
                                                                                                                                    ttGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGt                                     <  4:406017/97‑1 (MQ=255)
                                                                                                                                        gaaCAACGCGCACGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATaa                          <  3:2129/102‑1 (MQ=255)
                                                                                                                                        gaaCAACGCGCACGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATaa                          >  4:2129/3‑104 (MQ=255)
                                                                                                                                          aCAACGCGCACGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCaa                                                        <  4:1022191/72‑1 (MQ=18)
                                                                                                                                          aCAACGCGCACGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCaa                                                        >  3:1022191/1‑72 (MQ=18)
                                                                                                                                            aaCGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGa                             <  4:410572/97‑1 (MQ=255)
                                                                                                                                            aaCGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGa                             >  3:410572/1‑97 (MQ=255)
                                                                                                                                               gTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAagcg                       >  1:908517/1‑100 (MQ=255)
                                                                                                                                               gTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAagcg                       <  2:908517/100‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                    
TACCTGATAAATCAATCGCAAACAGCAAGAACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTAC  >  NC_000913/4079862‑4080125

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: